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Abstract 



Highly correlated electron materials are systems in which many new states of matter 
can emerge. A particular situation which favours the formation of exotic phases of the 
electron liquid in complex materials is that where a quantum critical point (QCP) is 
present the phase diagram. Neighbouring regions in parameter space reveal unusual 
physical properties, described as non-Fermi liquid behaviour. 

One of the important problems in quantum criticality is to find out how the Fermi 
surface (FS) of a material evolves near a QCP. The traditional method for studying 
the FS of materials is the de Haas van Alphen effect (dHvA). A quantum critical 
end point (QCEP) has been reported in the highly correlated metal SraPu^Cv, which 
is tuned using a magnetic field high enough to perform the dHvA experiment. It 
moreover features a new emergent phase in the vicinity of the QCEP, a nematic type 
of electron ordering. 

The subject of this thesis is the study of the FS of Sr 3 Ru 2 7 using the dHvA 
effect. Three aspects were explored. The first was the determination of the FS at 
fields both above and below that where the QCEP arises. The second was the search 
for quantum oscillations inside the nematic phase. The third was a reinvestigation of 
the behaviour of the quasiparticle effective masses near the FS. 

In collaboration with angle resolved photoemission spectroscopy experimental- 
ists, a complete robust model for the FS of Sr 3 Ru207 at zero fields was determined. 
Moreover, the new measurements of the quasiparticle masses revealed that no mass 
enhancements exist anywhere around the QCEP, in contradiction with previous spe- 
cific heat data and measurements of the A coefficient of the power law of the resistivity. 
Finally, we report dHvA oscillations inside the nematic phase, and the temperature 
dependence of their amplitude suggests strongly that the carriers consist of Landau 
quasiparticles. 
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Introduction 



Correlated electron materials are systems in which many exotic quantum states can 
emerge. These are phases of the electronic complex among many others that are pre- 
dicted by the fundamental expression of quantum mechanics, the Schrodinger equa- 
tion which, when one considers a specific arrangement of atoms, is difficult to solve. 
This equation gives rise to complexity from which these phases emerge. It is im- 
perative for the condensed matter physicist to find other ways to understand this 
phenomenon than starting from the Schrodinger equation, and one can, for instance, 
identify conditions that are favorable to their formation. Such phases of the electron 
liquid have important technological potential, for instance in the case of high tem- 
perature superconductivity (HTS), where room temperature superconductivity can 
potentially be achieved by understanding better how it is formed. But they are also 
of considerable intrinsic fundamental interest, and better knowledge regarding their 
nature could lead to the discovery of new types of emergent phenomena. 

One particular situation brings together the right conditions for the formation of 
many exotic phenomena, and is called quantum criticality. It corresponds to a region 
in parameter space, surrounding a so-called quantum critical point (QCP), which 
is dominated by quantum fluctuations and where in many cases superconductivity 
was found. Such a point may arise when the critical temperature of a second order 
phase transition is varied using a specific parameter, for example pressure or doping, 
and is made to reach zero temperature. Neighbouring regions in parameter space 
reveal unusual physical properties, described as non-Fermi liquid behaviour. Theory 
predicts several of these properties, notably an enhanced specific heat, which may 
also be expressed as an increase of the quasiparticle masses. 

The work of Jaccard et al. in 1992 and that of Mathur et al. in 1998 demonstrated 
the existence of superconductivity surrounding a QCP in heavy fermion materials, 
CeCu 2 Ge 2 PQ, CePd 2 Si 2 and Celn 3 [2j, where a Neel transition was suppressed to zero 
temperature using pressure. Figure [U left, presents the phase diagram of the first of 
these systems, the second being very similar. This research has become very influen- 
tial in the field of condensed matter, and many more systems have exhibited similar 
behaviour [3J H] . It has moreover been thought that quantum criticality could poten- 
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Figure 1: Le/i Work of Mathur et al. reporting superconductivity near a QCP in 
CePd 2 Si 2 p]. Right Typical phase diagram of high temperature superconductors 
(taken from Damascelli et al. [6]). 



tially be part of an explanation of the phenomenon of HTS, where a QCP would lie at 
the centre of the superconducting dome [5]. Figured! right, shows a typical phase dia- 
gram for high temperature superconductors, where superconductivity is found either 
by doping electrons or holes. 

An important problem in quantum criticality in metals is to find out how the Fermi 
surface (FS) evolves near a QCP. There are two main methods to study the FS of a 
material, which are by using angle resolved photoemission spectroscopy (ARPES) or 
the de Haas van Alphen effect (dHvA). It is however difficult to perform any of these 
near a QCP for the following reasons. ARPES is an experiment that is performed at 
ambient pressure and zero magnetic field, excluding these two parameters for tuning 
a QCP, and one is required to find a QCP that exists at zero pressure and field, which 



is not commonLU. dHvA is a probe of the FS that can be used more easily with tuning 
parameters, but requires that a magnetic field be applied, which may drive the system 
away from criticality, and very few QCPs tuned with magnetic field have been found 
at values high enough to observe the dHvA effect However, using such a material is 
an attractive option as it allows the analysis to be performed continuously near the 
QCP, an aspect not possible when using pressure or chemical doping. 

A quantum critical end point (QCEP) was reported in Sr 3 Ru207 [10], located at a 
field high enough to perform comfortably the dHvA experiment on both low and high 
field sides of the non-Fermi liquid region. Such a situation provides a unique oppor- 



1 Chemical doping can still be used with ARPES, and high temperature superconductors have 
been extensively studied in this regard. It is, however, not yet demonstrated whether a QCP exists 
in those materials. 

2 Most known QCPs in clean, undoped materials arise at low magnetic field values, for instance 
in YbRh 2 Si 2 , CeCoIn 5 and CeAuSb 2 0IS[9]. 
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tunity for the observation of the changes occuring in a FS that evolves from a normal 
Fermi liquid towards a non-Fermi liquid, and is the subject of this thesis. Moreover, 
this QCEP features, similarly to the cases reported by Mathur and co-workers, an 
unusual emergent phase of the electron liquid which is not superconductivity 
It was reported more recently to feature a nematic type of electron ordering, which 
produced an additional wave of interest in the community, motivating many theor- 
etical studies. Understanding the nature of this new phase also required information 
about the FS of S^Rx^C^, and part of our attention in this project was given to this 
issue. Finally, the dHvA effect also measures the quasiparticle effective mass, which 
in turn can be expressed as the electronic specific heat, and is expected to be highly 
enhanced near the QCEP. 

The goals of this project were very clear. We were determined to study and 
resolve three issues related to the FS of Sr 3 Ru20y. The first was to determine the 
three dimensional topology of the FS, which has been achieved in collaboration with 
ARPES experts. The second was to study the FS of Sr3Ru207inside the new nematic 
phase, and was successfully performed. The last was to revisit and analyse in detail 
the evolution of the quasiparticle properties near the QCEP, a study that revealed 
an unexpected result. 

This thesis is divided in four chapters. We present in the first chapter a review 
of the current knowledge surrounding Sr 3 Ru 2 7 , its QCEP and nematic phase, but 
also the standard theory of dHvA. Emphasis is given to all theoretical aspects of 
relevance for the correct analysis and interpretation of the measured dHvA data, as 
well as to the previous research that makes up what is known about the electronic 
structure of S^R^Oy. The second chapter introduces all the experimental methods 
that were used in this project along with the numerical methods for the analysis of 
dHvA data that were employed. The dHvA experiment has required the isolation of 
ultra-pure samples of Sr 3 Ru 2 7 , and this procedure and its results are described in 
detail in this chapter as well. Chapter three presents all the dHvA data measured by 
the author, raw and processed, with a basic analysis and interpretation. A detailed 
analysis follows in chapter four, where additional information is incorporated from 
other experiments, notably ARPES. A complete model for the FS is introduced, 
and the consequence of the nematic phase and quasiparticle mass measurements are 
discussed. The properties of the zero, low and high field FS are described, as well 
as those of the nematic phase region. Finally, we conclude with a brief review of the 
information gained in this project. 



Chapter 1 

Scientific Background 



1.1 Strongly correlated metals 

We introduce in this section the physical properties of strongly correlated electron 
systems, relevant to this work. Electrons in a correlated system are usually well 
described by Landau's Fermi liquid theory, a framework which enables condensed 
matter physicists to understand and predict the properties of most metals. However, 
we are interested in this project by a case which features characteristics that diverge 
from the normal properties predicted by the Fermi liquid theory, and we are therefore 
required to review the theory of non-Fermi liquids. This introduction to the subject 
is very brief, as this work is mainly of experimental nature, but we provide references 
for further reading. 

1.1.1 The Fermi Liquid 

Electrons in simple metals are usually well described by the Sommerfeld theory of 
the electron gas. Electron states are represented by an orthogonal set of plane waves. 
Consideration of the effects of an underlying periodic potential on the electron gas 
leads to dispersion relations, featuring band gaps, that can be more or less complex, 
and to a wide variety of different Fermi surface shapes. However, the basic physics 
of the Sommerfeld model remain relatively simple: the Pauli exclusion principle is 
respected and at zero temperature, electrons fill the lowest states in k space up to 
a surface of constant energy, the Fermi energy. For instance, in simple metals like 
the alkali elements, the Fermi surface is close to a perfect sphere. In this case, the 
dispersion relation is close to quadratic, like that of a gas of free electrons. In other 
cases, the Fermi surface can be very complex. 

In this framework a basic physical quantity is usually derived from the dispersion 
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relation e(k) of an electron gas, the velocity [12], 



v(k) 



1 



V k e(k). 



This expression leads to the concept of the effective mass m* of the electron, 



The Sommerfeld model has been extremely successful in the study of condensed 
matter, but it neglects a fundamental property of the electron gas, the Coulomb 
interaction and hence electronic correlations. Consequently, the question arises as to 
why the model works so well. In the 1950s, Landau constructed a theory that became 
the standard model for the physics of metals, called Landau's Fermi liquid theory. It is 
based on a critical assumption, now called adiabatic continuity, that when turning on 
the interaction between electrons in a non-interacting gas, the low energy excitations 
evolve continuously from the Fermi gas to the Fermi liquid, preserving their energy 
ordering and a one to one correspondence with those of the non-interacting gas. 
Consequently, in this model the Fermi surface picture of the Sommerfeld model is 
still correct. The wave functions of the new states are different, but they can be 
approximated by orthogonal linear combinations of plane waves, and the fundamental 
unit acting as an electron is called a quasiparticle. 

It is surprising that even with strong interactions, the excitations remain relatively 
long lived. One of the reasons why such a scheme works is because of the screening 
of the electrons. The coulomb potential has a very long range, but in an electron 
gas, any fluctuation of charge in space should in equilibrium be compensated by a 
variation of the density of electrons. The interaction is actually of very short range. 
An electron moving will affect many others, but the global movement of the charge will 
be well defined. The wave functions of the quasiparticles are not exact eigenstates of 
the Hamiltonian featuring the interactions and, even at zero temperature, excitations 
above the Fermi surface will be scattered. It can be shown [13] that the scattering 
rate of an excitation due to electron-electron interaction depends on its energy as 



which is much smaller than the excitation energy e itself, when e is sufficiently small. 
A quasiparticle is better defined the closer it is to the Fermi surface. 

For isotropic Fermi surfaces, Fermi liquid theory makes the following predictions 
for the resistivity due to interactions, the specific heat and the magnetic susceptibility 



hk 




v = — 



m* 



1 
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of a one band electron liquid at low temperatures: 



*~ = a °^3| t2 ' (l2) 

^ = ^IT, (1.3) 



3ft 3 
1 

l + Fq n 2 h 



x = , '/ Bk r\ (i.4) 



where i 7 ^ is one of the so called Landau F parameters, A is a dimensionless constant 
and n is the density of free carriers [T21 [T3J [H]. These quantities all involve the 
effective mass m*, into which was put part of the effect of interactions. The theory 
predicts that it should be a constant, determined using specific heat or, as we will see 
later on, de Haas van Alphen measurements. 



1.1.2 Non- Fermi liquids and quantum criticality 

In 1986, high temperature superconductivity was discovered by Bednorz and Miiller 
in layered materials involving copper oxide [T3] . It was found to be a type of uncon- 
ventional superconductivity that could not be understood with previously developed 
BCS (Bardeen-Cooper-Schrieffer) theory which applies to simple metals. The metal- 
lic state in transition metal oxides is very anisotropic and unusual and it was naturally 
thought that the superconducting state of a material would be understood only with 
a correct description of its normal metallic state. It was found that these metallic 
states are not very well described by Fermi liquid theory. For example, power laws 
for the temperature dependence of the resistivity were observed with exponents lower 
than 2, contrary to the Fermi liquid theory prediction for electron-electron interaction 
in highly correlated metals [13J. Examples of non- Fermi liquid behaviour include the 
fractional quantum Hall effect in 2D metallic systems, the Luttinger Liquid in one 
dimensional metals, the Kondo effect and quantum criticality, which is of interest 
here. 

A way towards non- Fermi liquid behaviour is to find a state of the electron liquid 
where the scattering rate is so great that the quasiparticles cannot exist for long 
enough to be well defined. One of the models proposed is the proximity to a QCP 
[TBlfTT] . or similarly to a QCEP [1O_0. Proximity to such a point in the phase diagram 
is dominated by quantum fluctuations. A quantum critical point arises when a second 
order phase transition approaches absolute zero in any phase diagram. For example, if 
the critical temperature of a second order transition can be varied using a non thermal 
control parameter like pressure or chemical doping such that it is depressed to absolute 

1 Described in section H .2 .41 
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Figure 1.1: Generic phase diagram for a heavy fermion system exhibiting a quantum 
critical point (image taken from the work of Coleman et al. [IB])- The Neel temper- 
ature is tuned to zero Kelvin using the control parameter p, to a QCP at P c . The 
renormalized Fermi temperature Tp is shown under which Fermi liquid behaviour is 
observed. The region between the line of transition and Tp exhibits non-Fermi liquid 
behaviour. 

zero, the end point constitutes a QCP. As the entropy of a system approaches zero 
when lowering the temperature to absolute zero, the transition then occurs between 
two perfectly ordered states and is dominated by quantum fluctuations. Figure 11.11 
shows the typical phase diagram for heavy fermion quantum critical systems, where 
the transition in question is the Neel transition, at which the system moves from an 
antiferromagnetic to a paramagnetic metal, and is tuned with control parameter p. 

Discovery of a possible close link between QCPs and high T c superconductivity 
was made in a variety of systems, for example, by the observation of a linear res- 
istivity at all temperatures outside of the superconducting phase in optimally doped 
LaSr 2 _ 1 .Cu 1 .04 [19], but not at any other doping values [20]. Conversely, in a few 
heavy fermion systems, unconventional superconductivity was found using pressure 
to tune a transition from antiferromagnetism to paramagnetism towards K, cover- 
ing the region where a QCP was expected to exist J2J. In this view, the proximity 
to a QCP and the domination of quantum fluctuations are thought to provide a 
mechanism for types of superconductivity not mediated by phonons. 
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Regions near a quantum critical point in a phase diagram are characterised by un- 
usual metallic properties. A metallic system containing a QCP to which the itinerant 
electrons have some coupling usually exhibits the following properties [T3] HH] : 

1. Fermi liquid behaviour away from the QCP at zero temperature, 

2. Close to linear temperature dependence of the resistivity, 

p(T)=p + AT a , (1.5) 

with a near 1, 

3. Divergent specific heat as a function of temperature at the QCP, 

-jr- =7oM-y): ( L6 ) 
which suggests, according to equation II .3[ a diverging effective mass m* . 



1.2 The strontium ruthenate oxide Sr3Ru207 

We present in this section the properties of the material of study, Sr 3 Ru 2 07, that we 
feel are relevant to this thesis. Numerous papers, both experimental and theoretical, 
have been written on the subject, and our aim is not to give a complete review on 
the subject. We rather insist on most aspects related to its crystalline and electronic 
structure, phase diagram and quantum critical properties. We therefore begin by list- 
ing the most important general physical characteristics of Sr3Ru207. This is followed 
by a description of its chemical structure, in order for the reader order become aware 
of the dimensionality of its electron liquid and the space group symmetry. We then 
present the electronic structure of Sr 3 Ru207, where using the space group symmetry, 
we review the shape and size of its Brillouin zone. We moreover construct an ap- 
proximate Fermi surface starting from arguments related to the various <i-bands that 
the quasiparticles populate. Finally, we present a review of the discoveries that were 
made in the past related to the phase diagram of the electron liquid in Sr 3 Ru207, 
and the existence of a QCEP. 

1.2.1 General physical properties 

Sr 3 Ru207 is a paramagnetic metal, thought to be close to a ferromagnetic instability. 
It features a peak in the susceptibility, as a function of temperature, near 16 K (see 
fig. 1 of the work of Ikeda and co-workers [21], or equivalently, figure [27161 section [2.4.51 
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n = 1n = 2n = 3 n = oo 

Figure 1.2: Schematic representation of the Ruddlesden-Popper series Sr„, + iRu n 03 n+ i 
with increasing n, showing from left to right Sr 2 Ru04, Sr 3 Ru207, Sr 4 Ru30io and 
SrRuOs (n = oo). The octahedra represent the Ru ions with six O at their corners, 
and the spheres correspond to the Sr ions. The crystallographic axes are given in the 
top right corner. 

of this thesis), and Curie- Weiss behaviour above 200 K (Qw — -39 K). It can be pushed 
into the ferromagnetic state using pressure, where under lGPa, magnetic ordering 
appears around 70 K [2T]. It possesses a highly two dimensional character, and the 
resistivity anisotropy is of a factor 300 at 0.3 K. It is nevertheless metallic, with strong 
e — e correlations, as the low temperature resistivity follows a T 2 behaviour. The 
specific heat coefficient 7 was found to be large, 110 mJ/mol Ru K 2 , much higher 
than in some of the other ruthenates, where in Sr2Ru04, 7 = 38 mJ/mol Ru K 2 
and in SrRuC>3, 7 = 30 mJ/mol Ru K 2 , but similar to that of Sr 4 Ru30io, where 
7 = 109 mJ/mol Ru K 2 . The Wilson ratio of Sr 3 Ru 2 7 is also high, R w = 10, due to 
strong ferromagnetic correlations. Consequently, this material is viewed as a strongly 
correlated metal on the verge of ferromagnetism. 
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1.2.2 Crystalline structure 

Sr3Rii207 is member of a family which has been extensively studied, the strontium 
ruthenate oxide family. Its crystallographic structure has been identified by Ruddles- 
den and Popper [22] as a perovskite series of the form Sr n+1 Ru n 3n+ i. The effective 
dimensionality of the electron liquid changes as n is varied: it defines the number 
of consecutive planes of ruthenium oxide, leading to more or less anisotropy of the 
metallic state. Figure 11.21 presents the crystal structure of a few members of this 
family schematically. The n = 1 member Sr 2 Ru0 4 , famous for its unconventional 
superconductivity, probably of the p-wave type, has been extensively described in a 
review by A. P. Mackenzie and Y. Maeno [23], and is isostructural to high T c su- 
perconductor LaSr 2 _ x Cu x 04. It is the most strongly two dimensional member of 
the family, n = oo SrRu0 3 is an itinerant three dimensional ferromagnet with a 
Curie temperature of around 160 K [2H [25j [26]. The n = 2 compound Sr 3 Ru207 is 
an exchange-enhanced paramagnet exhibiting a metamagnetic transition [27]. It has 
one bilayer per formula unit and consequently a dimensionality intermediate between 
that of Sr2RuC>4 and SrRu0 3 . Finally, the members with n > 2 possess physical 
properties similar to those of SrRu0 3 , with a dimensionality that varies with n. Of 
interest here, Sr 4 Ru30io, the n = 3 member is a ferromagnet with Curie temperature 
of 105 K [281 EH]- It is important to note that stacking faults can transform part of 
a crystal of one member of the ruthenate family into another member. For instance, 
suppressing RuO planes from Sr 3 Ru 2 07 leads to inclusions of Sr 2 Ru0 4 . In a similar 
way, duplicating planes in Sr 3 Ru207 will produce Sr 4 Ru30io. One can also deduce 
from the similarity between the lattice cell of high n and SrRu0 3 that such systems 
should have similar properties. 

Strontium ruthenate crystals are usually grown using the floating zone technique. 
In the case of Sr 3 Ru207 the crystals used in this project were produced by R. S. 
Perry using the method designed by Perry and Maeno [30], and has yielded crystals 
of exceptionally high purity and low disorder. It is well known that each member of 
the Ruddlesden-Popper family, during growth, requires specific stochiometric ratios 
of Ru and Sr. A growth performed with intermediate ratios results in a mix of 
different members (see, for instance, the work of Fittipaldi and co-workers [3"T]). 
Moreover, if during growth the stochiometry inside the floating zone region drifts 
slightly, stacking faults appear in order to correct the relative amounts of Ru to Sr, 
producing inclusions of other compounds. Consequently, the stochiometry needs to 
be very accurate in order to grow crystals of high phase purity. We will, in this 
project, present a method for identifying volume fractions of the various phases of 
the ruthenate family in samples (section 12.41) . 

Sr 2 Ru0 4 possesses a body-centred tetragonal crystal structure, of space group 
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Figure 1.3: Schematic representation of the octahedron rotation in SraRx^CV Left 
In-plane representation of the octahedra before rotation, top, and after, bottom. The 
basal unit cell is shown in blue. Middle Undistorted lattice unit cells. Right Distorted 
unit cell. Yellow octahedra are not rotated, the green are rotated clockwise in the 
afr-plane, and the red counter-clockwise. The respective crystallographic axes are 
represented, the ones for the distorted lattice rotated by 45° around the c-axis. 

symmetry IAmmm [23J. Consequently, it has fourfold rotational and inversion sym- 
metry. Early work on Sr 3 Ru207 reported it to possess the same space group [32] . 
but it was later found using neutron diffraction that this symmetry is broken by a 
7° rotation of the RuO octahedra [33], [31]. The resulting space group is Bbcb, corres- 
ponding to single face centred orthorhombic structure with a square base of sides \/2 
larger than those of the undistorted crystal. In this new space group, one no longer 
has fourfold rotation, even though the lattice parameters a and b are still equal and 
impossible to distinguish in a Laue experiment. Figure ITT31 presents schematically the 
rotation of the octahedra. On the left, one can see its effect on the basal plane, where 
blue lines represent the base of the unit cell. The middle part of figure 11.31 shows the 
unit cell before rotation, with space group IAmmm. The right part shows the effect 
of the distortion where each neighbouring octahedron is rotated in opposite direc- 
tion, green shapes indicate those rotated clockwise and red ones counter-clockwise. 
The primitive cell contains twice the number of atoms compared to the undistorted 
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I4mmm 



Bbcb 




Figure 1.4: Schematic representation of the Brillouin zone for IAmmm , left, and Bbcb, 
right, for the case where the lattice parameters are such that c > a = b (inspired from 
Bergemann [35J). 



lattice 2 



1.2.3 Electronic structure 

The rotation of the octahedra in Sr 3 Ru 2 07 has a significant effect on the electronic 
system. The first and obvious modification is the shape of the BZ, due to the change 
in lattice type. The BZ possesses the same symmetries as the crystalline structure, 
plus inversion through the centre. Consequently, the new BZ should possess a lower 



symmetry than the initial one. Figure 11.41 presents BZs for space groups IAmmm, 
left, and Bbcb, right, when the c lattice parameter is longer than both a and b, which 
are equals For IAmmm, we have the same BZ as for Sr 2 Ru04, which is fourfold 
rotational symmetric. In the case of Bbcb, the symmetry is not fourfold and the 
k x and k y axes are not equivalent, even though the in-plane lattice parameters are 
equal. Finally, for Sr 3 Ru 2 07, the base of the BZ has an in-plane area that is twice 
that of Sr 2 Ru0 4 . Table [TTT1 shows these areas for Sr 2 Ru0 4 and Sr 3 Ru 2 7 , in units of 
fc-space and dHvA frequencies (defined later in section 11.3. ip . which will be useful in 
the discussion, chapter HJ 

The FS of Sr 3 Ru 2 07 can be constructed qualitatively from hand- waving argu- 
ments, in the spirit of Bergemann et al. [35] (see figure 113]) . This does not involve 
any band structure calculation, but only the simple rule that bands cannot cross one 
another, and hybridise instead. We start from the situation in Sr 2 Ru04, where three 
bands are present, a, (3 and 7. These originate from the hybridisation of the various 



2 The undistorted primitive cell contains two octahedra, while the distorted one has four. 
3 Note that for Sr 2 Ru04 and Sr 3 Ru 2 07, the shape of the BZ will be flatter and wider than appears 
in figure \TM 
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BZ Area BZ Area 
A- 2 kT 


Sr 2 Ru04 


2.61 27.74 


Sr 3 Ru 2 07 


1.31 13.69 



Table 1.1: In-plane area of the BZ of Sr 2 Ru04 and Sr 3 Ru 2 07, in units of fc-space and 
dHvA frequencies. 

d bands of the Ru atoms, the d xz , d yz and d xy , which, due to crystal field splitting, 
are the ones crossing the Fermi level, all the others lying at higher energy values. 
The d xz and d yz give rise to almost no dispersion either in the z and, respectively, 
y and x directions, but rather to quasi one-dimensional hopping. Consequently, the 
Fermi surface corresponding to these bands should be close to planar in the k x k z and 
k y k z direction. These will cross in certain regions of the BZ (see figure [T31 a.), where 
hybridisation gaps will appear, and the sheets will reconnect into closed surfaces. The 
d xy orbital, however, allows hopping in all directions, and the corresponding Fermi 
surface is close to a perfect circle in the k x k y plane. The resulting FS with all bands is 
shown in figure 11.51 b. This construction is consistent the fourfold rotation symmetry 
of the IAmmm space group. 

In Sr 3 Ru207, one expects each band to duplicate, due to bilayer splitting 0. Con- 
sequently, the result is slightly different than for Sr 2 Ru04. One starts from six sur- 
faces (figure fl~5| a), four that originate from the d xz , d yz orbitals, which reconnect 
at the points of crossing, and one should obtain a result similar to that shown in 
figure 11.51 d. This is what the FS could be without the \/2 x a/2 reconstruction due 
to the octahedral rotation. But since the BZ reconstructs into a square twice smal- 
ler, back-folding of the bands occurs, shown in figure 11.51 e - 111 this last plot, many 
band crossings appear, and the way by which the surfaces reconnect is very complex. 
However, one can obtain hints from recent ARPES measurements Five orbits 
are expected, shown in figure 11.51 f, which take the form of square and cross shaped 
hole pockets in the centre, originating from the d xz and d yz orbitals, two lens shaped 
electron pockets at the M point, and a small pocket near the X point. The complete 
result from ARPES is more complex and will be discussed in section 14.11 Note that 
the result possesses fourfold rotation symmetry although the space group Bbcb does 
not. This is due to the fact that within a bilayer, the structure is fourfold symmetric, 
and only the stacking of the layers is not, seen in the right side of figure 11.31 Since 
Sr 3 Ru207 is quasi two-dimensional, we expect the FS to be very close to fourfold 
symmetric. 

4 After distortion, Sr3Ru2 0y possesses four instead of two non-equivalent Ru atoms, bringing the 
number of d bands crossing the Fermi level from three to six. 
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Figure 1.5: Schematic construction of the Fermi surface of Sr2Ru04, a and b, and of 
Sr 3 Ru207, c, d, e and /, where in the last two, the \/2 x V2 reconstruction of the 
BZ accompanied by back-folding was included, and the new BZ is represented by the 
grey diamond. All axes were scaled by the parameters of the undistorted lattice. 
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Figure 1.6: Le/t Fermi surface structure as calculated by Singh and Mazin 
represented in the \[2 x \[2 reconstructed zone with its centre at the corners of the 
figure. The left part of the figure corresponds to a cut of the FS at k z = 0, while the 
right part shows the FS at the top (or bottom) of the BZ. Right Three dimensional 
phase diagram for SrsRu207, with axes of temperature, magnetic field and magnetic 
field angle. The green surface corresponds to the metamagnetic transition, and the 
black line its critical end point. 



Singh and Mazin calculated the band structure using general potential linearized 
augmented plane wave method and, taking into account the distortion, suggested the 
Fermi surface shown in figure [T^l [37] . where the FS is plotted with the X point in the 
centre. In such a representation, the left side is located within a cut at k z = 0, while 
due to stacking of the BZ (see figure U~M . the right side corresponds to a cut at the 
top of the zone. The calculation features, in contrast to the ideal tetragonal case, lots 
of small electron and hole pockets, and a FS that breaks slightly the fourfold rotation 
symmetry. One can see that identifying all the different FS sheets can potentially be 
a difficult task without ARPES data. 



1.2.4 Metamagnetism and quantum criticality 

Metamagnetism, for a paramagnetic metal, is defined as a sudden superlinear rise in 
magnetisation as a function of applied magnetic field. Such phenomena fall into two 
categories, the first corresponding to an antiferromagnetic transition, and the second 
to a rise in magnetisation in an itinerant paramagnet, the latter being the one of 
interest here. In Sr 3 Ru207, it has the characteristics of a first order phase transition, 
as was inferred from real and imaginary AC magnetic susceptibility, measured by 
Grigera and co-workers [38]. In these experiments, the real part of the susceptibility 
(x r ) showed a peak at a critical field, resulting from a jump in magnetisation, while the 
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Figure 1.7: Left Exponent of the temperature dependence of the resistivity (eq. 
II. 5p as a function of magnetic field and temperature, obtained using the equation 
a = d\n(p — p res )/dlnT. Right Field dependence of the A coefficient of eq. 11.51 
obtained by calculating (p(T) — p res )/T 2 , and that of the residual resistivity p res . 
Both plots were reproduced from the work of Grigera et al. [10J. 

imaginary part (x") also showed a peak, indicating dissipation and signalling the first 
order nature of the transition. A field angle study of the complex AC susceptibility 
revealed a sheet of first order transitions in the 3D phase space defined by magnetic 
field, angle and temperature (H,9,T), terminated by a line of critical points, at which 
the peak in x" disappeared (shown in figure fL6l right). At higher temperatures, the 
peak in x" was n °t detected. This line was seen to approach absolute zero at a field 
of 7.85 T and an angle close to 90° with respect to the ab plane. 

This field angle analysis of the susceptibility in Sr 3 Ru207 led to the introduction 
of a new idea, a quantum critical end point (QCEP) arising where the line of critical 
end points crosses the T = plane. As described in section ll.l.2[ a QCP usually 
arises when a continuous (second order) phase transition reaches K, but a first order 
phase transition does not lead to a QCP. However, the critical end point terminat- 
ing a first order transition line in a phase diagram exhibits all the properties of a 
continuous transition except, in this case, spontaneous symmetry breaking. Grigera 
et.al. suspected a line of critical end points reaching K to have quantum critical 
properties, that is, to form a QCEP with field angle as a tuning parameter. 

A careful study of both the specific heat and resistivity in the (H, T) phase space 
revealed strong evidence of the existence of such a critical point and of non-Fermi 
liquid behaviour (TUJ U7\ • Figure 11.71 shows the exponent a of the T dependence of 
the resistivity (eq. I1.5P for a field aligned in the direction of the c-axis, using the 
equation 

_ dlnjp - p res ) 
a ~ d\nT ' 
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Near 8 T, linear T dependence was seen for all temperatures, but away from this 
field, T 2 behaviour was recovered, suggesting that a Fermi liquid existed on both 
sides of the QCEP. As the resistivity became linear near at 8 T, the parameter A of 
equation II .5\ obtained by calculating (p(T) — p res )/T 2 using data measured between 
0.2 and 0.9 K, was also seen to increase sharply. This suggested an enhancement of 
the quasiparticle mass of at least more than one band, as can be seen from equation 

off 

The specific heat of Sr 3 Ru 2 7 was also measured, and reproduced one of the 
predictions of quantum critical theory. Figure 11.81 shows curves of the electronic 
specific heat C e i/T as a function of temperature. Logarithmic divergence as a function 
of temperature, in accord with eq ll.6[ was seen near the QCEP, which again indicated 
a divergence of the quasiparticle effective mass m* (see eq. I1.3p . It became desirable 
to perform an independent measurement of the quasiparticle mass, and quantum 
oscillations experiments were performed by Perry et.al. [39], an d later by Borzi and 
co-workers near the QCEP |40j, where both showed a slight transformation of the FS 
across the metamagnetic transition and the latter claimed a divergence of the mass 
of two FS sheets near the QCEP. The work by Borzi et.al. will be described in more 
details in section 11.3.61 

The system under study met all the requirements (see the list, section I1.1.2P 
for quantum critical behaviour, as was demonstrated by Grigera et al. [TO]. However, 
anomalies were found very near the metamagnetic transition, between fields of 7.8 and 
7.9 T, where a decrease in A was observed and a resistivity temperature exponent of 3 
was measured between fields of 7.82 and 7.86 T. As new samples were later generated 
with lower residual resistivities, even more features were discovered near the QCEP, 
described in the next section. 

Finally, we add a brief note about the nature of the magnetic fluctuations present 
in Sr 3 Ru 2 7 , determined with inelastic neutron scattering by Capogna et al. [31] . 
Although at high temperature, the fluctuations are of ferromagnetic nature, as one 
expects, incommensurate antiferromagnetic fluctuations develop at zero fields as the 
temperature is lowered. Effectively, a resonance at q = was observed at a tem- 
perature of 150 K at an energy of 3.1 meV, which vanishes at around 15 K, and 
corresponds to ferromagnetism. Alongside this excitation were observed resonances 
at q = [±0.25 0] and [±0.09 0], which correspond to antiferromagnetism, at tem- 
peratures of around 15 K and below. The system is thought to exhibit a competition 
between the two types of magnetic ground states. 

5 Note that A was obtained from fitting the data only in regions where it follows a T 2 dependence, 
a region that shrinks in length as one comes closer to the metamagnetic transition. 
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Figure 1.8: Left Electronic specific heat as a function of temperature [27]. Log- 
arithmic divergence is observed at 7.7 T. Right New phase reported in Sr3Ru207 
surrounding the QCEP, taken from Grigera and co-workers [llj. Symbols represent 
features in data taken using various experimental techniques: coloured triangles are 
the centre of the peaks in the AC susceptibility, white triangles are the loci of peaks 
in the thermal expansion, orange circles are taken from the DC magnetisation, yellow 
diamonds were obtained from dp/dH, green diamonds are from d 2 p/dT 2 and the blue 
circles were extracted from maxima in the magnetostriction. The red arrows indicate 
the critical points at the end of each of the two lines of first-order transitions. 



1.2.5 Disorder sensitive phase formation 



As a new generation of samples were grown [SU], new experiments were performed, 
where the metamagnetic transition was revealed to be double [39]. From various phys- 
ical properties^ of Sr3Ru207 near 8T, strong evidence was found indicating a new 
ordered phase of the electron liquid surrounding the QCEP [11]. This phenomenon 
was reminiscent of famous cases where superconductivity was found surrounding a 



QCP, for instance in the work on CePd2Si2 and Celn3[2j. Figure fL8| right, shows 
the phase diagram of Sr 3 Ru 2 7 near 8 T, at low temperatures. In this plot, one can 
see the bounded phase in the centre, and two first order transition lines extending to 
higher temperatures away from the phase region. The nature of the new phase was 
not known but Fermi liquid behaviour was confirmed both to higher and lower mag- 
netic fields. Note that the sides of the phase feature much a stronger thermodynamic 
signature than the so-called "roof'L. 

The boundaries of the phase were studied in three dimensional space, and it was 
found to extend in angle to approximately 30° from c-axis, shown in figure [LT9"| left. 
In such a new phase, it was expected that symmetry should be broken, and that 



6 These included resistivity, magnetic susceptibility, magnetisation, thermal expansion and mag- 
netostriction 

7 This subject has been extensively studied by Rost et al., in preparation. 
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Figure 1.9: Left Three dimensional representation of the various phase transitions 
delimiting the new phase (see text), from Green et al. [44] . Centre and Right Nematic 
nature of the new phase, taken from Borzi et al. 



domains could be formed. Magnetoresistance experiments at different magnetic field 
angles were performed, comparing situations where it was aligned or perpendicular 



to the current, which showed anisotropy, shown in figure 11.9] centre and right. It 
was concluded that the electron liquid possessed nematic properties [32]. The ori- 
ginal article suggested that this could correspond to a deformation of the FS which 
produces domains that impede the transport of current in a preferential direction 
that depends on the orientation of the magnetic field. Other studies claim that it 
may originate from Condon domains [33], or from the appearance of a spontaneous 
transverse magnetisatio However, the first of these two was ruled out by a careful 
study of the the system with various crystal geometries, which showed no variation 
of the critical fields. We found, during this project using dHvA, that the system still 
possesses metallic properties, by the observation of quantum oscillations inside the 
phase (see section [375]) . 

One of the very interesting properties of Sr 3 Ru 2 7 is the dependence of its prop- 
erties on disorder. This has motivated us in performing more systematic growth and 
characterisation of crystals, in the hope of finding samples of outstanding quality. 
This has been one of the main aspects of this project, and is described in section [2741 



1.3 The de Haas van Alphen Effect 

The dHvA experiment is the traditional method for studying the FS of a material. 
It is normally used for determining the volume and shape of the FS, as well as 
quasiparticle masses. Accessing a quantum critical point requires a control parameter 

8 Work from Green et al., to be published. 
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that is usually not the magnetic field, and when it is, low field values have usually 
been required (see for example |7|). Applying a field thus tunes the system away from 
quantum criticality, so that the dHvA experiment is not possible near a QCP. It is 
not the case for Sr 3 Ru207, since the QCEP is field tuned between 5 and 8 T with 
field angle 9. It is the ideal system for exploring a FS near a QCP. Since the dHvA 
experiment yields the quasiparticle mass for each branch of the FS, it is possible to 
obtain them as a function of both field and field angle, but moreover, it can show 
whether there are changes in the FS topology as well. This section thus reviews the 
basics of the dHvA experiment, along with methods for calculating the quasiparticle 
masses as a function of magnetic field. For more details, the reader is referred to the 
main reference on dHvA, the book by Shoenberg |45j . and a review on Sr 2 Ru04 by 
Bergemann et al. [35] . 

This section is divided into six parts. We first derive the basic equations and ex- 
plain the origin of the effect at zero temperature in a perfect lattice. Next, we include 
the effects of finite temperature and quasiparticle mean free path. Then, a short 
digression is given about the effect of the spin. We discuss how the two dimensional 
nature of a material affects dHvA and finally, we review the first measurements that 
were performed on SrsR^Cv prior to this project. 

1.3.1 Oscillations at zero temperature 

We present here a semi-classical derivation of the quantum oscillations at zero tem- 
perature. The source of the quantum oscillations is the Landau quantisation of the 
orbital movement of electrons in a magnetic field. When a magnetic field is applied to 
a free electron gas in a box, in the z direction, electrons undergo an orbital or helical 
movement that corresponds to quantised orbits in k space, the k z quantum number 
remaining unchanged, with energies 



where u c = eH/m e is the cyclotron frequency. In such a system, without field, 
electrons would, at T = 0, take the lowest k values and form a sphere in k space, the 
highest energy being the Fermi energy ep, whereas with the field turned on, they will 
fill so-called Landau tubes, concentric hollow cylinders up to the surface of the same 
sphere. Effectively, for high quantum numbers n, the degeneracy of the Landau levels 
is such that within a fc-space area delimited by a contour of constant energy, without 
magnetic field, the number of states equals approximately the same as the sum of 
the number of states available in all the Landau levels included within that area for 
a specific field, and that degeneracy is linear in field. Consequently, the Fermi level 
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lies at the same energy value with or without a magnetic fielcj^. 

More often than not, the electronic system forms a Fermi surface that is not 
isotropic. The orbital movement in k is not circular but follows paths along surfaces 
of constant energy in k space. In a magnetic field, the rate of change of momentum 
of the quasiparticles is 



where 



V t €(k). 



v is always normal to a surface of constant energy, and the rate of change of mo- 
mentum is in the direction of constant energy, in a plane perpendicular to H. It can 
be integrated in time to give 



k-kn 



h 



r — r ) x H. 



;i.7) 



The Bohr Sommerfeld quantisation rule for periodic motion states that, for canonic- 
ally conjugate operators p and q, 



p • dq = 2nh{n + 7), 



where n is the orbital quantum number and 7 is a constant, p is the canonical 
momentum Kk — eA, A being the vector potential and q is the component of the 
position operator r in the plane perpendicular to H, that we will denote r'. Using 
(11. 7p for the first term and Stokes' theorem for the second, we obtain: 



/p.*, = e /[«k-«A].*, 



H • (V x dr') 



H • dS 



where S is the surface enclosed by the orbit in real space. In this expression, r should 
be decomposed into its components parallel and perpendicular to H, and only the 
first will contribute. The second integral is equal to |H| times the area S, while the 
first is twice that: 

h 2HS -HS = HS 



so that 



27r(n + 7) — 

e 



S(n,H) = 27r(n + 7 ) 



A 

eH' 



9 High quantum numbers n arise when the Fermi energy is high compared to nhuj c , which is true 
in most metals; exceptions to the rule are materials with very small Fermi surfaces, for instance 
doped semiconductor layers. 
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This means that it is the area in real space that is quantised, and it is also inversely 
proportional to the field H. We are interested in the area in k space, though, and 
from eq. fll.7jl we note that 



k- k = — — 

h 



fo- 



ri. 



It follows from the scaling factor, that the area in k space is proportional to that 
in real space through 

feH\ 



so that 



S k (n,H)=l—j S(n,H), 



eH 

S k (n,H) = 27r(n + 1 ) — , 

a 



which is proportional to the field H. Consequently, for any shape of orbital motion, 
the Landau tubes have a cross-sectional shape defined by constant energy paths in 
planes perpendicular to the magnetic field. If the value of H increases or decreases, 
then the tubes, respectively, grow or shrink in size. When a tube crosses a region of 
the Fermi surface parallel to H (an extremal region of the Fermi surface with respect 
to the field direction), a sudden change in the DOS at the Fermi surface arises, with 
a period that is a constant of the inverse field X = 1/H. This can be more easily 
understood by looking at the area of the nth tube crossing the Fermi surface at a 
field Hi, and that of the next, n — 1 at H 2 , as H increases: 



S k {n, Hi) 
S k (n - l,H 2 ) 



, \27re TT 
{n + ^—Hi, 



n 



1+7)^2- 



h 



These areas are equal, and by isolating n, we calculate 

S k h s k H 



2-neHi 



7 



2neH 2 



7 + 1, 



The period in inverse field Tyn is a constant of the system, and corresponds to a 
dHvA frequency F: 



1 

Hi 



Ti/h — — — — 



1 



1 

H 2 
-i -l 



K H> 



2ne 
hSk 
hSk_ 
2vre' 



;i.9) 



The frequency F is proportional to the cross-sectional area S k of the Fermi surface. 
The critical consequence of this phenomenon is that all physical quantities that 
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depend on the DOS at the Fermi level exhibit oscillations as the magnetic field is 
swept o, of the form 

p 

where p stands for harmonics of the fundamental frequency F, due to the fact that the 
oscillations are not sinusoidal. The name de Haas- van Alphen was given to oscillations 
in the magnetisation, by the name of their discoverers [IB] , while Shubnikov-de Haas 
(SdH) was given to those in the resistivity |47j . and these experiments followed a 
prediction by Landau in his work on diamagnetism [48J . Applying a Fourier transform 
to quantum oscillation data as a function of X = 1/B, one extracts the area of all 
extremal orbits that exist in a field direction, and by rotating the sample, one can in 
principle construct the full 3D shape of the FS. In 2D systems, frequencies of different 
extremal orbits may be very close together and give rise to beat patterns, and this is 
discussed in section 11.3.41 

For later use, we rewrite the energy e n of the quasiparticle as a function of the 
area of the nth Landau level and its k z wave vector: 

h 2 ^ h 2 k 2 

en = A^ + ^£. (i.io) 

2nm* 2m* 

Finally, we note here that electron densities can be calculated from FS volumes 
determined with dHvA, and the resulting number should respect the stochiometry of 
the material under study. This is called Luttinger's theorem [49J. For Sr 3 Ru207, we 
will use this in two dimensions, where this translates to the sum of in-plane areas, 
with dHvA frequencies F n of n different sheets, 

"e = 2^£F„ = ^I>, (1.11) 

with S^ z and Fbz the in-plane area of the BZ and corresponding dHvA frequency, 
and the factor two stands for the two spin species. Note that for hole-like FS pockets, 
one should use the area of electrons, which is the area of the BZ minus that of the 
cyclotron orbit, or alternatively, F n = Fbz — Fhoie- The result of this sum should 
correspond to the number of electrons in the material minus the number of electrons 
in filled bands. 



10 For instance, oscillations have been measured in the magnetisation, resistivity, susceptibility, 
magneto-caloric effect and even sample size. 
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1.3.2 Oscillations at finite temperature 

We discuss in this section the role of temperature. At T = 0, the jumps in the 
DOS are produced by the crossing of the Landau levels through a perfectly sharp 
Fermi surface, which broadens proportionally to k B T as temperature increases. This 
results in smearing of the quantum oscillations, through an average weighted by a 
probability distribution, g(e) that corresponds to the negative derivative of the Fermi- 
Dirac distribution with respect to energy e, 

ts = d f( e ) = 1 (1 u) 

{e) de 2fc B T[l + cosh((e-//)/fcT)]" ' 

The oscillations reduce in amplitude, through a convolution of the signal xP0> X = 
1/B, by g(e) as a function of dHvA phase 0, which we normalise by an unknown 
factor A: 

X (X,T) = ^ X T^pFX + (f>)g(^d<f>. (1.13) 

The relation between e and <fi comes from slight differences in frequency, or dephasing 
of 



2npm* 

-Ae 



ehH 

with Ae = e — fi, fi being the chemical potential, and the effective cyclotron mass m* 
is taken from its definition, 

m = — . 1-14 

2vr de K 1 

We rewrite g(e) as a function of 0/A, with <f>/\ = Ae/^T: 

1 



9\ T 



\) 2^T[l + cosh(^-^0)] 



and obtain 

2Trpm*k B T 
e~hH ' 

which corresponds to the half width of the function g, and is specific to one Fermi 
surface sheet with quasiparticle mass m* . Eq. 11.131 is the Fourier transform of g: 



Re 



j A p exp \2-KipF n x\ exp(i0)#f^ 
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Figure 1.10: The factor LK(X,T) as a function of temperature, using typical values 
for the mass and magnetic field (m* = 7m e and H = 12T respectively). 



It is a calculation done using Cauchy's residue theorem, and is given in appendix [S] 
The result is a coefficient to the oscillations of the form 

~f\\ n ^ 
9W 



sinh(7rA) 

This function is called the Lifshitz-Kosevich relation (LK) [SHI EI], which we can 
rewrite as a function of the variables of interest, X and T: 

LK ( X,T)= c= 2j^ (1.15) 

sinn(Cm*Ai J e/i 

The LK prefactor is different for dHvA frequencies with different masses, but also for 
all harmonics p. For an orbit of frequency F with mass m*, its harmonics at pF will 
appear to possess a mass of pm* . 

Figure II. 101 shows the shape of LK(T). One can see that the oscillations die off 
at a temperature of approximately 0.148 m ^ m , the half-width of the function, which 
involves the mass m*. Moreover, it is possible to extract m* from dHvA using a non- 
linear fit of the amplitude of one frequency as a function of the temperature, which 
will be shown in section 12.2.41 The masses are usually different for each FS sheet 
since these originate from the hybridisation of different atomic orbitals with different 
dispersions, but can be extracted independently. 

One last word should be said about quasiparticle masses. These can be compared 
with the electronic specific heat of the material under study, in order to help evaluate 
whether a complete set of Fermi surface sheets has been identified in a dHvA exper- 
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iment. For two-dimensional materials, one can calculate that the linear specific heat 
coefficient 7 relates to the masses through a simple sum, 

7 = ^^-J2< = 1.48mJ/molK 2 5>;, (1.16) 

n n 

where a is the in-plane lattice parameter, and the factor 1.48 mJ/molK 2 was obtained 
using the lattice parameter a of Sr 3 Ru 2 O70. The values of m* are the masses of each 
FS sheet. The result of this equation possesses units of mJ/mol K 2 per formula 
unit. In our case, we will want to quote the result per mole of Ru, expressed as 
mJ/mol Ru K 2 , and a factor two will be taken out of the coefficient. This relation 
has been used with success in various systems, notably in Sr 2 Ru04 [35] and Sr 2 Rh04 



1.3.3 The role of disorder and impurities 

Crystal disorder and impurities also have an effect on the amplitude of the dHvA 
oscillations in a very similar way to thermal phase smearing. Disorder is temperature 
independent, and only affects the field dependence of the oscillations. Scattering of 
electrons on impurities or defects has the effect of a finite lifetime for Landau levels. 
If the scattering centres are distributed uniformly and fixed in space, the rate of 
scattering F is constant, and the probability of having an electron travelling without 
scattering during time t is 

P{t)dt = e~ rt dt, T = -, 

T 

In frequency (or energy) space, this corresponds to a lorentzian probability, 

r 

P{u)du = — — —dw. 

T 2 + (u - u ) 2 

The Landau levels are broadened, and the probability of having an electron in level 
n with energy between e and e + de is 

P{e)de = i + (e-L) 2 M 2 r 2 *- 

The effect is the same as for the Fermi surface thermal broadening, it corresponds to 
a convolution of the signal with a Lorentzian function, where the phase <p is defined 

11 For Si'3Ru207, the in-plane lattice parameters are equal, a = b. 
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as 

_ Ae 
A ~ ftf' 

where this time A is 

27rpcm*r 
A = ^~ • 

The prefactor to the oscillations is the Fourier transform of the Lorentzian, 

g(X) — exp(— |A|). 

This reduction function is called the Dingle factor, after its discoverer |53j. Shoenberg 
expresses the function as [15] 

D(X) = exp(-Cm*T D \X\), (1.17) 

where C is the same constant as in the LK function, and the Dingle temperature 
T D , assuming all the missing quantities, is expressed as a temperature but relates to 
disorder and impurities: 

k B 

A subtlety arises from the fact that the scattering rate does not depend exclusively on 
the mean free path £, but also on the size of the Fermi surface and the quasiparticle 
mass: 

i m*f 

where vp and kp correspond to the average Fermi velocity and kp vector along the 
cyclotron orbit. The Fermi wave vector in this equation is an average relating to one 
cyclotron orbit. The Dingle factor can be rewritten in the following ways: 

D(X) = exp ( _ = exp( _|^) = exp( _^), (1 ,8) 

where F is the dHvA frequency for one cyclotron orbit and Cf the circumference 
of that orbit in k space. The damping is consequently stronger for larger FS cross 
sections but does not depend on the quasiparticle mass. 

1.3.4 2D systems and the angular dependence of dHvA 

Quasi two dimensional metallic systems have very special characteristics in terms of 
dHvA. When metals possess a very anisotropic resistivity, the Fermi surface usually 
connects with itself through the Brillouin zone in the direction defined as the k z 
direction, as it is the case in Sr 2 Ru0 4 and Sr 3 Ru 2 7 . For Sr 2 Ru0 4 , the Fermi surface 
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has been determined with extremely high accuracy and a good review about how 
it was done was written by C. Bergemann et al. [35]. We reproduce here the points 
that are relevant to this thesis, and the main calculation, which is rather arduous, is 
presented in detail in appendix [B] The main characteristic that we observe is that 
when a material is close to being perfectly two dimensional, the Fermi surface takes 
the shape of one or many cylinders, oriented in the k z direction and connecting with 
the top and bottom of the BZ, with an arbitrary shaped base in the k x , k y plane (not 
necessarily a circle), and only small deviations along k z . This situation is generally 
accompanied by a real space lattice c parameter that is much larger than the other 
two, a and b, such that the overlap of the orbitals in that direction is small, resulting 
in a much lower conductivity in the c direction. In such a case, the Brillouin zone is 
shorter in the k z direction than in the other two. 

When one performs a dHvA experiment with the magnetic field aligned with the 
c-axis, the frequencies obtained correspond to the area of the base of each cylinder, 
and they are generally easy to identify to, for instance, k space areas calculated 
from angle resolved photoemission spectroscopy (ARPES) data. Moreover, if one 
tries to perform dHvA with the magnetic field aligned in the basal k x , k y plane, no 
oscillations are detected, since no closed Landau orbits exist in that orientation. dHvA 
frequencies are at a minimum in the c-axis direction and, as one rotates the magnetic 
field from the c-axis towards the plane, they gradually increase to infinity, and their 
amplitude decreases to zero. In quantitative terms, as one rotates the magnetic field 
the cross- sectional shape becomes stretched in one direction, by a scaling factor which 
corresponds to the inverse of a cosine: 



where F is the dHvA frequency in the c-axis direction. 

Since metallic materials are never truly two-dimensional, and the various Fermi 
cylinders have to connect to the top of the Brillouin zone at an angle of 90° for 
perdiodicity, it follows that they must possess at least two extremal areas in the 
k z direction. If the area difference between those is small compared to the spacing 
between the Landau levels at the magnetic field that is used, then the extremal area 
treatment of dHvA is not appropriate, and one has to use a slightly more elaborate 
approach. Moreover, the oscillations produced possess very special beat patterns, 
which evolve with field angle. We will now show that these patterns can be identified 
in order to work back and deduce the deviations of the true FS from simple cylinders. 

When the extremal orbit approximation is not applicable, one must calculate the 
following relation in order to obtain the fundamental oscillatory part in the magnet- 
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isation, for one Fermi surface pocket: 

Moc J C dk z sm( ^ nF ^ kz ^ oc J dK z sin (^A(k z )XJ , (1.20) 

where X is the inverse of the field, k z corresponds to k z x 27r//i& z , hf, z being the height 
of the Brillouin zone 4ir/c, c the lattice parameter, and A(k z ) is the cross-sectional 
area (not necessarily an extremal area) of the Fermi pocket at k z . For the complete 
dHvA signal, with the harmonics and the whole Fermi surface, one needs to sum one 
such term per harmonic, using p times the area A(k z ), and one per Fermi surface piece, 
using an appropriate area A n (K z ). This equation produces the beat pattern that is the 
result of the deviations from a cylinder in the k z direction, roughly the interference 
pattern from the difference between the extremal areas. When we rotate the magnetic 
field away from c-axis, these areas change, and so does the beating period. This 
produces a pattern in the field-angle plane that is unique to the corrugation of the 
Fermi surface pocket. It can thus be used to identify accurately those deviations. 
Consequently, what is required is to calculate this function of k z and field 

angle 9, and putting it into the integral, one should obtain the full dHvA signal as a 
function of X and 9. 

The deviation of A(n z ) from a cylinder is arbitrary, but it is a periodic function 
of k z , and one can expand it in a Fourier series. Moreover, since the in-plane shape 
of the Fermi surface is not known, it can be expanded in polar harmonics. Since this 
in-plane shape has a specific orientation in k space, one more parameter is required, 
which we express as an in-plane angle </>q. We then write this cylindrical expansion 
with subscripts, \x for the polar expansion and v for the n z Fourier series : 

oo 

k F (n z ,(f), O ) = k 00 + ^2 ^cosz/K z cos(/i0 + o ), (1-21) 

fl=l,U=l 



where the parts sinusoidal in vk z were dropped for reasons of symmetry 12 !. One is 
required to calculate the area difference between the real area and that averaged along 
k z , A A ft,,, as a function of k z and 9. This is a difficult calculation that may be of 
interest to some readers, and it is presented in appendix [B] The result is 

27ik 00 kf tu Jf l (vK F taii9) , 

A/l^ = — — cos/i0o x ( — 1) 2 cos vk, z , even /i, 

cos 9 

2irko kf lu Jf l (vK F ta>n9) , ^it^il . 

■cos/i0o x ( — 1) 2 smvK x , odd //, 



cos# 



1.22) 



12 



In the Brillouin zone, electrons at k and — k possess the same energy. 
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where J M is the /xth Bessel function of the first kind. For a Fermi surface pocket that 
possesses such a complex structure that more than one warping parameter k^ u are 
significant, then one must sum all the AA^ V and put them into the integral of eq. 11.201 
in order to calculate the magnetisation. This cannot be performed analytically but 
is done with a computer. Figure fLTTl right, shows an example of a pattern produced 
by a simple warping parameter hot which has 0.25% of the size of the in-plane Fermi 
wave vector kp. 

Crystal symmetry restricts the possible values of the warping parameters. Contri- 
butions with fi = are allowed in all 2D crystal symmetries, since they are isotropic 
in the plane. But parameters with fi > must respect the in-plane symmetries since 
they provide the Fermi surface its in-plane shape, where for example a square involves 
parameters with /i = 4. Consequently, in Sr 3 Ru20y, for electron Fermi pockets in 
the centre of the zone, odd parameters or parameters that do not respect twofold 
rotation symmetry are forbidden. Effectively, /i = 1 produces a surface that has 
no symmetry at all, and since the Brillouin zone must possess inversion symmetry, 
E(k) = E(— k), it follows that such a surface cannot be in the centre of the BZ. I 13 l 
Warping parameters of /i = 3, 6, 9... may only appear in hexagonal systems, hence not 
in any of the ruthenates, and higher odd numbers cannot be accommodated in any 
periodic system. Finally, the terms with v = 0, /i > are invisible to dHvA, where 
effectively, the in-plane shape cannot be determined. For Sr 3 Ru 2 07, we will discuss 
only parameters with v > 0, \i — 0, 2, 4.... 

Furthermore, while the warping situations where fi = are isotropic, for all warp- 
ing within the class with u > 0, /i > 0, some directions of rotation will produce no 
beat patterns at all, due to the factor cos/i0 o - For instance, with k 2 i, depending 
on the orientation of the Fermi surface pocket, one of the situations, 0o — 45° or 
0o = 0°, should produce an area AA^ V that vanishes. This property may be used for 
the identification of such warping 14 !. Any difference in patterns between directions of 
rotation should be attributed to such parameters. 

If only one warping parameter is relevant, as it is sometimes the case, then one 
can go further and calculate the integral of eq. 11.201 with single warping k^, which 



13 If such a surface is not at the centre, it should appear at least twice, with inversion symmetry 
through the centre. Then, with dHvA and rotation, one should measure both at the same time, but 
with a difference in cf>Q of 180°, and due to the factor of cos n<po, it should cancel exactly, and cannot 
appear in the data. 

14 For example, in the case of Si'2Ru04, the a pocket possesses a term in &2i, which has the 
symmetry of a screw, or as C. Bergemann calls it, the shape of a "snake that swallowed a chain" 
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Figure 1.11: Left Angular dependence of a dHvA frequency produced by a simple 
corrugation of type fcoi- F° r both figures, red corresponds to high values, green to 
intermediate and dark blue to low values. Right Field-angle beat pattern produced 
by a simple corrugation of type fc i- 

is given in appendix [Bj and the result is 
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where Jo is the average dHvA frequency, and AF^ the difference. The modulation of 
the dHvA oscillation is not a cosine but a Bessel function. When one uses small fields, 
though (when the argument of Jo is larger than about 4), the modulation becomes 
very close to a cosine function of the form 



{2wAF^X 



7CX 



COS 



I" 



cos 9 



-J M (z//c otan( 



7T 



;i.24) 



There exists an angle at which the argument of the modulation factor becomes zero for 
all values of X, which arises when the angle 9 is such that the function J A1 (i//c otan^) 
vanishes. In that case, there is no modulation and the oscillatory signal is maximal. 
This is called the Yamaji angle, and for one specific type of warping, it is only function 
of fcoo, the average in-plane hp. Figure [T.ll[ right, shows an example of a beat pattern 
produced by a simple corrugation of type k 01 of 0.25% of k F = 2.33 x 10 9 m _1 , 
corresponding to a frequency of 1.8 kT, in a Brillouin zone of height 606xl0 7 m _1 , 
identical to that of Sr 3 Ru207. One can see the first Yamaji angle at 45°, and the 
second one at 65°. 
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Figure 1.12: Left Real dHvA frequency, proportional to the respective sizes of the 
Zeeman split FS. Right Observed dHvA frequency, after the back projection. 



In dHvA frequency space, peaks evolve with angle with the form 
F{9) 



_Fq ± tan( 

cos 9 COS 9 



where, since AF^ V is small, one observes two peaks that oscillate around Fo/ cos 6. At 
the Yamaji angles, the second term becomes zero and both peaks cross, the amplitude 
becoming large. Figure |1.11[ left, shows a simulated example of this phenomenon, 
where dHvA data was generated in the same way as for the beat pattern shown on 
the right, but using a larger warping of 2%, for clarity of presentation. 

Finally, we show one more calculation that will be used later, which aims to find 
curves of B(9) where the amplitude of the modulation is either maximum or zero. 
Using the cosine approximation for the Bessel function of eq. I1.24[ the condition for 
the modulation to vanish is 



2nAF, u X 7T 

r — J^i/zcootanfl) - 7 

cos 9 4 



±.mr 



n 



1,2,3... 



which translates as 



B{6) 



8AF^ U J M (z/ k 00 tan 9) 



(1.25) 



cos 9 1 ± 4n 

where the ± corresponds to positive for maxima, and negative for minima. Using 
n = 1, one obtains the first zero near the Yamaji angle, and with larger n, all the 
other zeros going to lower and lower field values as they intercept the 9 = axis. 
Note that in the case of a single warping parameter, beat pattern nodes cannot cross 
one another. 



Chapter [I]: Scientific Background 



33 



1.3.5 Spin splitting 

The previous sections have neglected the spin of the quasiparticle and Zeeman split- 
ting of the FS leads to an additional amplitude reduction factor that requires at- 
tention. As one ramps the magnetic field, the FS splits into two, one for each s pin 



flavour, one growing and the other shrinking in k space. In normal spin splitting | 15 | . 
when a paramagnetic metal possesses a constant magnetic susceptibility, the value of 
the average Fermi wave vector changes linearly in field, 

kf = k F ± xiB, 

where Xu the rate at which hp changes with field, is proportional to the spin suscept- 
ibility, k space cross sectional areas increase with B, and so do the dHvA frequencies. 
In 2D systems, such areas are also be rescaled by the cosine of the angle 9 between 
the sample and the field. Neglecting the second order term, for one cyclotron orbit, 
we obtain 

F n (B, 9) = ± — ^B = ± — ^— — B. 1.26 

cos 9 e cos 9 cos 9 2 cos 9 m e 

When put into the oscillatory term of dHvA, we obtain an additional phase <fi = 

cos 8 m e ' 

This has no effect on the dHvA frequency spectrum, since only the phase changes. 
However, it affects the amplitude if the term cos <fi vanishes, which in 2D systems is 
called a spin zero, and the condition for that happening is 

ng in* 7i m* 
' cos 9 = 2g — , 



cos 9 m P 2 m. 



e 



One can thus find in the data a set of angles where field independent zeros appear, and 
obtain the g factor in this waj^. This additional amplitude modulation corresponds 
to interference of dHvA oscillations between the spin species. 

It is not the same situation when, in a certain field region, B c , a metal undergoes 
a superlinear rise in magnetisation, as in a metamagnetic transition. In this case, F^ 



15 For details of anomalous spin splitting, see the review by Bergemann [35] pp. 658. 

16 This phase of the spin zero, proportional to gm*, refers to non- interacting systems. In Fermi 
liquid systems, the spin susceptibility is enhanced by interactions, and one should use an effective 
g* factor, which is equal to . g pa , i 7 !? being one of the Landau Fermi liquid parameters. The phase 

then becomes <p — -^-g "^pi" , where m* remains the unrenormalised (band) effective mass. 
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will increase faster than linearly 1 7 I. 

F n (0, B) = J- ± -^(xiB + X ,B% 
cos 6 e cos 

and one is left with terms periodic in B (or higher orders of B), 

cos — — cos + 2tc -B . 

\ B J \ e cos J 

In such a case, the observed momentary frequency, measured from various Fourier 
transforms taken as a function of 1/B, corresponds to the real frequency minus its 
linear part jSH [55] , 

dF 

F obs = F — B-. (1.28) 
dB 

For instance, if in a system, one of the spin FS sheets evolved with field faster than 
linearly, one would observe the behaviour described in figure 11.12^ where the real 
frequency is depicted on the left and what is observed on the right. Below the crit- 
ical field B c , where simple paramagnetism is present, no difference in frequency is 
seen between both spins, but near B c , one might detect splitting of the peak in the 
spectrum, which becomes very large at B c (shown in red), and decreases again at 
high fields (shown in blue), but never disappears. This means that in metamagnetic 
systems like CeRu 2 Si 2 [SH], UPt 3 [SS], Sr 3 Ru 2 7 [10] and others, splitting is expected 
on the high field of the superlinear increase in magnetisation. 



1.3.6 First dHvA measurements in SraRu^Oj 

dHvA was measured in Sr 3 Ru207 previously by Borzi and co-workers [40J and the 
main properties of the FS have been explored. Quantum oscillations were obtained 
with the magnetic field parallel to the c-axis and five quasiparticle orbits were reported 
in the low field side of the metamagnetic transition, and four in the high field side. 
Figure [T. 131 left, shows both low and high field side spectra calculated from data in 
the ranges 5.5 to 6.5 T for the low field side, and 10 to 15 T for the high field side. 



Table [T~2l presents the value of the observed frequencies. Peak splitting was observed 
in the high field side, as one expects in a system with non-linear susceptibility (see 
section 11.3.51) . but the picture that emerged was one that is similar in both field 
sides, with peaks located at similar frequency values, which possessed similar mass 
values^£|. Moreover, closer to the metamagnetic transition on both field sides, some 
of the frequency peaks were reported to increase sharply, shown in the inset to figure 



17 Where \i is much larger than the prefactor to the second order term that was dropped in eq. 
18 The mass values for the peaks of the high field side were not specifically quoted in the text. 
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frequency (kT) (T) 

Figure 1.13: Left dHvA spectra from the measurements by Borzi and co-workers 
[3D] . The red trace refers to the field range 5.5 to 6.5 T, and the blue one to 10 to 
15 T. Right Mass divergence of the 1.8 and 4.2 kT peaks, as reported by Borzi. 



5.5 to 6.5 T 


10 to 15 T 


Peak (kT) 


m* /m e 


Peak (kT) 


m* /m e 


0.2 


6±3 


< 0.4 




0.5 


13 ±4 






0.8 


10 ±1 


1.1 - 1.5 




1.8 


8± 1 


1.6 




4.2 


9±2 


4.0 





Table 1.2: Table of the frequencies observed by Borzi, below and above the meta- 
magnetic transition, with associated quasiparticle masses. 



11.131 left. The frequencies had no field dependence in the field ranges that were 
quoted as used for the FFTs of the plotted spectra. 

Moreover, away from the metamagnetic transition, quasiparticle masses were cal- 
culated from the temperature dependence of the quantum oscillations, in the same 
field ranges used for the plotted spectra. Those were fairly high, in the range of 7 to 
12m e , and are presented in table [TT2"1 Borzi also observed that they have a field de- 
pendence, and increase as one moves towards the QCEP from both field sides. Figure 
11.131 right, shows the field dependence of the mass of two of the frequencies, where 
the vertical line represents the metamagnetic transition, and masses close to 30m e 
were observed. This observation has attracted a fair amount of attention. It is an 
important measurement for the quantum critical picture in SrsPa^CV 



Chapter 2 

Experimental methods and 
analysis procedures 

This chapter aims at describing in detail the experimental techniques and numerical 
methods that were used during this research project. These can be divided into two 
categories; those that relate to the de Haas van Alphen experiment and those that 
are part of the procedure for sample characterisation and sorting of single crystals of 
Sr 3 Ru207. 

Firstly, the dHvA experiment is a standard one, but some of the fine details are 
often overlooked. Thus, the author found it pertinent to present the relevant parts 
with a lot of detail. In particular, some aspects of the data analysis procedures are 
somewhat obscure in the standard textbooks, but also, special steps can be taken in 
the numerical methods to tackle our specific aims. We will emphasise particularly 
practices that may lead to systematic errors. 

Secondly, the most important part of this chapter may be the description of the 
main experiment, but it could not have been carried out successfully on average 
Sr 3 Ru207 samples. The purity requirements of dHvA are extremely high, and they 
had not yet been achieved by material scientists in the early stages of this project. 
Indeed, the highest purity was reached by R. S. Perry in a collaboration with the 
author over several months, the information that was provided by the characterisation 
and sorting procedure of section [27^1 being used to optimise crystal growth conditions. 
Moreover, the best qualities were only reached in some spatial parts of single crystals, 
which required to be isolated. In consequence, an extensive string of measurements 
were applied to each sample cut individually in order to search for the best ones. 

This chapter first presents the dHvA experiment using the field modulation method. 
Next will be given the complete description of the various data analysis procedures 
and the rigorous checks that were made to ensure proper thermal equilibrium of the 
system. The second category of techniques will be introduced at this point, first de- 
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scribing the AC susceptibility probe that was built for the adiabatic demagnetisation 
refrigerator. The complete procedure for characterisation and sorting of samples will 
follow with the other already established techniques, along with data relating to the 
purest samples. 



2.1 AC susceptibility in a 3 He/ 4 He dilution refri- 
gerator 



We introduce in this section the main method that has been used in this experimental 
project, AC susceptibility. This technique can be used for several purposes, and in 
respect to this fact we begin with a general presentation of its operation through 
equations. We were, however, interested in measuring quantum oscillations, and carry 
on with a description of the case where the magnetisation is an oscillatory function 
of the magnetic field. Finally, we provide details of the specific experimental system 
that was used to obtain the data at the core of this thesis, along with quantities 
related to signal and noise levels. 

2.1.1 Basics of AC susceptibility with equations 

Magnetic measurements using AC susceptibility is an established technique. The 
principle is to induce oscillations in the magnetisation of a sample and to measure 
the change in the magnetic flux using a coil. In order to improve sensitivity, one 
usually connects another empty identical coil to the first in opposition, in order to 
make the measurement differential. When measuring susceptibility as a function of 
magnetic field, one uses a small oscillating field Hac on top of a large slowly varying 



field H 



DC- 



H(t) 



H DC (t) + H AC (t), 
at + HAcsin{uit). 



The magnetic induction in the first coil is 



B(t) 



H(t) + M(t), 



and the voltage on its terminals is 




dM 



dt 
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The voltage on the second, empty coil is equal to the first term only, and the sub- 
straction will leave a differential voltage proportional to the differential magnetic 
susceptibility 

dM dH A ■ . . 

AV = A »°dW IT = A ^ ^ 

The magnetisation as a function of magnetic field is usually a non- linear function and 
the excitation field Hac will produce harmonics, which can be detected with a phase 
sensitive detection device. Expanding the magnetisation as a function of Hac-, one 
obtains 




One requires its derivative as a function of time. Assuming that the time derivative 
of the DC component is small compared to the oscillating terms, one finds 



dM , JdM\ H\ c , , ,Jd 2 M\ 



H ac uj cos(wt) + ^ £ ^( 1 - sin(2wt)) 



dt v '\dH)„ 2 v y JJ \dH 2 



+ -|^w(cos(3a;t) + cos(a;t))f + ... (2.2) 

One may observe two facts in this equation: first, the nth term includes the nth 
harmonic and those are multiplied by the oscillating field to the power n. Second, 
the phase changes by 90° at each harmonic. 



2.1.2 Oscillatory signal in the field modulation method 

This section presents the calculation performed by Shoenberg ([15], pp 103-105), 
with emphasis on aspects that will be used later. During dHvA experiments, the 
magnetisation follows oscillations as a function of inverse magnetic field, of frequency 
F, which are assumed to be sinusoidal 0: 

M(H) = M (H) + A sin ( ) , 

\at + Hac sm{ut) J 

where M represents the magnetic background and is not oscillatory as a function 
of magnetic field. This function is quite awkward to use, and a form that involves 
harmonics for phase sensitive detection is required. One uses the fact that the os- 
cillating field amplitude is smaller than the magnitude of the DC field (usually by 

^^When they are not, each oscillation harmonic, not to be confused with detection harmonic, can 
be treated independently. 
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factor between 10 3 to 10 5 ), 

1 1 H AC (t) 

H DC (t) + H AC (t) ~ H DC (t) H 2 DC (t) ' 

so that 



M(H) ~ M (H) + Asm ( J cos (Asin(wt)) 

\HdcJ 

—A cos ( — — J sin (Asin(cjt)), 
V Hdc J 



with 

27rF# AC 



A = (2.3) 



11 DC 



One then faces the complicated factors cos(A sin(a;t)) and sin(Asin(u;t)). These are 
periodic functions which can be expanded in Fourier series. This calculation is given 
in appendix O and the result involves Bessel functions of the first kind Jjt(A), plotted 
in figure I2.lt 



— - J ( J (A) + 2 V J 2fe (A) cos(2A;wt) 

Hoc J v t^i 

-A cos J ^2 ^ J 2fc -i(A) sin ((2A; - . 

Each harmonic has as a prefactor, a Bessel function of order k and the magnetic 
oscillations are in phase with Hac for the odd harmonics and ? out of phase for the 
even ones. 

In order to calculate the differential voltage, one requires the time derivative of 
M(t). In real systems, M usually features much slower field variations than the 
oscillations, and all higher derivatives are negligible, except, for instance, very close 
to a phase transition. Neglecting those, one obtains the harmonic expansion: 

dM dM 0/ TT . , N „ /2ttF\^ 

= [ a + riAcU cos(ut) j + 2aMsm I — — > 2/cJ 2 fc(A) sin(2reu;t) 

a# \HdcJ ~ 



dt 



( '2irF\ °° 

-2uA cos — - V(2A; - 1) J 2 fc_i(A) cos ((2Jfe - l)wt). 



(2.4) 



One may conclude from this result that the slowly varying magnetic background 
is only present in the first harmonic, and it can be useful to use second or higher 
harmonic detection in order to remove it from the oscillatory signal. 

It is an interesting fact that the dHvA amplitude is proportional to a function that 
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X Magnetic field ( T ) 



Figure 2.1: (Left) Bessel functions of the first kind. The case of n — 1 is approx- 
imately linear at low values of A, and the n = 2 case is quadratic for small values of 
A. (Right) Bessel functions of the first kind Jfc(A) as a function of the DC magnetic 
field He>c f° r k — 1,2, for a frequency F = 1.8 kT and an oscillating field value of 
H AC = 1QG. 

oscillates with the AC field value and it is worth examining why. It corresponds to a 
smearing effect of the dHvA oscillations when the peak to peak value of the AC field 
is close to the DC field interval it takes for a specific dHvA frequency to complete a 
cycle. The condition for zero amplitude on this frequency is: 



2H AC = AH DC} 

and the period corresponds to 

1 1 1 AH 



DC 



F H 2 Hi Hp C 



The condition for zero amplitude is then 

2H AC F A 

This is approximately true, since the bessel function Ji(A) has its first zero at A ~ 
3.81. The reason is that the DC field continuously changes, and it takes a value for 
the oscillating field to be just a little over half the DC field interval. 

Finally, also shown in appendix [Cj when A < 1, the time derivative of M reduces 
to a simpler expression, since 

- 2^T (2-5) 

indicating that each harmonic k has an amplitude proportional to H\ c / H D k Cl a fact 
that will be used later. 
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Figure I2TTI shows Bessel functions of the first kind J*. (A). From the left panel, one 
can observe that for low values of A, the first harmonic has a much higher amplitude 
than the first, but above A ~ 2.2, the second harmonic possesses a higher intensity. 
Thus, when the magnetic background is important, one may gain significantly in 
using a high modulation field and second harmonic, assuming that no eddy current 
heating is present (an effect that will be discussed in section |2~!6|) . The right panel of 
fig- 12. II shows the behaviour of the k = 1, 2 bessel functions for fixed dHvA frequency 
and oscillating field, as a function of DC field. One can see that as A is function of 
the DC field, the system can pass through a zero of the bessel function. It is then 
important to set the oscillating field such that for all DC field and dHvA frequency 
values, the system does not cross a zero of the nth bessel function when using the 
nth harmonic. 

2.1.3 The dHvA apparatus 

The dHvA apparatus that was mainly used in this project is located in the Cavendish 
Laboratory of the University of Cambridge and operated by the Quantum Matter 
researchgroup. The system is composed of a dilution refrigerator of base temperature 
of 6 mKf, which features an 18 T superconducting magnet @. The cryostat is located 
on a vibration isolation base, made of several tons of concrete placed on top of a 
rubber layer, contributing to very low vibrational noise levels. 

The samples were placed inside coils situated on a rotation mechanism, inside the 
magnet bore. Figure l2T2| left, shows a photo this system. It is composed of a frame, 
small bobbins onto which are fixed the compensated coil pairs, which can rotate 
around, and a rod that, by being pulled or pushed, controls the orientation of the 
bobbins, and spans around 90°. The rotation was operated outside the cryostat by an 
electrical motor, the procedure and calibration of which is explained in appendix [El 
Figure I2.2[ left, presents a sketch of the sample mounting, heat sinking and coil 
system. The pick-up coils were wound in Cambridge by Swee K. Goh. The pairs 
were well compensated, with coil resistances of 361 Q and 361 O for the first pair 
used with sample C698I, and 350 Q and 348 Q for the second coil used with sample 
C698A, and all had approximately 1000 turns. The samples were heat sunk using 
three gold wires each, bonded to the material using sintered Dupont 6838 silver paste, 
10 cm long, connected to a 1mm thick high purity silver wire tightly screwed against 
the mixing chamber of the dilution refrigerator. As we will see in section 13.1. 3[ the 



2 Due to an additional factor of 4 in ea 12.41 for k — 2 compared to k = 1. 

3 The base temperature is only reached in the best conditions, which was not the case in this 
project. 

4 Details of the dilution refrigeration cycle can be found in the book by Pobell [57] . 
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Figure 2.2: Left Photo of the plastic probe used in the Cambridge dilution refriger- 
ator. Top Right Schematic representation of a heat sunk sample of Sr 3 Ru 2 7 . Right 
Schematic representation of a compensated coil with a sample as was used in the 
susceptibility system. The rotation axis is shown with a curved arrow. 
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thermalisation of the samples might not have been perfect at the lowest temperatures, 
but good above around 90 mK. 

The modulation field was produced using a high quality power audio amplifier, 
to which was fed the oscillatory output voltage signal from a SR830 lock-in amplifier 
(LIA). A low modulation frequency of 7.9 Hz was used in order to obtain penetration 
depths longer than the sample size, but also, as we will show in section I2.6[ when 
using second harmonic, for constant sample heating by eddy currents, gains were 
made in signal amplitude by using as low a modulation frequency as was possible. A 
careful study of eddy current heating, presented in sections 12.61 and T3.ll revealed that 
the optimum modulation field to use in this specific case was of around 20 G. 

Each coil in a pair was connected in parallel with a variable resistor, using closely 
twisted wire triplets. Adjusting the resistors effectively changed the magnitude of 
the signal from each coil, and allowed improvements to the compensation of the coil 
pairs. The resulting signal was fed to low temperature transformers which were main- 
tained at 1.5 K, and amplifed the signal by a factor 100. In order to damp noise at 
higher frequencies than the measurement one, 0.1 uF capacitors were added in par- 
allel with the low temperature transformers. The resulting signal was fed to floating 
preamplifiers of amplification factor 1000, which contributed to matching the imped- 
ance between the low temperature transformers and the LIA, and an improvement 
in signal amplitude of approximately six was obtained in this way, normalised by the 
amplification factor. 

Throughout the measurements, an experimental procedure designed to reduce 
noise levels associated with the operation of the magnet power supply was used. By 
adding a large power resistor of small resistance value, the superconducting magnet 
system was transformed into an RL circuit (see figure l2~3| Top left). When a current 
is present in such a system, along with a voltage at the power supply output, as 
one changes abruptly the voltage, the current adjusts exponentially to the new value, 
the field energy being absorbed by the power resistor. One may express this with a 
differential equation, where the current is /, the resistance R and the inductance of 
the magnet L, with solution I(t): 

V = RI + L d i 
dt 

I(t) = (Ii-If)e-fr + I f , 

with Jj and If = Vf/R the initial and final values for the current, Vf being the new 
voltage value. Figure |2T3"| Top right, sketches the profile of current when the voltage 
is abruptly changed from one value to another. With such an analog current sweep 
method, the magnetic field produced was more stable than with the normal mode and 
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Figure 2.3: Top left. RL circuit that represents the magnet system. Top Right. Time 
dependent profile for the magnet current when one abruptly changes the voltage on 
the leads of the magnet. Bottom. Effect of the voltage limited mode on the noise 
levels in a dHvA measurement. The system mode is changed at 10 T from the voltage 
limited to the linear one, when sweeping downwards. The change of noise levels, when 
normalised by the amplification factors, at 10 T is from 100 to 500 pV/y/Hz. 

the noise levels picked up by the measurement coils became at least five times lower. 
It was not usable with the current system at low magnetic fields, since the sweep 
rates became too slow. Consequently, in the work presented in chapter [3], the long 
field sweeps from 18 to 2 T were performed in three parts, two at high fields, from 18 
to 15 T and 15.1 T to 10 T, using the voltage limited mode at two different voltage 
values, resulting in sweep rates of around 0.04 T/min, and from 10.1 T to 2 T using 
the normal mode with a constant sweep rate of 0.02 T/min. Due to time constraints 
and the fact that the oscillations have longer periods in field at high fields than at 
low fields, a higher sweep rate was used in the high field sweeps. A careful procedure 
was created in order to join the resulting data files using the regions of data overlap. 
Figure I2T3"1 Bottom, shows a strip of data between 9 and 11 T, just near the region 
where the change from voltage limited to linear modes occurs, at 10 T, and a huge 
change in noise amplitude can be seen. 

As a brief conclusion, we quote the noise levels that were obtained. When not using 
the voltage limited sweep mode, the noise levels normalised by the total amplification 
factor of 100 000 were typically of 500 pV/V Hz (peak to peak) and very stable. 
At high fields, when the voltage limited mode was used, the noise levels were of 
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100 pV/ V Hz (peak to peak). In the high field data, the noise was five times lower 
than in the quietest measurements performed previously with the current St Andrews 
system, but combined with the improvement to the signal due to a few more electronic 
components (the variable resistors, the capacitors and the preamplifiers), the total 
improvement in signal to noise, normalised by the square of the modulation field, was 
even higher. 

2.2 De Haas van Alphen data analysis 

This section aims to present the details of the numerical analysis methods of dHvA 
data that were used in this project. This experiment usually produces huge amounts 
of data, which contains a lot of information, but complex numerical methods are 
required to be used in order to extract the parameters one is normally looking for. 
Moreover, since this experiment is at the core of this project, the author felt it relevant 
to describe also problems and pitfalls that exist in the analysis of quantum oscillations, 
which are not always well known by many condensed matter physicists. Such problems 
can lead, and have led in the past, to the production of systematic errors. A good 
understanding of this subject has been imperative throughout this project, and is 
hopefully well reviewed in this section. 

We first review the analytical expressions to the various terms involved, which 
have been derived in the last chapter, in order to set a convention for the language 
used in this thesis, and sketch how these vary with temperature and field. We then 
describe how discrete Fourier transforms are performed using windowing. Since Four- 
ier transforms presented in other work are not always appropriately normalised, but 
often feature arbitrary units, we felt that a description of how this was done in this 
project was not superfluous. One of the important pieces of information provided by 
dHvA is the quasiparticle mass, and we therefore carry on with a description of the 
non-linear fits involved in its extraction. The use of large field windows may lead in 
some case to systematic errors in the extraction of the quasiparticle mass, and we fur- 
thermore present how this can be avoided. Finally, the oscillatory amplitude contains 
valuable information which can be obtained by a method that we introduce. In par- 
ticular, such data can provide knowledge of the mean free path of the quasiparticles 
through the Dingle factor. 

2.2.1 Analytical expression for the oscillations 

We review in this section the analytical expressions for the different terms in the 
equation for the quantum oscillations, which were derived in the first chapter, and 
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Figure 2.4: Three dimensional representation of the Lifshitz-Kosevitch factor as 
a function of temperature T and inverse magnetic field X, for a constant mass 
m* = 7 m e . The shaded region illustrates a typical large X measurement inter- 
val, where the magnetic field extends between 6 and 15 T. The white line indicates 
Xq, the average X value for this interval. 



1 





Figure 2.5: Illustration of the Dingle factor as a function of X. It is symmetric around 
infinite magnetic field (X = 0) and zero at zero field (X — > oo). 
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their dependence upon the different physical parameters that come into play in these 
experiments. This is done in a compact and consistent form that makes later descrip- 
tions of the numerical methods easier. 

We write the expression for the oscillations in the form of four functions, J*., LK, 
D and E, relating to the Bessel function modulation (eq. 12.4)) . the Lifshitz-Kosevitch 
(LK) factor (eq. 11.15)) . the Dingle factor (eq. II . 17)) and the finite size of a real data 
set. We express with the variables X and T the inverse field and the temperature, 
respectively, and with the integers n and k, the various FS sheets and the detection 
harmonics. Moreover, we use F n as the reciprocal variable to X, the oscillatory 
frequency of FS sheet n, X as the average inverse field and m* as the quasiparticle 
mass of the nth FS sheet. The complete expression for the oscillatory part of the 
susceptibility in a dHvA experiment is 

oo 

Xoc Jk(X,F n )LK n (X,T,m* n )D n (X)E(X - X )cos(2irF n X). (2.6) 

n,k=l 

As seen in section 12.1 .21 (and appendix [C]), the oscillations are modulated by Bessel 
functions of the first kind 

which when A n < 1 reduces approximately to 

This factor is present in dHvA but absent in Shubnikov-de Haas oscillations. 



The second part of expression 12.61 is the LK(X, T) function, which features the 
only temperature dependence of that equation, 



sinh Cm* XT' 



with variable C containing various physical constants. When using m* as a quasi- 
particle mass enhancement without units, 

C = ^^14.7T/K. 

en 

Figure E3] illustrates the dependence on T and X of the LK factor. It is identical 
for both variables, but values that can be reached in an experiment are very different. 
Often, a small region in X will be used, over which the LK factor varies very little, 
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and X may be safely replaced by its average X , shown as a white curve. In a case 
where a large X interval is chosen, depicted with a shaded area, the dependence on X 
has to be considered. The implications of this will be shown in the next two sections. 

The third factor modulating the oscillation amplitude is the Dingle factor D(X), 
plotted in figure 12.51 

D(X) = e~ XnW , 

with the constant A n proportional to the scattering rate of quasiparticles in FS sheet 
n. 

The last coefficient in eq. I2.6l is the function E(X). Its origin lies with the interval 
of X over which data is available, where the data is assumed to vanish outside of that 
region. In order to improve the clarity of the Fourier spectra, this function is usually 
chosen such that is vanishes smoothly at each ends, a practice described in the next 
section. Without the function E(X), the oscillations would extend to infinity and be 
symmetric around zero. 

2.2.2 Fourier transforms and windowing 

In order to study the Fermi surface structure of a material, one is required to sep- 
arate the different contributions at different frequencies from its various FS sheets, 
using Fourier transforms. Unfortunately, it is not obvious what the analytical Four- 
ier transform of equation 12.61 is. We may nevertheless show here why windowing is 
needed, which type we chose and its effect on the Fourier spectra. 

The convolution theorem states that the transform of a product of functions is 
equal to their convolution in reciprocal space, 

J r [f(x)g(x)] = F[f(x)\ * F[g{x)\ = f(x) * ~g(x). 

Applied to eq. 12.61 expression, one obtains 

F[ X (X, T)} cc E{F) * (d u (F) * (LK(F, T) * (5(F - F n ) + 5(F + F n )))\ 

where the Fourier transform is taken as a function of X and yields functions of the 
frequency F . 

In a dHvA experiment involving short X intervals, one assumes D and LK con- 
stant, and the Fourier transform is proportional to 

n ± 



Chapter Experimental methods and analysis procedures 



49 



the Fourier transform of the envelope E centred at ±F n . If the data of a dHvA 
set extended from zero to infinite magnetic field, E(F ± F n ) would be a Dirac delta 
function centred at ±F n , but with finite size data, E(F ± F n ) possesses a width and 
a shape defined by the nature of E(X). If one takes E(X) as a top hat function of 
width 2a centred at Xq, equal to one where dHvA data exists, and zero everywhere 
else, then E(F ± F n ) is 



E(F) = / e~ iFX dX = -^=ae- iFXo sinc(aF). 

V27T JXn-a V27T 



"X a +a 
'X Q -a 

The sine function oscillates around a central peak, and taking modulus of this result, 
the satellite peaks could be misinterpreted as additional dHvA frequencies. One 
generally uses another form of windowing for E(X) that decreases smoothly at the 
edges of the data set. An article by Harris reviews most types of windowing techniques 
[58] , but in general, it is the Hann or the Hamming windows that are the most often 
used. The first is defined as 



E(X) = .5 + .5 cos 

= otherwise 



-(X-X ) 
a 



—a < X — Xn < a 



Its Fourier transform is 



E(F) = -Le iFX ° sinc(aF)- """ 



/27T 7r 2 - a 2 F 2 



and decreases faster than the sine function away from the main peak, having much 
smaller satellite peaks. The Hamming window, on the other hand, is tuned with a 
parameter a to suppress to zero the first satellite peak on each side of the main peak 
and is obtained by adding a very small top hat window with a Hann window: 



E(X) = a + (1 - a) cos 
= otherwise. 



-(X-X ) 
a 



-a < X - X < a 



The value of a ~ 0.54 gives the best results. Its Fourier transform is 

E(F) = -|=e^ sinc(aF) (a + (1 - «)^^) , 

The Hamming window was the one applied in this work. 

Along with windowing, one generally also uses the technique called zero padding. 
It consists in artificially increasing the length of the windowing function E(X) with 
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values of zero in order to reduce the coarseness of the Fourier spectrum. It is equivalent 
to an appropriate interpolation of data points between existing ones, such that the 
spectrum is smoother and more continuous. 

Finally, we may add a comment about the width of the Fourier transform of the 
windowing functions. The half width at half maximum of the sine function is found 
by equating 

2a sin(aAF) 
v/2^r aAF 

sin(aAF) 

which yields numerically aAF = 0.3017, where AF is measured from the centre of the 
function. The same can be done with the Hann and the Hamming windows, which 
give, respectively, 

AF = 0.3017/a Top Hat 

AF = 0.5/ a Hann 
AF = 0.4538/a Hamming 

These widths are function only of that of the window E(X), regardless of, for in- 
stance, the number of cycles of oscillation within that set. For one to be able to 
compare spectral information at differing frequencies, the Fourier transforms have 
to be performed with equal inverse field widths, not the same number of cycles of 
oscillation. 



2 

aAF 



2.2.3 Appropriate normalisation of Fourier transforms 

Throughout the data analysis of this thesis, it has been important to be able to 
compare the amplitude of various dHvA signals measured on different samples of 
Sr3Ru207, which has required a rigorous understanding of the normalisation proced- 
ure of dHvA spectra, with appropriate units. We present in this section the exact 
analytical method with which one should normalise and present dHvA spectra, and 
how to extract information from them. 

As one measures voltages and takes Fourier transforms as a function of inverse 
field X in T" 1 , the frequency units are T and one expects that the units of the Fourier 
transforms should be V/y/T. The Fourier transform AV of the voltage measured 
AV is generally a complex number data set. The power spectrum W(F) represents 
the power radiated per unit frequency, and is the square of the modulus divided by 

5 In reference to time signal processing, where the spectral amplitude has units of V/VHz. 
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the interval (that we denote 2a): 

W{F) = Ya^ V{m2 = 2^ AVAV ' 

which has units of V 2 /T. In general, it is not the power spectrum that is plotted but 
its square root, with V/v^T. In order to know the voltage amplitude of the signal for 
a certain frequency interval, the power spectrum is integrated within that interval, 
which yields units of V 2 , and a square root of the result is taken. Due to spectral 
leakage, in order to know the voltage amplitude of a certain dHvA frequency, one is 
required in theory to integrate the power spectrum from minus infinity to infinity, 
since E(F) is non-zero for all F. In practice, the presence of multiple frequencies 
makes this impossible and one may integrate a peak using its full width at half height 
as boundaries. This yields a fraction of the real amplitude, that depends on the 
integration limits but also on the windowing function used in the Fourier transform. 
The voltage measured on a coil may be expressed as 



AV(X) = J2^nCos(F n X), 



with amplitude coefficients A n of units of V. The Fourier transform one performs is 
AV(F) = -= / E(X)AV(X)e iFX dX, 



oc 



/2tt J-c 
Y^A n E(F±F n ), 



and 



Y^A n E{F±F n ) 

n,± 

One usually requires the amplitudes A n , assuming that the frequencies are well enough 
separated such that the various dHvA peaks E(F — F n ) do not overlap. Since the 
coefficients A n are real, one requires the value of the integral of \E(F)\ 2 around zero 
in order to find A n . In the case of a top hat windowing function, the integral is 

/ \E(F)\ 2 dF = / sinc 2 (aF)dF = 2a, 

and the total power of the nth frequency is W n = A 2 n . Since one generally integrates 
only the peak at positive frequency F n , an additional factor of 2 is multiplied to the 
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integrated power spectrum, and the formula to use when extracting A n from AV is 



A,= v^=j2£i-|AV(F)|MF. 



In the case of a Hann window, the integral is 



and the formula is 



A "=f\L^ Kv(F)ViF - 

Finally, with a Hamming window, the integral is 



OO 1 

\E(F)\ 2 dF = 2a(a 2 + -(l-a) 2 ) 

-oo ^ 



and with the usual parameter a = 0.54, the result is t^cl = 0.7948a. The formula 
to use is 



^= W2^ ri|AV(F)|W. 



1987 2a' 

These formulae enables one to obtain average dHvA amplitude values over spe- 
cific field ranges, and do not depend on the length of the data set if the amplitude 
remains constant. Thus, amplitudes from sets of differing lengths may be compared. 
In this thesis, quoted amplitudes were always obtained using this method, with a 
Hamming windowing scheme, and was proved correct by feeding simulated data to 
the calculation and obtaining the right result. The voltages thus calculated were 
directly comparable with voltage noise levels, measured in V/vTIz. 



2.2.4 Mass analysis 

One of the parameters of primary interest in dHvA studies is the quasiparticle cyclo- 
tron mass. This value is extracted from the temperature dependence of the amplitude 
of the oscillations, and one such mass can be obtained from each frequency in the 
dHvA spectrum, and can consequently be assigned to a specific branch of the FS. 
In this section, we explain the numerical method used to extract the quasiparticle 
masses. 

One assumes once again that the functions LK(X,T) and D(X) vary slowly 
enough and can be taken out of the Fourier transform, by replacing X by X . This 
is true most of the time and we will show in the next section the situations where 
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it is not possible to make this approximation, in other words, when the inverse field 
window is too large compared to D and LK. One calculates the Fourier transform of 



eq. UM 

X(F, T) ex LK n (X , T)D n (X )e iFXo E(F ± F n ) 

n,± 

The temperature dependence is only present in the function LK n (X ,T), 

where m* is the mass enhancement normalised by m e . In order to extract m*, one 
is required to perform a non-linear fit of the LK function. This fit is carried out on 
values of the square root of integrals calculated over peaks in the power spectrum, as 
described in the last section, as a function of temperature. The fit possesses only two 
adjustable parameters, A and m: 

If in some cases a dHvA spectrum possesses a constant additive offset, identical for all 
temperatures, it is possible to use a three parameter non-linear fit, with parameters 
A, m and C, 

We will show in section 14.4.11 that, in contrast to the two parameter LK fit, the use 
of this scheme for the extraction of the mass can lead to large systematic errors, and 
was possibly responsible for the mass enhancement measured by Borzi and co-workers 
[4Ti] which has not been reproduced by the author. 



2.2.5 Fourier transforms using large field windows 

There are cases where the inverse field window one needs to use in order to calculate 
spectra with a desired resolution is too large compared to the width of the functions 
LK(X) or D(X). In this section, we demonstrate that in such a situation, the normal 
procedure of fitting LK(X ,T), the LK function at the average inverse field, always 
underestimates the real quasiparticle mass. This problem is generally overlooked, but 
the solution presented here has been critical for the analysis of dHvA data from one 
material, Sr 2 Rh04 [59], where the total quasiparticle mass calculated with the stand- 
ard method did not agree with the electronic specific heat. The analysis presented 
in this section was created by the author and contributed to finding of the "missing 
mass" of Sr 2 Rh0 4 . 
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When one works with well separated frequency peaks in the dHvA spectrum, one 
is not very limited in possible sizes of the inverse field window with which to perform 
a Fourier transform. Normally, noise is the limiting factor, as enlarging the inverse 
field width averages more data and reduces the noise floor in comparison to the signal. 
The situation where one is required to use large field windows is if dHvA peaks are 
too close to one another to be integrated separately when using suitably small field 
windows. 

Parseval's theorem states that for an amplitude function and its Fourier transform, 
the area under the power spectrum is equal to the area under the square of the 
function. In the case of dHvA, the integrated power spectrum is the average of the 
square of the product of the amplitude functions, 

W = 2a(\E(X - X )LK(X,T)D(X)\ 2 ) 
= 2a(E 2 (X - X )LK 2 (X,T)D 2 (X)). 

Supposing for a moment that D(X) is constant and that E(X — X ) is a top-hat win- 
dow, one observes that the T dependence of \/W is proportional to 2a-\J (LK 2 (X } T)^, 
which is not equal to LK(Xq, T), since the width in temperature of LK(X, T) varies 
with inverse field X (depicted in figure 1^41) . Consequently, this averaging produces an 
Li^-like function that cannot be adjusted adequately by a non-linear fit of LK(X , T). 
Moreover, this average is always wider in T than LK(Xq, T), the LK function at the 
average field X . Consequently, the fit of LK(X , T) always generates a mass that is 
smaller than the real value. 

Assuming that D and E are not constant worsens matters, since these functions 
influence the average of the LK function by weighing differently the amplitude at 
different values of X, E towards the centre of the window and D towards higher 
values. Nevertheless, as values at higher X produce a higher contribution, the mass 
extracted by a fit of LK(X ,T) still underestimates the value. 

Figure |2"1)1 shows LK fits taken on simulated dHvA data, using the normal method 
of eq. 12. 7\ as a function of the position of the lower field boundary of the inverse field 
window, the higher one fixed. We observe that as the lower boundary moves towards 
the higher one, the mass approaches the real value. 

Two methods of analysis may be used for such circumstances. The first is is to 
use the appropriate fit function, which is 

^2a(E 2 {X - X )LK 2 {X, T)D 2 (X)). 

This requires knowledge of the Dingle temperature To for the function D(X), the 
method to extract will be shown in the next section. The second is to find the mass 
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Figure 2.6: Relative value of the quasiparticle mass extracted with the normal fit of 
the the function LK(Xq,T) applied to simulated dHvA data, to the real value. The 
relative mass is plotted as a function of the position in B of the lower boundary of the 
field window, the higher boundary fixed at 15 T, where the oscillations are the highest. 
The frequency used was of 1 kT and the mass of 3 electron masses, parameters similar 
to those of Sr2Rh04. For a vanishingly small window, the calculation error is nil, and 
increases with the size of the window. The error can be higher than 50%. 



in a self-consistent way using a simulation. In this procedure, one uses data simulated 
with the real Dingle temperature and a chosen quasiparticle mass, performs a Fourier 
transform using the same large field window as with the experimental data. One then 
carries out an LK fit and adjusts the mass of the simulation until the LK fit yields 
the value obtained with the experimental data. The adjusted mass of the simulation 
corresponds to the real quasiparticle mass. This is the method that was used by Perry 
and co-workers [59], as a result of the insight given in the present section. 

2.2.6 The envelope extraction method 

We describe here the method which will be used in section 13.3.31 in order to extract 
the field dependence of the amplitude of the oscillations, taken from the review on 
Sr 2 Ru0 4 by Bergemann et a/.j3S], the envelope extraction. The envelope corresponds 
to the product of the functions Jk, D, LK and some interference patterns, as discussed 
in section 11.3.41 Using this technique, which consists of filtering, one can obtain 
continuous functions of the modulation amplitude for all individual dHvA frequencies. 
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The dHvA oscillations take the form 

X (X) = Y t MX)cos(F n X), 

n 

The functions A n (X) are the envelopes one wishes to extract. Carrying out a Fourier 
transform yields 

X(F) = W MX)(e i(F - F)X ~ e«W)dX, 

= ^7^Y} A ^ F - ~ + F n )] 

v n 

One requires to isolate one of the various A n (F), denoted A m (F), filtering the spec- 
trum by using a frequency window that leaves only one peak and sets the rest to zero. 
If one also sets the negative frequencies, which appear symmetrically to the positive 
ones, to zero, one obtains, 

Xflt(F) = -±=A m (F - F n ). 
An inverse Fourier transform is then used to calculate the desired A m (X): 
XMX) = ±- j_J m (F')e^ X dF' = ^A m {Xy 



)e lF - x . 



By taking the modulus, one is left with the envelope, without oscillations, 

\Xfit(X)\ = ±-\A m (X)\. 

It is possible to use this method even if the frequency peaks are split or vary with 
inverse field, where one includes everything in the filter window, where it is effectively 
the case with Sr 3 Ru207 close to the metamagnetic transition. The shape of the filter 
window is moreover important for the same considerations as for the function E, 
where essentially the envelope extracted by this method consists of the real envelope 
convolved with the Fourier transform of the filter function. In particular, using a 
top-hat filter leads to oscillations added to the extracted envelopes. One can damp 
this effect by using a top-hat filter with "rounded corners''^: 



V2a 



(( F-a ,F + a s 

erf( — — — ) — err(- 



V2a y/2t 



a 



'In exact terms, a top-hat function convolved with a gaussian of width a. 
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where a is the roundness of the filter and erf is the error function. The rounder the 
filter, the fewer oscillations are present in the result. 

The usefulness of this method is multiple. One can in principle fit the LK relation 
at each value of the inverse field X. This makes the calculated quasiparticle mass 
absolute for a value of X, instead of an average over the window E used in a Fourier 
transform. If the mass is not a constant of the field, one can, in principle, calculate 
its field dependence using this method. The inconvenience is that it produces noisy 
temperature distributions and fits with large error bars. More importantly, the main 
interest in this method is in the analysis of modulating factors to the oscillations 
other than the LK function. As we will show in the next section, it is possible to 
extract the Dingle temperature with this method. But mainly, it will be used, section 
I3.3.3[ in order to obtain interference patterns as a function of magnetic field and field 
angle. 



2.2.7 The Dingle analysis 

When the temperature is strictly zero, the LK function is a constant of the inverse 
field. In an experiment, depending on the quasiparticle mass, it is often possible 
to find a temperature for which there will be an accessible region of inverse field 
where LK(X,T) is very near or equal to one (High fields, low temperatures). When 
one works with a system without beat patterns or Bessel function coefficients (as 
in Shubnikov-de Haas measurements), using the envelope extraction method of the 
previous section, one can extract the Dingle temperature using a simple exponential 
fit of the form 

D(X) = e- xm , X = U.7m*T D} 

and extract the Dingle temperature To- 

In general, the task is slightly more complex, as it is not possible to remove the 
Basse, or beat pattern coefficients! Thus one is compelled to fit a function which 
features all the contributions to the modulation of the oscillations. For dHvA data 
with a beat pattern due to simple warping (section [1.3.44 eq. 11.24)) . this takes the 
form 

AD(X)LK(X, T)J k {X) cos(uX + 0), 
CXTm* 

= Ae~ xlxl ————— | cos( W X + 4>) J k {2 * pi * H AC FX 2 )\, 

where the variables A and a are free parameters. The constants u, <f), Hac, F, m*, 
C and T are known values. The mean free path is extracted from A through the 



7 Thcse cross zero in certain field regions, and taking them out leads to divisions by noise. 
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equation 



2.3 AC susceptibility in an adiabatic demagnetisa- 
tion refrigerator 

The de Haas van Alphen experiment, as we know from section [I.3.3[ depends on the 
impurity scattering rate of the electrons exponentially. As we explain later in section 
12.44 obtaining good samples in this respect was not trivial, but improved dHvA data 
immensely compared to using samples of average quality. Consequently, a sample 
sorting procedure was developed by the author. We will show that for high purity 
Sr 3 Ru207 samples, mean free paths determined using the residual resistivity is not 
representative of bulk samples required for dHvA. The dHvA experiment itself may 
be used for that purpose, but the use of the apparatus described in section 12.1.31 is 
time consuming^. This motivated the development in this project of an alternative 
AC susceptibility system for an adiabatic demagnetisation refrigerator (ADR), which 
could accommodate large unprepared samples. It was not intended for precision 
measurements but for ease of use and time of operation. 

This system has served two different purposes, discussed in section 12.41 The 
first was to determine the amplitude of the dHvA signal as a measure of the mean 
free path, and the second was to determine the superconducting volume fraction, 
which corresponds to that of Sr2Ru04 inter-growths, by detecting zero DC field AC 
susceptibility as a function of temperature. This section describes the characteristics 
and testing of this new experimental probe, along with a short introduction to the 
ADR system. 



2.3.1 The adiabatic demagnetisation refrigerator 

The adiabatic demagnetisation refrigeratoi^ (ADR) relies on the very high difference 
in entropy of a paramagnetic system between its states in high and low magnetic 
fields. The entropy is lowered by applying a high field, and when the system is put 
in thermal isolation, reducing its magnitude lowers the temperature by a transfer of 
entropy from the sample and measurement apparatus to the paramagnetic material. 

Figure 12.71 illustrates the main parts of the ADR. It is put into a dewar filled 
with liquid He 4 . It features two independent stages that can be thermally connected 
or isolated, and possesses two cooling mechanisms, one on each part. The first is 



8 Performed in a dilution refrigerator, it has a turnaround time of 2-3 days. 
9 Cambridge Magnetic Refrigeration model mFridge mF- ADR/50. 
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Figure 2.7: Sketch of the ADR components. The left side of the figure shows the top 
part, featuring the 1-K stage, and the right side the LTS. 

a standard 1-K pot, which refrigerates the apparatus by pumping 4 He gas from its 
liquid surface, lowering its vapour pressure and temperature down to 1.5 K (left part 
of the figure). The second corresponds to a pill of paramagnetic salt situated in the 
bore of a 6 T magnet. 

Connecting both stages, by cooling down the system to 1.5 K using the 1-K 
pot while applying the maximum field of 6 T to the paramagnetic salt pill, the low 
temperature stage (right part of figure 12. 7p is put in its state of low spin entropy and 
high temperature. Then, opening the thermal connection between both stages, the 
low temperature stage (LTS) may be stabilised at any temperature between 0.1 and 
1.5 K by adjusting approximately linearly the value of the demagnetisation field. A 
heat leak is present, though, and the field is required to continuously decrease in order 
to compensate the heat input. Cooling is not possible anymore when the field reaches 
zero, and all temperature values will possess a different hold time, of zero at 0.1 K up 
to several hours near 0.5 K. Appendix iDl presents a calculation of the cooling power 
of an ideal ADR, where one finds that, effectively, the constant entropy is a function 
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of H/T. 

The ADR can moreover perform experiments with a sample magnetic field, inde- 
pendent of the demagnetisation field, produced by an 8 T superconducting magnet 
(lower right in the figure) at the sample region, at constant temperature. The LTS 
also features a compensated field region, where the thermometry is located, between 
the two magnets. The temperature is measured using two different thermometers, 
one for high and the other for low temperatures. For temperatures between 300 
and 1 K, a Cernox semiconductor resistive thermometer is used, and for temperat- 
ures between 1 K and 50 mK, the system switches to a Ge semiconductor resistive 
thermometer. Next to these is also situated a heater that can be used to obtain tem- 
peratures higher than the 1-K pot temperature, such that the system may be used 
continuously between 15 and 0.1 K. 

2.3.2 AC susceptibility apparatus for large crystals 

Figure [2781 shows the schematics of the AC susceptibility probe for the ADF, designed 
and built by the author in order to possess as many pick-up coils, along with a 
modulation coil, as the space inside the sample magnet on the LTS allowed. 

All parts were made of tufnoj^j or copper, except for supporting stainless steel 
threaded rods. The pick-up coils were made with a large bore diameter of 4 mm in 
order to accommodate large unprepared crystals, and the table in figure 12.81 lists all 
their parameters, where resistance values were obtained at room temperature. They 
were wound using copper wire, of 125//m for the drive coil and 50/im for the six pick-up 
coils. The coils were numbered 1 to 6 from the top to the bottom of the probe, and the 
pairs 1,2; 3,4; 5,6 were connected in opposition. The sensitivity refers to the voltage 
on a coil in a field of 1 G at 1 Hz. When the probe is installed on the LTS, the pairs 
are normally connected to low temperature transformers located on the 1 K stage of 
the ADF, with an amplification factor of 100. This increases the signal to noise level 
of the measurement apparatus, and the noise levels range between 100 pV/v / Hz and 
1 nV/vTlz (peak to peak and normalised by the amplification factor). 

Up to three samples may be placed inside the pick-up coils by fitting them inside 
a copper braid, which is inserted from the bottom of the apparatus, and subsequently 
thermally anchored to the low temperature stage. High vacuum grease is used to 
improve thermal contact with the braid. The thermal equilibrium of the samples 
with the system was verified, shown in the next section. 

1 A hard plastic used for cryogenic applications. 
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Figure 2.8: top Schematics of the AC susceptibility probe for the adiabatic demag- 
netisation refrigerator. Bottom Technical data for the various coils of this apparatus. 
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2.3.3 Thermal equilibrium tests and calibration 

In order for the AC susceptibility measurement probe to be used reliably, we have 
characterised its performance. Two properties were investigated, the thermal equi- 
librium during experiments and the calibration of voltage to magnetic susceptibility. 
Both measurements were performed using the diamagnetic signal of superconductors. 
We present in this section the method and data which were used for these purposes. 

In order to verify the thermal equilibrium of the samples during an experiment 
that uses temperature as a variable, we chose the measurement of the superconducting 
phase transition of different metals as an absolute thermometer. We used three metals 
of different critical temperatures, aluminium, zinc and Sr 2 Ru04, which possess critical 
temperatures (T c ) of 1.196, 0.875 [12] and 1.430 Kr 1 ! respectively and critical fields H c 
of 105 G and 53 G for Al and Zn, and 750 G (B || c-axis) and 1.5 T (B || ab-plane) 
for Sr 2 Ru04 (601 E3]. AC susceptibility was measured on samples of these materials, 
which had shapes of rolled foil aligned with the oscillating field, for Al and Zn and 
1 mm 3 cubes for Sr 2 Ru04, at zero sample field as a function of temperature. 

Figure 12.91 presents the in-phase (red curves) and out-of-phase (blue curves) data 
measured on these samples. With Al and Zn, susceptibility was measured with two 
different temperature sweep rates, and we obtained curves with hysteresis, indicating 
that a small residual magnetic field was present. We obtained critical temperatures of 
1.175 K for Al and 0.880 K for Zn, within 20 mK of literature values^, with hysteresis 
of 10-15 mK. In the case of Sr 2 Ru04, the critical temperature was of 1.53 K ( B | 
c-axis) and 1.55 K (B || ab-plane), with h yst eresis of 25 mK, 120 mK higher than the 
value measured by N. Kikugawa in Kyotd_£. 

From these values, we demonstrated that no significant temperature gradient was 
present, in the region from 0.8 to above 1.5 K, between the samples and the thermo- 
meter. Effectively, any gradient would result in the samples being warmer than the 
thermometer, since the thermometer is situated between the cooling mechanism and 
the samples. The temperature at the thermometer position, from the superconduct- 
ing critical temperatures, was either very near or higher than expected, suggesting a 
possible temperature calibration or thermometer self heating problem. 

For further tests at lower temperatures, we measured the temperature dependence 



of dHvA in Sr3Ru 2 07 (figure [279| bottom two plots) and in Sr 2 Ru04 (not shown), 



n Private communication from N. Kikugawa, who measured T c of the same Sr2Ru04 samples in 
Kyoto. 

12 In a magnetic field, the superconducting transition is of first order. Hysteresis is present, and T c 
measured using both upwards and downwards field sweeps are lower than the zero field value (see, 
for instance, [57], PP 205-207). The value measured with an upwards field ramp is always closest to 
the zero field value. 

13 Later measurements on samples of Sr2Ru04 grown in St Andrews reproduced these results. 
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Figure 2.9: Top and middle left In-phase and out of phase susceptibility as a function 
of temperature for rolled foil shaped samples of high purity aluminium (top) and zinc 
(bottom) (99.999%). The data was taken using two different sweeping rates, 0.11 
K/min and 0.04 K/min, upwards and downwards in temperature. Top and middle 
Right Susceptibility as a function of temperature for samples of Sr 2 Ru04, one oriented 
with H || c-axis (top) and the other with H || afr-plane (bottom). Bottom left dHvA 
oscillations and spectrum in Sr 3 Ru207, sample C698A. Bottom right Amplitude of 
dHvA as a function of temperature, stars, and LK fit, solid line. 
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between 0.25 and 0.45 K (Sr 3 Ru 2 7 ) or 0.8 K (Sr 2 Ru0 4 ). Using LK fits, we obtained 
quasiparticle mass values for the 1.8 kT peak of Sr3Ru207 and the 3.0 kT peak of 
Sr 2 Ru0 4 of 8.4±0.1 and 3.1 ±0.1 respectively. Again here, had the temperature of the 
samples been higher than that of the thermometers, the half- width of the LK function 
may have been larger, and the mass, proportional to the inverse of the width, would 
have been smaller than the real value. Since we know from the literature [4"0~t |35] that 
the masses are of 8 ± 1 and 3.3 ± 0.3, for Sr 2 Ru0 4 and Sr 3 Ru 2 7 respectively, we 
conclude that there was no significant temperature gradient below 0.8 K. 

Finally, from the diamagnetic signal of Al, we calibrated the voltages measured 
onto the coils into approximate real susceptibility values. We mainly used this inform- 
ation for calculating volume fractions of Sr 2 Ru0 4 into Sr 3 Ru 2 7 samples (see section 
12.4. 3p . in the small fraction limit. From Figure [279]. we obtained the amplitude of 







the diamagnetic signal of two samples with a demagnetisation factor close to zerq_ 
The voltage measured on a pair of compensated coils, with a sample placed into one 
of them, with an oscillating magnetic field B is AV = ARBx, where A is the total 
area within the turns of the coil containing the sample and R the fraction of this area 
intersecting the sample, which is also equal to the volume fraction of the sample with 
respect to the coil, R = v s /v c , and x is the susceptibility. 

The real coil volume is not an obvious quantity to measure, and for that we used 
the fact that when empty, the total coil area A is equal to the measured sensitivity 
S given in the table of figure 12.81 



15 



SfB AC v s 



Using a voltage difference of 1.9 mV for Al from figure 12.91 which includes an 
amplification factor of 100, an average sensitivity of 11.1 /iVG _1 Hz _1 (or, in area, 
11.1 xl0~ 2 m 2 ), a frequency of 71 Hz, a modulation field of 0.428 G, a mass of 
52.6 mg and a density of 2.702 g/cm 3 , we obtained a value for the volume of the 
coil of v c = 0.346 cm 3 . This was used to simplify our equation and we deduced an 
approximate calibration, 



4/ cm 3 HzG N AV 
x = 3.12 x 10 4 — — -— — 

V B ac Hk 



V ] B AC i 

We conclude by noting that this was not used as an accurate calibration, but only as 
a measure the order of magnitude of volume fractions. 



14 Effectively, their shape was of cylinders of rolled Al or Zn foil, which allows magnetic flux to 
penetrate between the layers. 



15 Since AV" = $ = AB AC , we have A = AV/B AC f. 
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2.4 Systematic procedure for the search for high 
purity samples 
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Figure 2.10: Standard recipe for cutting crystal rods grown in the image furnace 
for the various experiments in the sorting procedure. The rods were cut into slices 
1.8 mm thick. These were then further cut to obtain square prisms of dimensions 
1.8 x 1.8 x 3 mm. These parts were intended for the final dHvA experiment and were 
subjected to the susceptibility and magnetisation experiments. The off-cuts were 
used for measuring resistivity and the X-ray spectrum. The residual resistivity and 
the X-ray spectrum of the square prism were extrapolated from those of the off-cuts. 
The best two crystal prisms were finally polished into octagonal prisms, heat sunk 
and used in detailed dHvA experiments. 



We mentioned in the introduction for this chapter that the material purity require- 
ments for dHvA experiments are very high. For this work, I was fortunate enough to 
interact with the materials scientist who was growing the crystals of SrsRu207, R.S. 
Perry, since this was done in the laboratory next door. The strategy we adopted was 
one of feedback style, where I carried out a chain of characterisation experiments on 
grown samples and provided the information back to R. S. Perry, in order to improve 
the growth parameters and restart the process. The sample sorting procedure was 
carried out using two types of characterisation methods. First, we will describe those 
used to detect different phases of the ruthenate family. Second, we will present the 
measurements of scattering from impurities and crystalline disorder. 

In the first dHvA experiments of this project, we used an average quality sample 
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of Sr 3 Ru 2 7 . We obtained relatively good quantum oscillation data, but used a very 
high modulation field in order to obtain a good signal to noise ratio, which led to eddy 
current sample heating. This sample also led us to erroneous interpretations when 
we discovered a dHvA frequency at 12.7 kT, which belongs to Sr 2 Ru04. It became 
clear that the sample was not purely of the bilayer type but had small inclusions of 
Sr 2 Ru04. This led us to suspect that it might also contain Sr 4 Ru30io, SrRuC>3, etc. 
Moreover, these inclusions of Sr 2 Ru04 possessed a purity comparable to that of the 
Sr 3 Ru 2 07, since we detected the oscillations well. Consequently, any good sample 
for dHvA should be thoroughly checked for all possible impurity phases before its 
oscillation spectrum can be interpreted. 

We used the magnetic properties of the different ruthenate phases to identify them 
quantitatively. The reason being that only Sr3Ru 2 C>7 is paramagnetic at all temperat- 
ures. The n = oo, SrRu0 3 and n = 3, Sr 4 Ru 3 Oi , members of the Ruddlesden-Popper 
series are ferromagnetic with Curie temperatures below room temperature, of 165 and 
105 K respectively (281 ESI [26], EI]. We also expect the phases with n > 3 to be fer- 
romagnetic, but no clear evidence of their presence was observed in this project. We 
used a Superconducting Quantum Interference Device (SQUID) magnetometer at 
small fields for detecting the spontaneous ferromagnetic moment of our samples. And 
finally, the n = 1 phase, Sr 2 RuOzi, is a superconductor with critical temperature of 
1.48K or less ^ We used AC susceptibility at zero field to detect the diamagnetic 
volume fraction of the samples. 

The next subsections will describe all the characterisations that were applied to 
the samples, in the same order that they were performed. The order was important 
for time saving. The first and quicker methods determined whether the samples 
were good enough to be submitted to the next. The first experiment performed 
on a newly grown single crystal rod was powder X-ray diffraction, which identified 
all the crystalline components of the material, but was also destructive, performed 
with a small piece ground to a fine powder, and the content of the crystal rod was 
extrapolated from that result, to first order. The second experiment was that of 
resistivity as a function of temperature, where one used a long and thin plate of 
crystal mounted in the way that will be described. When extrapolated to absolute 
zero, one obtains a value proportional to impurity and crystalline disorder scattering. 
Again, disorder in the larger crystal rod was extrapolated from the result. Crystal 
rods that had good properties were divided into as many samples as possible with 
the right dimensions for the dHvA apparatus of section l2.1.3l Next, the samples that 
appeared promising were put into the ADR, for the measurement of AC susceptibility 



16 Quantum Design MPMS model. 

17 For complete information on Sr2Ru04, see [23] . 
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and dHvA, where we measured the Sr 2 Ru04 volume fraction and the scattering, 
respectively. Finally, the ferromagnetic content of the samples was determined with 
the SQUID magnetometer. 

This procedure was performed on more than 50 samples from around 10 crystal 
rods. 23 of them made it through the five experiments and were kept for this project. 
Figure 12.101 illustrates the way the crystal rods were divided for the various experi- 
ments. Data from the sample that we used in the early stages of this thesis project 
for dHvA experiments was kept as a standard with which to compare and against 
which to improve. The last subsection will present a summary of the data obtained 
with these samples. We will also demonstrate that the X-ray spectra and the re- 
sidual resistivities agree only approximately with the other measurements, although 
they were still good starting points. Finally, the best two samples were retained for 
the dHvA experiments described in chapter [31 where unprecedented signal quality 
was achieved. Moreover, this sorting procedure turned out extremely useful to other 
experimentalists. In particular, five of the crystals of tables 12.11 and 12.21 were given 
to Angle Resolved Photoemission Spectroscopy (ARPES) experimentalists A. Tamai, 
F. Baumberger and co-workers, of St Andrews, who produced results that have been 
essential to this project through experiments performed at Stanford, presented in 
section I4.1[ and two were used for detailed specific heat and magnetocaloric measure- 
ments by A. Rost, in St Andrews. Others were distributed for various projects, and 
their characterisation data has proven a critical asset. 



2.4.1 Powder X-ray diffraction 

Powder X-ray diffraction was measured first on all newly grown crystal batches. It 
took half an hour to perform and detected most phases of strontium ruthenate when 
present in large volume fractions. This method is of course destructive, so it was 
carried out on a small amounts of single crystal ground to powder. This measurement 
was always performed first in order to ascertain that the main phase produced during 
crystal growth was predominantly S^R^Oy. 

The measurement was produced with Philipps copper X-ray source and diffracto- 
meter, which emits at a wavelength of 1.54 Aand detects the reflection intensity as 
a function of angle. The various atomic planes of the crystal reflect the X-rays at 
specific angles, which are tabulated in standard databases 18 !. Figure 12.111 shows two 
examples of spectra, offset for clarity, the first featuring only Bragg peaks associ- 
ated with Sr 3 Ru207 (in blue), and the other one exhibiting additional peaks, from 
Sr 2 Ru04 and Sr 4 Ru30io (in black). The chosen angle range was between 5° and 35°. 



18 



From the International Centre for Diffraction Data. 
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Figure 2.11: Two examples of powder X-ray diffraction measurements, offset for 
clarity. The top spectrum is of a crystal containing some detectable Sr 2 Ru04 and 
Sr 4 Ru 3 Oi intergrowths, while the bottom one was measured on a crystal where these 
phases were not detected. 

The reason behind this choice was that higher angles are populated with too many 
peaks from Bragg planes of Sr3Ru207 and their identification is difficult. The main 
peaks that belong to Sr 3 Ru 2 07 between 5° and 35° were located at 8.6°, 17.2°, 23.3°, 
31.6° and 32.6°. When present with Sr3Ru 2 C>7, Sr2Ru04 had only one visible peak 
which was easily distinguished from the others, at 13.9°, and Sr4Ru 3 Oio had two, 
at 12.4° and 18.6°, the second of higher intensity. Unfortunately, SrRuC>3 featured 
peaks that were too close to those of Sr 3 Ru 2 7 or not intense enough to be used for 
its identification and none of the X-ray peaks of this phase were clearly observed. 

Since each X-ray spectrum featured an arbitrary intensity factor, only amplitude 
ratios between various peak intensities were meaningful. Consequently, when meas- 
uring the amount of Sr 2 Ru04 or Sr 4 Ru30io with this method, we took the ratio R, 
in the first case, of the peak at 8.6° of Sr3Ru 2 C>7 with that of the one at 13.9° of 
Sr 2 RuC>4, and in the second case, the peak at 8.6° again with the peak at 18.6° of 
Sr4Ru30io- Note that among the 23 samples that were kept for later experiments, 
only one featured a content of Sr 4 Ru30io detectable with powder X-ray diffraction, 
C6102A, shown in figure I2TTT1 It thus meant that all the measured samples contained 
relatively small amounts of this type of intergrowths. 

2.4.2 Residual resistivity 

The standard method to determine disorder in a layered metallic sample is to measure 
its residual resistivity, in the aft-plane direction. This configuration was used since the 
ao-plane is the direction in which Sr3Ru 2 C>7 cleaves readily into thin plates, making it 
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Figure 2.12: Left Illustration of a mounted resistivity sample in a four point resistivity 
measurement configuration. The sample is oriented with the afr-plane in the horizontal 
direction. Right Typical measurement of the resistivity as a function of temperature 
(blue), plot as a function of T 2 , with a linear regression of expression 12.111 (broken 
line). Its extrapolation to zero temperature yields the residual resistivity po- 



possible to obtain samples of very small cross-section 19 !. As was discussed in section 
ll.l.l[ in a highly correlated metal, resistivity as a function of temperature at low 
temperatures may be expressed as 



p(T) = Po + AT 2 . 



(2.11) 



The T 2 term is fundamental to the Fermi liquid theory, the factor A depending on 
parameters of the metallic state like the quasiparticle mass m* 2 and the Fermi wave 
vector fcp. In this case, we were interested in the constant po, which depends on static 
lattice disorder, 

m*r Hkp 
^° ne £ne ' 

with n the density of carriers and T the scattering rate of the electrons due to lattice 
imperfections and impurities, t the associated mean free path and kp the average 
Fermi vector. It is standard to measure disorder in units of pfl cm, and one simply 
needs to adjust the low temperature data with a linear regression of T 2 , extrapolating 
the data to absolute zero. Previous crystal growth experience led us to believe that 
only Sr 3 Ru207 samples that featured a po of lpQ cm or less were likely to exhibit 
good quantum oscillations. The best known Sr 3 Ru207 samples had po — 0.25pX2 cm. 

An illustration of a sample mounted in a four point resistivity measurement in 
the afe-plane configuration is given in the left side of figure 12.121 Very thin single 
crystal were cleaved from bulk samples, and then mounted on quartz substrates such 
that they were held in mid-air by thin gold wires. 50 pm thick gold wire was used 
for current and voltage contacts, fixed with high temperature cured silver epoxy, 
shown in light grejj 2 ^. This system allows for the relaxation of thermal stress during 



19 Plates of thickness 20-100 /jm were obtained 
20 Type Dupont 6838. 
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cooldown through the flexibility of the gold wires. In order to obtain a good signal to 
noise ratio, the sample geometry needed a very high aspect ratio. Typical dimensions 
were of 2 mm x 0.2 mm x 0.02 mm. Six contacts were used and two resistivity 
measurement were made per sample, one on each side of the crystal. The quality of 
the electrical contacts was determined by comparing the measurement performed on 
both sides, which normally yielded the same po- The most common problem that 
arose was a c-axis contribution to the resistivity, where bad quality electrical contacts 
at the current leads force the current to travel in the c-axis direction, for which the 
resistivity tensor component is around 300 times higher at 0.3 K j2TJ, and leads to 
erroneous residual resistivities. The way by which this was avoided was to make sure 
all contacts covered the thickness of the sample. 

The measurements were performed in a simple He flow cryostat, which can cool 
down from room temperature to 3.5 K. For Sr 3 Ru 2 0y, this base temperature was 
low enough to obtain a reasonable amount of data in the T 2 dependence region of 
p(T) that enabled us to extrapolate it to absolute zero. The measured resistance was 
converted to resistivity using the room temperature value of 232 uVl cm. The right 
side of figure 12.121 shows a typical resistivity measurement and T 2 extrapolation. 

2.4.3 Zero field AC susceptibility 

The most sensitive method to determine the presence of the superconducting phase 
Sr 2 Ru04 in samples that we intend to use for dHvA experiments is to measure AC 
susceptibility at zero magnetic field, as a function of temperature below 2 K. The 
transition from the metallic state to superconductivity is observed by a jump in the 
in-phase AC susceptibility (x') and a peak in the out-of-phase AC susceptibility (x") 
(see figure [2791 for AC susceptibility data of pure Sr 2 Ru0 4 ). In the limit where the 
demagnetisation factor of these inclusions is near one, the amplitude of the jump in 
x' corresponds to the volume fraction of superconducting material. In our case here, 
since we were searching for crystals without inclusions of Sr 2 Ru04, we were only 
looking for the order of magnitude of the volume fraction. 

The measurements were done with the opportunity that it was made in the same 
apparatus as for the dHvA experiments of the next section, and consequently, with 
the same cooldown of the ADR. It also implied that the samples were positioned with 
their c-axis parallel to the axis of the coils, as is required for dHvA. Due to time 
constraints, we were not able to measure each sample with B parallel to the afr-plane. 
We expected, from polarised light optical microscopy experiments done by Fittipaldi 
et al. [31] . Sr 2 Ru04 inclusions to have a geometry of flat islands perpendicular to the 
c-axis, and the c-axis of Sr 2 Ru0 4 to be in the same direction as for the Sr 3 Ru 2 07. Due 
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Figure 2.13: Typical zero field AC susceptibility measurement. The top curve (red) 
is proportional to the out-of-phase susceptibility x"> while the lower one (blue) is 
proportional to the in-phase susceptibility x' ■ The data was divided by the mass of 
the sample. The peak in x" yields the critical temperature of the superconducting 
transition, while the amplitude of the diamagnetic signal, the value at which x' satur- 
ates towards low temperature, is, to first order, proportional to the superconducting 
volume fraction of the sample. The arrows indicate the direction of the temperature 
variation: it was first lowered from 1.6 K to .1 K, and then raised again to 1.6 K. 



to demagnetisation effects, this geometry made our calculation of the superconducting 
fraction difficult, but in order to obtain the order of magnitude of the content of 
Sr 2 Ru04, we assumed that the diamagnetic signal was proportional to the volume 
fraction to first order. 

The parameters of these experiments were the following. The magnetic field was 
set to zero, which meant that only a remanent field of less than 50 G was present. 
The temperature was swept from the 1-K pot temperature to the base temperature 
and back to the 1-K pot temperature of the ADR at a rate of 25 mK/min. Finally, 
the samples were excited with a modulation field of 0.428 G at a frequency of 71 Hz. 

In order to estimate an order of magnitude of the volume fraction from the amp- 
litude of the diamagnetic signal measured in the c-axis direction, some assumption 

21 The remanent field of the magnet was due to the use of field for the dHvA experiments that 
were also performed at each cooldown. 
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had to be made about the demagnetisation factor n. We supposed the following: the 
average demagnetisation factor in the c-axis direction is to first order roughly the 
same in all samples of interest or, in other words, that the shape of the inclusions 
had a distribution that was statistically the same. This is reasonable for samples 
of small Sr 2 Ru04 content, according to the polarised light experiments. The second 
assumption we used was that the demagnetisation factor in the afe-plane direction is 
n — 0, such that in that direction we measured the real susceptibility, and hence the 
real volume fraction. Then, measuring the AC susceptibility of only one sample in 
both directions should provide us with an approximate proportionality factor for all 
the other samples, and thus with an upper boundary value for the volume fraction of 
superconducting material. 

We measured the AC susceptibility in both directions of sample C698G (figure 
I2.14p . which yielded a ratio at base temperature of n + 1 = 8.85 for the c-axis value 
with respect to the a6-plane value. Taking the normalisation factor of section 12.3.31 
eq. 12.101 the values of AC field and frequency, and transforming the volume into a 
mass of Sr 3 Ru20y using the density of 6.117g/cm 3 , we obtained 

AV AV 

[n + l)m m 

which ranges between and -1. This corresponds effectively to a fraction of the real 
susceptibility of a sample entirely superconducting, and consequently to the volume 
fraction of Sr 2 Ru0 4 into Sr 3 Ru 2 07. Since a large error is associated with our assump- 
tions regarding the demagnetisation factor, these could be wrong by up to a factor 2 
or more, but we considered reasonable their order of magnitude. 

2.4.4 The de Haas van Alphen effect as a measure of disorder 

The best way to determine in a comparative way the disorder in large crystals of 
high quality is to perform a dHvA experiment on them. Using the AC susceptibility 
system described in section |2".3. 2\ we systematically measured dHvA on samples using 
the same parameters such that the amplitude of dHvA was comparable between runs. 
In order to be able to compare the quality of the various samples, we used the square 
root of the integral of the power spectrum, divided by the mass, as a measure of the 
amplitude, which according to eq. 11.181 depends exponentially on the inverse of the 
mean free path. We obtained values between zero and 100 nV/g. Note that these 
values were always measured with the c-axis parallel to the magnetic field. 

Intensity values of 10 nV/g were usually close to the noise levels of our suscept- 
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Figure 2.14: Susceptibility measurements as a function of temperature performed on 
sample C698G, of Sr3Ru207, with the excitation field aligned with the c-axis and the 
afr-plane of the crystal. The demagnetisation effect produced a larger difference in 
voltage when the field was aligned with the c-axis by a factor 8.85. 



ibility syste Effectively, as we always used a time constant for detection of one 
second, transformers with amplification factors of 100 and sample masses of about 
50 mg, this corresponded to intensities between .5 and 10 nV/g. When only noise was 
observed we quoted the intensity value with a less than symbol but still retained the 
integrated intensity in the usual frequency region. The parameters used for measure- 
ment were the following. We applied a modulation field of 0.428 G at a frequency of 
71 Hz. We used the field region between 6 and 7 T, where, for Sr3Ru20y, the oscil- 
lations have the highest amplitude, and varied the DC field at a rate of 20 mT/min. 
We verified carefully that these parameters did not cause eddy current heating, as 
was described in section 12.3.31 (see, in particular, figure l279|) . 



2.4.5 Magnetisation 

Measuring the magnetisation offered the easiest way, in our laboratory, to detect 
ferromagnetic signals in samples with extremely high sensitivity. We had the op- 
portunity to use a Quantum Design MPMS Superconducting Quantum Interference 
Device magnetometer (SQUID). This system can reach temperatures between room 

22 See section 1273.21 for details 



Chapter^ Experimental methods and analysis procedures 



74 



Magnetic field ( T ) 




2 4 6 8 10 12 14 



Frequency ( kT ) 

Figure 2.15: Typical dHvA measurement performed in the ADR. The top curves, 
with axes at the right and top, feature the oscillations in the measured voltage, which 
is proportional to the in-phase susceptibility. The lower curve is the absolute value 
of the Fourier spectrum, with axes at the left and bottom. The integral of the peak 
at 1.8 kT in the spectrum divided by the sample masses yields an intensity for the 
oscillations of 100 nV/g in this case, the mass being of 80 mg. 

temperature and 2 K, and fields of up to 5 T. It allows, among other things, to 
measure the DC magnetisation as a function of temperature or magnetic field. 

When measuring the magnetisation of our samples as a function of temperature, 
starting at room temperature and decreasing to 5 K, at small fields, we generally 
observed two super-linear rises of the total magnetic moment, at around 105 K and 
165 K (see figure |27L6| top right, solid line). The first is the Curie point of Sr 4 Ru30io 
[28] and the second, that of SrRuOs [26]. These were superimposed on the paramag- 
netic signal of Sr 3 Ru 2 7 , which features a peak at 16 K [21] (see figure 12.161 top 
right, dashed line). This moment being smallest at low fields, it often had lower mag- 
nitude than the spontaneous ferromagnetic moment of small inclusions SrRuOs and 
Sr 4 Ru 3 Oi . It was thus relatively easy in this way to detect their presence quantitat- 
ively, and for that we only needed to calibrate the amplitude of the two ferromagnetic 
signals into volume fractions of samples. Figure I2.16[ bottom right, shows the field 
cooled spontaneous magnetisation of sample C6102A, the one that exhibited the fer- 
romagnetic signals of highest amplitude, and which was used for the calibration of 
the volume fractions. 

To further demonstrate the presence of ferromagnetic materials in our samples, 
we performed hysteresis loops. The left side of figure 12.161 shows the one performed 
on sample C6102A, for magnetic fields aligned with the c-axis (top) and the ab- 
plane (bottom). In the first case, the ferromagnetic signal for both S^R^Oy and 
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Figure 2.16: Magnetisation measurements performed with the SQUID magnetometer. 
Left Hysteresis loops measured on sample C6102A, which featured large amounts of 
SrRuC>3 and S^RusOio- The top graph is for a measurement with the magnetic field 
along the c-axis, and the bottom one for the afe-plane. Arrows indicate regions of 
hysteresis. Top right Typical SQUID mass susceptibility measurements as a function 
of temperature, with the sample field cooled at a magnetic field of 20 G. The dashed 
curve is for a very pure sample of Sr 3 Ru207, while the solid one was measured on 
a sample featuring small amounts of SrRuC>3 and Sr 4 Ru30io- Bottom right Field 
cooled mass susceptibility as a function of temperature for sample C6102A, with a 
magnetic field of 20 G. 



Chapter Experimental methods and analysis procedures 



76 



3 

E 

<x> 



0.2 



0.15 



0.1 



r 0.05 



c 

CD 

E 
o 



-0.05 

c 

(0 



-0.1 



-0.15 



-0.2 



□ SrRu0 3 
„ Sr Ru 

o ^4 3 10 



o o 



o □ □ 

o 



B 
B 

n » ° 



-1000 -500 500 

Magnetic Field ( G ) 



1000 



Figure 2.17: Values of the extrapolation to K of the field cooled spontaneous mag- 
netisation performed at various magnetic fields on sample C6102A. The squares are 
for SrRuOs and the circles for Sr 4 Ru30io. 



SrRu0 3 were indistinguishable and saturated below IT, indicating that the easy axis 
of SrRu0 3 was probably not in the same direction as that of Sr 4 Ru 3 Oio, but lay in 
the direction of the bilayer of Sr 3 Ru207. In the second case, it was Sr 4 Ru30io, as we 
know from the work of Cao and co-workers [28], that featured hysteresis loops near 
2.5 T, making its identification easy. We saw three sharp rises in magnetic moment, 
at 0, 2.5 and above 4 T. They correspond, in order, to SrRu0 3 [26J, Sr 4 Ru 3 Oio and 
Sr 3 Ru 2 07, [27] , The first two produced hysteresis, but not the third one, as one 
expects for two ferromagnetic and a paramagnetic materials. 

In order to calibrate the measurement of the field cooled spontaneous magnetisa- 
tion as a function of temperature for all samples, we used sample C6102A, for which we 
could neglect the signal of Sr 3 Ru2C>7. We used a simple extrapolation to zero temper- 
ature of the total moment of each ferromagnetic phase. Figure 12 . 1 71 shows how these 
moments evolved with the magnetic field at which the samples were cooled, for both 
SrRu0 3 and Sr 4 Ru 3 Oio. One sees that the data are symmetric, indicating that the 
remanent field of the SQUID magnet is very small, and both data sets saturated above 
0.1 T. Using saturated moments of 0.112 emu (SrRu0 3 ) and 0.180 emu (Sr 4 Ru 3 Oi ), 
and those given by Cao ([SSI ESJ) of respectively approximately 1.0 /xg/Ru and 1.1 
/jb/Ru, which correspond to 23.6 emu/g and 22.7 emu/g, we obtained total masses 
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of 4.7 mg of SrRu0 3 and 7.9 mg of Sr 4 Ru 3 Oi . These corresponded to 9% and 16% 
of the total sample mass. 

From the lower right graph of figure \2A6\ we extrapolated the moments to OK, and 
obtained 0.049 emu for SrRu0 3 and 0.075 emu for Sr 4 Ru 3 Oi . With the masses pre- 
viously calculated, we obtained the values of magnetic moment at 20 G of 12.9 emu/g 
and 10.3 emu/g. These are the numbers we used to calibrate the volume fractions of 
SrRu0 3 and Sr 4 Ru 3 Oi in all the other samples. 



2.4.6 Complete quantitative results of the search 

In this section, we will present the complete results of the sorting procedure for 
23 new samples of Sr 3 Ru 2 7 . Table 12.11 presents all quantities measured using the 
characterisation methods described in this section: powder X-ray diffraction, resid- 
ual resistivity, dHvA, zero field AC susceptibility and low field magnetisation. The 
same measurements were performed on the sample that was used in preliminary ex- 
periment^, shown here in the first row of data named "Reference", for comparison 
purposes. 

Using the various calculations described in the previous sections, we estimated 
volume fractions of Sr 2 Ru0 4 , SrRu0 3 and Sr 4 Ru 3 Oi for each sample. Table 12.21 
summarises these values, along with the same amplitude values for dHvA as the pre- 
vious table, and the average residual resistivities. Reading carefully this table reveals 
the strategy to employ when choosing samples for extensive dHvA experiments, which 
we describe here. 

The very first column of interest is that of the amplitude of the quantum oscilla- 
tions at 1.8 kT. The reference sample provided us with a value for comparison: all 
samples with a value higher than 15.5 nV/g were worth using in dHvA experiments, 
but six of them had exceptional values, which were C697C, C698A, C698B, C698F, 
C698H and C698I. As we know from experience with crystal growth, it was hardly 
surprising that most of them originated from the same crystal growth experiment, of 
batch name C698. According to these data, sample C698I featured by far the highest 
amplitude, 100 nV/g, over 6.5 times that of the reference sample. 

The next column of interest is that of the mass fraction of Sr 2 Ru0 4 calculated 
from the amplitude of the superconducting diamagnetic signal. We observed that 
even within a single crystal batch, for instance C698, the amount of Sr 2 Ru0 4 varied 
by orders of magnitude. Among our exceptional samples, it was sample C698A that 
featured the smallest amount, 50 times less than the reference sample. It was also 
interesting to find that the highest value belonged to the reference sample, meaning 



23 The results of these experiments are not shown in this thesis, since equivalent higher quality 
measurements are presented in chapter [31 
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Table 2.1: Complete qualitative results of all characterisation methods for all the 
samples that were analysed. The second and third columns give the integrated dHvA 
amplitude for the 1.8 kT peak, related to Sr 3 Ru207, and the 3.0 kT peak, related 
to Sr 2 Ru04, normalised by the mass of the samples. The fifth column is the ratio 
R of the powder diffraction peaks at 8.6°, associated with Sr 3 Ru 2 7 , and 13.9°, of 
Sr 2 Ru0 4 . The next two show the residual resistivities measured on each side of the 
resistivity crystals. Next are shown the columns with data from AC susceptibility, the 
superconducting critical temperatures and the amplitude of the diamagnetic signal. 
Finally, in the last two columns we present the amplitude of the ferromagnetic signals 
of SrRu0 3 and Sr 4 Ru 3 Oi measured with the SQUID, extrapolated to absolute zero. 
The first row of data belongs to the sample that was used in preliminary experiments, 
and was the standard with which to compare the new samples. The empty elements 
of the table indicate measurements that were impossible or failed for various technical 
reasons. 
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Sample Name 


Mass 


Disorder 


Impurity fractional content 






dHvA Resistivity 


Sr 2 Ru0 4 


SrRuO, 

3 


fi a 

Sr 4 Ru 3 O 10 




mg 


nV/a 

1 1 v/ y 


ii Qcm 


% 


% 


% 


Reference 


51.7 


15 5 


1 3 


8 


0.009 


0.03 


C697A 


38 


< 5 5 

Jl mJ 


85 


2 


2 


5 


C697B 


55 


16.7 


0.25 


6 


0.03 


0.3 


C697C 


57 


56.1 


0.63 


6 


0.02 


0.4 


C6102A 


49.9 


< 6 




0.06 


9 


13 
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51.9 


16.5 




0.5 


0.03 


0.2 
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< 4 


0.6 
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5 
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< 11 


1 i 
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33.6 


24 
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0.7 
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64 


< 5 
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7 


9 
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60 


66 




0.2 


0.001 


0.002 


C698B 


42.4 


56 


0.98 


0.01 


0.002 


0.05 


C698C 


58.2 


19 


0.81 


0.07 


0.007 


0.04 
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53.9 


22 




4 


0.4 


0.01 
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47.6 


25 




4 


0.04 


0.01 


C698F 


62.9 


63.6 


0,76 


3 
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27.6 


0.53 
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0.01 


0.01 
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55.4 


52.3 


0.72 
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31 


0.68 


5 


0.002 
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0.45 
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0.03 


0.1 


C698J 
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0.48 


0.6 


0.008 
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31 


0.51 
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0.03 


0.05 
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49.7 


27 


0.59 


16 
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43.9 


< 7 




12 







Table 2.2: Fractional mass of the impurity phases Sr 2 Ru04, SrRuC>3 and Sr 4 Ru30io in 
ppm, for all the measured samples. These are quoted along with the amplitude values 
of dHvA and the average residual resistivity, for convenience. Five samples feature 
remarkable dHvA properties. Among them, one possesses by far a much higher dHvA 
signal, C698I, while one contains negligible amounts of the impurity phases, C698A. 
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that in this respect, all the new samples were cleaner. Moreover, the sample which 
was best for dHvA, C698I, featured a lot more Sr2Ru04 than C698A, by a factor 38. 

The third set of values to examine was the mass fraction of SrRu0 3 and Sr 4 Ru 3 Oio- 
The two values, for each sample, seemed correlated. This is not surprising, since we 
intuitively know that these inclusions correspond to stacking faults. The n > 2 mem- 
bers of the Ruddlesden-Popper series possess physical properties that resemble more 
those of the n = oo member as n increases, and it is probable that slight deviations 
in growth parameters can produce a whole variety of them, and members with very 
large n are essentially indistinguishable from n = oo. However, the fractional val- 
ues for SrRu0 3 and Sr 4 Ru 3 Oi were visibly not correlated with those for Sr 2 Ru0 4 . 
The sample with the smallest amounts of SrRu0 3 and Sr 4 Ru 3 Oi was C698H, for 
which the ferromagnetic signals were not at all observed. In this regard, the reference 
sample contained very small amounts of these, while for sample C698I the fractions 
are 4 (SrRu0 3 ) and 39 (Sr 4 Ru 3 Oi ) times larger. Sample C698A also contained very 
small amounts of ferromagnetic inclusions, being one of the samples with the best 
phase purity. 

Moreover, the residual resistivities were surprisingly not correlated with dHvA 
amplitude values. This effect can only be explained by the fact that the resistivity 
crystals were not representative of large dHvA samples, much less of whole crystal 
rods, at these levels of purity. One should only rely on residual resistivity to infer 
disorder for bulk crystals only in a larger scale of impurity content, where quantum 
oscillation experiments are not possible, or to a very small spatial scale within a 
crystal. 

Before we finish this section, we make a few more observations. Examining the 
data from powder X-ray diffraction, in table I2.1[ one finds that it also possessed 
no apparent correlation with the amplitude of the diamagnetic signal produced by 
Sr 2 Ru0 4 . These data were not representative of large crystals, except in the limit of 
large volume fractions. A signal was almost always observed at 13.9°, even in cases 
where very low amounts were present in the neighbouring dHvA sample. Also, in some 
instances, the quantum oscillation signal of the 3.0 kT orbit of Sr 2 Ru0 4 was observed, 
and its amplitude was not correlated with that of the diamagnetic signal. Since in 
such cases, the quality of the inclusions was very good, the disorder in Sr 2 Ru0 4 
intergrowths was not correlated with the fractional amount. However, the correlation 
could probably be made between the amplitude of the signal of the 3.0 kT dHvA 
frequency and the superconducting critical temperature, since these measurements 
were performed on the same crystal and, as was demonstrated by Mackenzie et al. 
[oT] . a correlation exists between disorder and the superconducting T c . 

In conclusion, we summarise by describing the strategy we have used for choos- 
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ing crystals for the extensive dHvA experiments on S^R^Oy that is at the core of 
this thesis project. In order to maximise the dHvA signal, we chose the sample with 
highest dHvA amplitude at 1.8 kT, C698I. But we know, from preliminary experi- 
ments and our sorting procedure, that it is possible that with higher resolution, dHvA 
peaks produced by Sr 2 Ru04, SrRuC>3 or Sr 4 Ru30io are detected using this sample 
since the volume fractions of these impurity phases are not negligible. Consequently, 
we chose to use as well a sample with a slightly higher disorder level but lower amounts 
of impurity phases, C698A, in order to make rule out interpretation errors due dHvA 
frequencies of other phases than Sr 3 Ru207. 

2.5 Laue orientation of the samples 

Before placing samples into pick-up coils of a dHvA system, in order to measure and 
interpret the angular dependence of the quantum oscillations, one is required to find 
the direction of the crystallographic axes using Laue X-ray diffraction. For Sr3Ru207, 
the c-axis is easy to determine, as the material readily cleaves along the ab-plane, but 
finding the orientation of the a and b directions is impossible by eye. One is required 
to use the standard procedure with Laue X-ray diffraction, which is presented here, 
with the resulting photos. 

We used the Laue X-ray source and camera situated in Cambridge, a model that 
employs polaroid photographs. We placed the samples such that their c-axis was 
parallel to the X-ray beam, and the photographic plate was put between the source 
and the sample, the sensitive side towards the sample. The beam then passed through 
the centre of the plate and was reflected by the sample back to the photographic plate. 
We obtained complex dot patterns along straight lines, shown in figure 12.18} top left 
and right, for samples C698I and C698A, respectively. 

In order to find out which of those lines of spots corresponded to which crys- 
tallographic direction, we simulated such patterns with a computer program, Orient 
Express, using the space group of I4 mm rH with the right .attics palters of 
a = b = 3.88A and c = 20. 7 A. We obtained the plot on the bottom right of figure 
12.181 In this simulation, one recognises the [100] from the [110] directions by using 
the sets of spots between the vertical or horizontal lines and the diagonal ones, which 
lie slightly closer in angle to the diagonals than to the vertical or horizontal lines. We 
chose to place sample C698I in a pick-up coil of the system described in section 12.1.31 
such that the magnetic field rotates from [001] towards [110], and towards [100] for 
sample C698A. 

24 The rotation of the RuO octahedra does not affect the Laue pattern in a significant way, and it 
looks qualitatively the same as that of Sr2Ru04. 
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Figure 2.18: Top left Scanned Laue photograph taken on sample C698I. Top right 
Scanned Laue X-ray photograph taken on sample C698A. Bottom left Laue pho- 
tograph taken on sample C698K by A. Rost with the system at the University of 
Edinburgh. Bottom right Simulation of the Laue pattern using the software Orient 
Express, where the a and b directions are along the figure edges, and equivalent. 
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The orientation procedure of samples C698I and C698A was verified independently 
by A. Rost, with a different Laue system located at the University of Edinburgh, using 
sample C698K, in preparation for specific heat experiments, shown in the lower left 
panel of figure 12.181 Sample C698K was a single crystal neighbour of C698I, and 
consequently possessed the same orientation. The sample was rotated until the Laue 
picture presented crystallographic directions aligned with the photo edges, and yielded 
the same rotation angle. 

2.6 Avoiding eddy current heating in dHvA ex- 
periments 

The effect of eddy current heating is to increase the sample temperature compared 
to that of the thermometer, which in a dilution refrigerator, is located outside of the 
magnetic field region, at the mixing chamber. In dHvA experiments, even for small 
eddy current heating, the systematic error on the thermometry produces dramatic 
effects on the calculation of the quasiparticle effective masses during LK fits. This 
can be explained by considering that the effective mass is inversely proportional to 
the width at half maximum of the LK function (see figure [TTTOl) . A positive systematic 
error on the measured sample temperature always leads the non-linear fit of the LK 
function to underestimate the quasiparticle mass. Moreover, for a same heat input, 
this effect is the more pronounced the higher the mass is. Overall, this may have the 
effect of suppressing, partly or completely, an enhancement of the quasiparticle mass 
as a function of magnetic field, and serious measures had to be taken in this project 
to ensure that such heating was not significant during dHvA experiments. 

The source of heat in this phenomenon is produced proportionally to the rate of 
change of magnetic field, and the modulation field normally possesses the largest time 
derivative. For a conducting cylinder of volume V and resistivity p in a varying field, 
the power dissipation is ^ 



Therefore, one is required to use the lowest possible excitation field, which also reduces 
the amplitude of the signal, and a compromise has to be found where the temperature 
difference between the thermometer and the sample is negligible. The quasiparticle 
mass of one dHvA frequency extracted with a non-linear LK fit is a good measure 
of sample heating. The mass saturates to a maximum value as the excitation field 
is decreased, and for high eddy current heating, significant deviations from LK are 

25 see the book by Pobell, [ST], p. 177 



P 



eddy 




(2.12) 
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Figure 2.19: Data for the temperature dependence of the oscillation amplitude of 
the 1.8 kT peak between 6 and 7 T for two values of the product of modulation 
field and frequency: ujHac = 676.5 G-Hz (left) and 12.5 G-Hz (right). The data 
was measured with different but equivalent experimental apparatuses, so the absolute 
values of amplitude cannot be compared directly. Solid lines correspond to LK fits 
and the results are given in the legends. 

measured. 

The guidelines we used stem from the following considerations. We showed in 
section 12.1.21 (see eq. 12.121 and 12 .5 j) that, in the limit where A < 2, 

1. Eddy current heating is proportional to H 2 oc u: 2 H\ c , 

2. The signal in the first harmonic is proportional to uHac, 

3. The signal in the second harmonic is proportional to uH\ c . 

We conclude that increasing by the same factor the AC field or its frequency produces 
the same increase in eddy current heating, as well as the same increase in amplitude 
of the first harmonic, and therefore increasing one of these parameters while keeping 
the product ujHac constant results in constant eddy currents and first harmonic 
amplitude. However, it is not the case with the second harmonic, where reducing the 
frequency and increasing the modulation field increases the signal while leaving the 
heating constant. Thus, when using the second harmonic, one should use the lowest 
frequency possible. 

Figure 12.191 presents two examples of the temperature dependence of the dHvA 
amplitude and LK fits at two extremes, the first with a product of modulation field 
times frequency high enough that eddy current heating is significant, left, and the 
second with such a product 54 times smaller, producing eddy current heating nearly 
3000 times smaller. We found in the first case that we could not properly fit the 
LK relation to the data. The adjusted curve overestimates intermediate temperature 
points. Since, with eddy current heating present, the error on temperature values 
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decreases as temperature increases, deviations to LK decrease at high temperature 
points. By excluding low temperature points, we found that better fits could be 
obtained, and that using only data above 149 mK produced the same LK fit result 
as in the low modulation field case. The quasiparticle mass obtained increased as 
we removed data points saturated to the value measured with low heating, and in 
doing so the fit quality improved at intermediate points. Therefore, for extended 
dHvA experiments in a material that respect LK temperature dependence, one should 
use the highest modulation field at which the LK distribution is not significantly 
disturbed. 



Chapter 3 



The de Haas van Alphen 
experiment on S^R^Oy 

This chapter presents the complete experimental data obtained from the de Haas 
van Alphen experiments that we performed on S^R^Oy. These measurements were 
carried out in three sets, performed in 2005, 2007 and 2008, the first two in St Andrews 
and the third in Cambridge. The St Andrews experiments did not yield a complete 
set of information for the completion of this project, but mainly served to generate 
the motivation for a collaboration with the dHvA group of Michael Sutherland in 
Cambridge. The combination of the high quality samples that we provided and the 
performance of the dHvA apparatus in Cambridge, described in section I2.1.3[ has 
been the key to the success of this project. 

dHvA has been measured in the Cavendish Laboratory in Cambridge for decades, 
mainly by the late professor David Shoenberg, who was also the author of the main 
reference on the subject, still used today [35] . and responsible for the determination 
of the Fermi surfaces of a large number of elements, notably the noble metals. The 
current Cambridge dHvA system was developed mainly by Stephen Julian, Gilbert 
Lonzarich and Christoph Bergemann, and complete Fermi surfaces of many com- 
plex systems were determined with it, notably those of Sr 2 Ru04, CaV03, URu2Si2, 
CeRu 2 Si 2 , UPt 3 and recently Ag 5 Pb 2 6 . Charge of this apparatus was recently taken 
by Michael Sutherland, under the direction of whom our experiments have taken place. 
The work was performed by myself in close collaboration with Cambridge postgradu- 
ate students Swee K. Goh and Eoin C. T. O'Farrell. All subsequent data processing 
and analysis were performed by myself only. 

All the samples that were used were grown in St Andrews by Robin Perry and 
thoroughly characterised by myself. The 2005 series of experiments were performed 
using the sample labelled "Reference" , in comparison to which we strove to improve 
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the purity of crystal^. The 2007 and 2008 series of experiments were carried out using 
simultaneously two samples of the highest quality recorded by our sorting procedure, 
labelled C698I and C698A, which featured very low disorder and high phase purity. 
Finally, since most measurements performed in St Andrews were repeated in Cam- 
bridge, where the signal quality was superior, only data from the 2008 experiment 
series in Cambridge are presented, except where stated otherwise. 

The first section compares detection harmonics, followed by the lengthy way in 
which we determined the modulation field to use in order to avoid eddy current 
heating. This involved measuring the temperature dependence of the dHvA signal in 
a certain field region for many different modulation fields. The next section shows 
the general description of the dHvA data, spectra and quasiparticle masses away 
from the metamagnetic transition. We then present the largest data set, the angular 
dependence of dHvA. Next, we reinvestigate the temperature dependence of dHvA in 
Sr 3 Ru 2 07 and the field dependence of the quasiparticle masses, in which we contradict 
one of the conclusions of the previous work by Borzi [40J. This is finally followed by 
another very successful experiment, the measurement of dHvA inside the nematic 
phase. 

3.1 Determination of experimental parameters 

This section is dedicated to the method used to select the optimal parameters with 
which we performed the dHvA measurements. There were many ways in which we 
could perform this experiment, and we desired to tailor it to our specific needs. 
In principle, one can use various values of modulation field and different detection 
harmonics, which enhance different aspects of the data, as was seen in section 12.1.21 
However, the modulation field can lead to significant heating through the process of 
eddy currents (see section [276]) . and this aspect also needed attention. 

The method with which we determined the amount of heat generated by eddy cur- 
rent heating relied on the principle that a systematic difference between the measured 
amplitude of dHvA for one orbit and that predicted by theory appears when strong 
heating occurs at the sample such that there is a temperature difference between the 
sample and the thermometer. It stems from the fact that, by the LK relation (eq. 
11.15] p. 124"]) . the dHvA amplitude decreases monotonically as temperature increases, 
and when heating occurs, the measured value is lower than the prediction. There are 
two ways by which such a difference can be measured. The first consists in measuring 
dHvA as a function of modulation field and compare the measured amplitude to the 
Bessel function it is expected to follow (eq. 12.4] p. 138]) . The second is to look at the 

1 See section 12711 for details. 
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Figure 3.1: dHvA measured using first (top) and second (bottom) detection harmonics. 
The data were measured on the sample labelled "Reference" in section 12.4} using a 
modulation field of 16 G at a frequency of 41 Hz. The total amplification factor was 
of 100, and the plots were not normalised. 



temperature dependence of dHvA for many modulation fields and compare with the 
LK relation. 

We first show why second harmonic detection is appropriate for measuring dHvA 
in Sr 3 Ru20"7, through comparison of the first and second harmonics. Then, we present 
the modulation field dependence of the second harmonic amplitude, which we also 
used to calibrate accurately the magnitude of the modulation field. Finally, we show 
which modulation field should be used in order to avoid eddy current heating. 



3.1.1 Comparison of detection harmonics 

We showed in section 12.1.21 that an important reduction of non-oscillatory magnetic 
background in AC susceptibility measurements can be obtained by using second har- 
monic detection. It is particularly the case in S^R^O^, where the magnetic signal 
produced by the metamagnetic transition is of several orders of magnitude larger 
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than the oscillations. These measures led to significant improvements in sensitivity 
compared to the experiments of Borzi [40J. 

Figure 13.11 presents a comparison of first and second detection harmonics when 
used with S^R^O^, where the top plot refers to the first and the bottom the second. 
The data was measured using the St Andrews system, with sample labelled "Refer- 
ence" in section 12.41 during preliminary experiments. The first harmonic corresponds 
to the in-phase component of the AC-susceptibility, the first field derivative of the 
magnetisation of Sr 3 Ru 2 7 and features three peaks. Those correspond to a crossover 
at 7.5 T, and two first order metamagnetic transitions at 7.9 and 8.1 T, and are the 
same as those in the data published by Grigera [TT]. With such a signal, setting the 
LIA sensitivity high enough that the oscillations away from the peaks are measured 
without digitisation results in the saturation of the LIA input in the peak region. 

The bottom plot of figure 13.11 shows the second harmonic, which corresponds to 
the out-of-phase component of the second derivative of the magnetisation. One can 
observe that the magnetic background is reduced to very sharp asymmetric peaks 
at 7.5, 7.9 and 8.1 T of smaller magnitude compared to the oscillations, and full 
sensitivity can be achieved on the LIA using a simple offset. With such data sets, 
dHvA can be analysed through the complete field range available. We conclude from 
this that in order to analyse properly dHvA near the metamagnetic transition, second 
harmonic detection was the appropriate method to use. 



3.1.2 Modulation field dependence of dHvA 

Figure 13.21 shows the modulation field dependence of the amplitude of the 1.8 kT 
frequency between 6 and 7 T, measured on both samples C698I (squares) and C698A 
(circles), using a modulation frequency of 7.9 Hz. The solid lines correspond to 
non-linear fits of Bessel functions. As we will see in section I3.3[ the alignment of 
the crystallographic c-axis of these was not identical, and a difference in field angle 
was present between them. The measurement was performed with sample C698I 
positioned with its c-axis parallel to the magnetic field, and simultaneously with 
C698A at an angle of 10°. 

Bessel functions were fitted to the data with scaling parameters a and b, using the 
form 

/(A) = aJ k (b\), k = 2, A= 2 ^f^ AC , (3.1) 

DC 

where F is the dHvA frequency and Hac the rms modulation fielcjf). When using the 
value of the voltage over a resistance of 1.7 Q in series with the modulation coil as a 

2 Which requires a factor y2 in order to be consistent with the calculation of section 12.1.21 (eq. 
12. 3p . in which Hac refers to half the peak-to-peak amplitude. 
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6 to 7 T, 7.8 Hz, T = 29± 5 mK 
9i ■ ■ ■ ■ 




Modulation Field ( G rms ) 

Figure 3.2: Amplitude of the dHvA signal for the 1.8 kT frequency between 6 and 
7 T, as a function of the rms amplitude of the modulation field at 7.9 Hz, compared 
with the predicted Bessel functions. The squares correspond to sample C698I, the 
circles to sample C698A, and the solid lines to fits of Bessel functions. 

variable, we obtained the values for b of 57.5 and 60.7 G/V using the coils containing 
C698I and C698A, respectively, corresponding to approximately 98 and 103 G/A for 
the modulation coil. Assuming it is possible that each pick-up coil pair experienced 
a slightly different modulation field amplitude due to field inhomogeneity in space, 
this allowed us to scale properly the abscissa of figure 13.21 for both data sets. 

It is not obvious, however whether there are deviations in the data due to eddy 
current heating. The fit reproduces extremely well the data below 60 G. There is 
a difference between the fit and the data near the peak of the function, where the 
fit overestimates the data. As we move to higher modulation, the fit joins the data 
again, suggesting that this is not the result of eddy current heating. Had it been the 
case, we would have expected the difference between the fit and the data to increase 
indefinitely with modulation field. Since there were no clear signs of eddy current 
heating in the modulation field dependence of the amplitude of the oscillations, we 
next looked at the temperature dependence and the fits of the LK relation. 
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3.1.3 Analysis of eddy current heating 

The ultimate test for eddy current heating lies with the non-linear fit of the LK func- 
tion. We expect the temperature dependence of the amplitude of dHvA in Sr3Ru207, 
away from the metamagnetic transition, to follow the LK function. Since the specific 
heat of most materials reduces with decreasing temperature, a constant heat input 
at the sample produces an error in temperature that becomes larger the lower the 
temperature. We conclude that any temperature gradient between the sample and 
the thermometer situated at the cooling apparatus will result in the low temperature 
part of the data saturating to a value lower than it would normally reach in normal 
conditions. This subtle change can make the data appear obviously non-LK, but it 
could also appear like an LK function with a lower quasiparticle mass than reality. 

We measured the temperature dependence of dHvA between 6 and 7 T, at seven 
different modulation fields between 5 and 80.4 G at 7.9 Hz, for both the 1.8 and 
the 4.2 kT peaks. Figure 13.31 shows the data for both frequencies, both samples for 
the modulation values of 5, 10.1, 20.1, 30.2 and 80.4 G. We observed a systematic 
deviation at low temperatures, which made the LK fits impossible when using all 
temperature points. By looking at the modulation field dependence of this deviation, 
we found there was none and that the deviation was identical at all modulation field 
values, which span over more than an order of magnitude. We concluded that this 
deviation was not related to the modulation field. By excluding temperature points 
below 90 mK, we found that the LK function fits the data well up to 80.4 G, where 
signs of eddy currents start to appeajf). The fits are shown in solid lines in figure 13.31 
and the results are presented in the legends. Slight deviations in the low temperature 
points generally force the LK fit to a lower mass, and one usually observes a deviation 
between fit and data at intermediate temperatures (as indicated by a black arrow in 
the lower left plot, figure 13. 3p . However, the value of the extracted mass was the 
best indication of heating: we observed a deviation from the normal on-axis value of 
7.15 ± 0.10m e for the 1.8 kT peak, and 11.0 ± 0.2m e for the 4.2 kT peak, at 30 G 
and above @. We concluded that eddy current heating is only significant at such 
a modulation field or higher, and that the appropriate value to use safely lies near 
20.1 G, at 7.9 Hz. 

We furthermore verified whether the deviation had its origin in the proximity 
to the metamagnetic transition, where effects like spin-dependent masses or non-LK 
behaviour could potentially be observed. We measured the temperature dependence 

3 Note that sample C698A lies with an angle of 10° between the field and its c-axis, leading to 
slightly higher masses and frequencies since, as seen from eans ll.l^l and flTTOl the extremal area and 
the mass, in 2D systems, scale with the inverse of cos 8 

4 The same can be said of the values measured on C698A, slightly off-axis, which produced higher 
mass values, but the same deviation. 
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Figure 3.3: Modulation field study of eddy current heating. The temperature de- 
pendence of the amplitude of the 1.8 and 4.2 kT peaks were measured for different 
modulation fields, from left to right, top to bottom: 5, 10.1, 20.1, 30.2, 80.4 G and 
20.1 again, at 7.9 Hz. The data were measured on samples C698I and C698A, where 
C698I was aligned with its c-axis parallel to the magnetic field, and C698A lied with 
an angle of 10°. All Fourier transforms were performed between 6 and 7 T, except 
for the bottom right graph, which was measured between 4.15 and 4.65 T. Solid lines 
correspond to LK fits, and the resulting masses are quoted in the legend. Error values 
are those given by the fit procedure. 
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of dHvA in another field region of similar inverse width, 4.15 to 4.65 T, much further 
away from the metamagnetic transition (see lower right graph of figure 13. 3j) . The 
same unusual behaviour was observed at low temperatures, excluding this possibility. 
Moreover, the calibration of the thermometer was later verified by comparison with 
another new RuO thermometer that was installed on the mixing chamber of the 
refrigerator, which revealed deviations too small to explain this problem. We conclude 
that if a temperature gradient was present in our system, it was not due to eddy 
current heating, and it is still not resolved at this point whether the deviation is part 
of the physics taking place in SrsR^CV or if it was an experimental problem. As 
a result, all LK fits performed in this thesis were done using data from 90 mK and 
above. 

3.2 dHvA oscillations and spectra 

We present in this section information that was extracted from the dHvA data for the 
c-axis direction. We first show raw data of second harmonic dHvA for both samples. 
We then reveal spectra that were calculated over large field ranges both in the low and 
high field sides of the metamagnetic transition, and provide a table of the extracted 
frequencies and quasiparticle masses. These are average values, mainly of interest for 
their comparison with, for instance, ARPES data (introduced in section |4~T]) . Finally, 
we show plots of the field dependence of the dHvA frequencies, which will be used 
later in section 14.2.21 in order to demonstrate that quasiparticle transfers occur in 
Sr 3 Ru207 at the metamagnetic transition. 

3.2.1 dHvA oscillations 

The oscillations measured in Cambridge were of exceptional quality compared to any 
of our previous experiments. The samples, C698I and C698A, were measured sim- 
ultaneously with second harmonic detection only, using the AC susceptibility probe 
described in section [2.1.31 with a modulation field of 20.1 G at 7.9 Hz, and a sweep 
rate of approximately 0.04 T/min above 10 T, and 0.02 T/min below, as described 
in section 12.1.31 

Temperatures of 25 mK or above were reached, and angles between -10° and 
54° (see appendix [E] for details) were probed. We found a discrepancy between the 
DC field positions of the metamagnetic transitions as measured in these experiments 
compared to previous ones performed in St Andrews, of approximately 3%. Since the 
values measured in St Andrews were previously confirmed by two other laboratories 
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Figure 3.4: Raw second harmonic data measured at 25 mK, using a modulation of 
20.1 G at 7.9 Hz, for both samples in the c-axis direction. The total amplification 
factor is of 10 5 . 



elsewhere in the worldj, we corrected all Cambridge data sets by the appropriate 
scaling factoid- 
Figure 13.41 shows typical dHvA data at 25 mK for both samples, where the amp- 
litude is approximately twice higher with C698I. These were not measured simultan- 
eously, but at two different rotator positions, such that the two crystals were aligned 
with their c-axis parallel to the magnetic field. Since the rotation study was per- 
formed before the mass analysis, we knew precisely which rotator positions aligned 
each sample. 



3.2.2 dHvA spectra and quasiparticle masses at c-axis 

The dHvA spectra away from the metamagnetic transition were obtained by using 
Fourier transforms applied to the data within the intervals of 5 to 6.5 and 10 to 18 T, 

5 At the Max Planck Institute in Dresden, and at the University of Kyoto 

6 The conversion number used was of 0.9743, the field values of Cambridge being too high. 
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Figure 3.5: Fourier transforms of the dHvA data for both samples C698I and C698A, 
in the low field side (top parts of the graphes), using data between 5 and 6.5 T, and 
the high field side (bottom parts of the graphs), using data between 10 and 18 T. 
The peak at 1.8 kT was cut off for clarity at 350 /iV/T and 120 /xV/T for C698I and 
C698A respectively. 
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5 to 6.5 T 


10 to 18 T 


Peak (kT) 


m*/m 


Peak (kT) 


m* /m 


0.15 


5.6 ±0.7 


0.14 


4.9 ±0.1 


0.43 


8.4 ±0.3 


1.01 


10.9 ±0.2 


0.82 


7.7 ±0.7 


1.12 


10.57 ±0.6 


0.91 


7.7 ±0.8 


1.20 


11.6 ±0.2 


0.98 


8.0 ±0.3 


1.30 


9.3 ±0.3 


1.78 


6.6 ±0.4 


1.35 


10.6 ±0.8 


4.13 


10.1 ±0.1 


1.45 to 1.73 


7.3 ±0.2 


4.30 


10.1 ±0.1 


1.75 


8.4 ±0.5 


4.40 




1.88 


12.2 ±0.1 






2.00 


14.4 ±0.1 






3.14 


13.3 ±0.2 






3.26 


11.7 ±1.5 






3.36 


12.5 ±0.1 






3.51 


14.9 ±0.1 






3.65 


12.6 ±0.1 






3.74 to 4.16 


11.3 ±0.3 






4.16 to 4.32 


13.9 ±0.2 






4.37 


13.0 ±0.1 






4.57 


15.4 ±0.2 






7.94 


17.2 ±0.3 



Table 3.1: Table of the observed frequencies in the Fourier transforms of the dHvA 
data from sample C698I, along with associated quasiparticle masses. The Fourier 
transforms were taken from 5 to 6.5 T, for the low field side, and from 10 to 18 T 
for the high field side. Some peaks in the high field side were difficult to integrate 
individually, where in such cases larger frequency intervals were used in the integrals 
for the non-linear LK fits. Such intervals are indicated. 
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Figure 3.6: Temperature dependence of the amplitude of the oscillations for the five 
main frequencies present in the low field side, along with LK fits, shown in solid lines. 
Fourier transforms were taken from 5 to 6.5 T in all cases except for F = 0.15 kT, of 
which the Fourier transform was taken from 6 to 7 T. 
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Figure 3.7: Field dependence of the dHvA spectra for samples C698I (top) and C698A 
(bottom), at 25 mK. The intensity of the data was cut off at 25 (top) and 8 uY/yT 
(bottom) in order to produce contrast for the low amplitude details. 
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which correspond to similar inverse field intervals. Since, as we know from the work 
of Borzi [40J, the dHvA frequencies vary close to 8 T, we could not include that region 
in the calculation. Moreover, field ranges below 5 T do not improve the clarity of the 
spectra. 

Figure 13.51 shows Fourier spectra for both samples calculated using these field 
regions. In the low field side, we obtained five peaks or groups of peaks, positionned 
at 0.15, 0.43, 0.9, 1.8 and 4.1 kT. These are essentially the same as those measured 
by Borzi, but we observed that the peak at 4.2 kT is split into three and the one 
near 0.9 kT into fouiQ. Moreover, both samples produced the same results, with 
a resolution slightly lower with C698A. In the high field side, the improved data 
quality revealed more complexity than in the work of Borzi. We observed four groups 
of peaks, each of which appears to correspond to a superposition of a great number 
of frequencies and lie near 1.6, 3.2, 4.0 and 7.9 kT. 

Table EH] presents two lists of frequencies and corresponding masses for both field 
sides of the metamagnetic transition in sample C698I. An identical set of values was 
also calculated for sample C698A, which turned out very similar, and is not shown. 
The table moreover lists the value of the quasiparticle mass associated with each of 
these frequency peaks calculated with non-linear LK fits applied to the temperature 
dependence of their amplitude, as described in section 12.2.41 For low field side data, 
figure 13.61 presents these data and LK fits for the five measured frequencies, where 
data from 90 mK and above were used, for the reason described in section 13.1.31 All 
masses calculated are more precise and accurate than in the previous work by Borzi, 
as error bars are smaller and only two parameter fits were used as opposed to threqf . 

The complexity of the high field side does not simply arise due to the additional 
intrinsic inner magnetic field produced by the sample itself, producing a complex 
dependence of the oscillations on B. Using the work by Perry et al. [27], we calculated 
its magnitude, and obtained around 0.015 T, which, when added to the measured 
field values, does not make any visible difference to the spectra; much higher offsets 
are required to produce broadening. Moreover, the multiple peaks of the high field 
side do not originate from simple spin splitting as one would normally expect in a 
metamagnetic syste In such a scenario, each peak of the low field side splits into 
only two as the system crosses the metamagnetic transition. Part of the complexity of 
the high field side spectrum may originate from spin splitting, but many more peaks 
are found. Finally, the quasiparticle masses associated with these peaks are not all 
consistent within each group; instead, we find different values ranging between 7 and 



7 Only three are visible in figure 13.51 with a larger field range, the central one splits into two. 
Larger field ranges reduce the clarity of the overall spectrum. 

8 We will show in section [4.4. II how three parameter fits induce systematic errors to LK fits. 
9 see section ll.3.5l for details 
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Figure 3.8: Dingle analysis of the amplitude of the 1.8 kT frequency in the low field 
side with the field in the c-axis direction. The experimental amplitude function is 
shown with a solid blue line, and the broken red curve corresponds to eq. 13.21 



17 m e . 

Since the Fermi surface undergoes an evolution across the metamagnetic trans- 
ition, it is interesting to look at the field dependence of the spectrum shown in figure 
13.71 This calculation was performed by taking Fourier transforms over 150 small over- 
lapping inverse field intervals of width AX =0.02 T -1 , and plotting the intensity of its 
power spectrum at each average field position. Resulting horizontal lines of intensity 
correspond to quasiparticle orbits, and the vertical feature near 8 T is produced from 
the Fourier transform of the metamagnetic peaks. 

One can see, in the low field side, that the dHvA frequencies at 0.45 1.8 and 
4.2 kT undergo a dramatic increase as they come towards 8 T. The peaks at 1.8 and 
4.2 kT decrease again in the high field side towards 0.16 and 4.0, respectively, while 
the peaks at 0.45 and 0.9 kT disappear. Finally, one observes that the complexity of 
the high field side gradually increases with field. 

3.2.3 Dingle analysis of c-axis dHvA 

The low field side data at c-axis, 25 mK, may be used in order to calculate the mean 
free path of the quasiparticles, using the method described in section 12.2.71 Only the 
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1.8 kT peak possesses a modulation that is simple enough to model, and thus will be 
the only one used. We employed the envelope extraction method, section I2.2.6[ and 
adjusted a function of the type 

CXTm* 

Ae~ x ^ . r~— . | cos(u;X + <j>)Bj{2 * pi * H AC FX 2 ) \ , (3.2) 
smh(GA _/ m*J ' 1 

with free variables A and A, assuming all the other parameters known. 



Figure 13781 presents the analysis, with the envelope of the 1.8 kT peak plotted in 
blue and eq. l3.2l overlaid in red. For most of the low field range, reasonable agreement 
was found. The data differs from the model near the metamagnetic transition, as a 
result of dephasing due to spin splitting, which will be discussed in section 14.21 We 



extracted the value A = 36 T, which, using eq. I2.9[ corresponds to a mean free path 
of around £ = 270 nm. 



3.3 Angular dependence of dHvA 

The angular dependence of dHvA in S^R^Oy was measured in Cambridge only 
during this project. Traditionally, the angular study of dHvA is used to obtain 
information about the three dimensional shape of a Fermi surface. In two dimensional 
layered materials, as we described extensively in section ll.3.4[ the angular study of 
dHvA allows one to determine in detail the corrugation of the various Fermi surface 
pockets. It can also help one to determine their symmetry and position inside the 
BZ. This type of analysis was not performed completely in the past in Sr3Ru207, and 
we present here the data resulting from the first attempt at describing the complete 
three dimensional shape of the Fermi surface in S^R^Cy 

The samples were cooled down and kept at a constant temperature of 30±5 mK. 



These were oriented along crystallographic directions, described in section 12.51 The 
rotation of the samples was performed using the system described in section 12.1.31 
and appendix |Ej such that sample C698I was rotated from [001] towards [110] and 
C698A towards [100]. We found that errors in the installation of one coil and in the 
cutting of sample C698A resulted in a difference of alignment between the samples, 
which lies in the plane of rotation o Over 35 angles at steps of 1.7°, second harmonic 
AC susceptibility was recorded along the full field range between 18 and 2 T, using 
the method described in section 12.1.31 The resulting data possessed an unexpected 
complexity which, as we will present in section [3.3.31 is n °t explained by the approach 
of Bergemann [35] , and a model will be presented in section 14.21 suggesting an origin 
to this phenomenon. Appendix |E] presents plots of the raw data at all angles, for 



10 o , 



see appendix lEl for details 
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the interest of the keen reader, where the complexity of the beat patterns can be 
appreciated, as well as the behaviour of the metamagnetic transition signal. 

Sr 3 Ru207 is close to tetragonal; although the a and b axes are not equivalent from 
a symmetry point of view, these lattice parameters are equal in length and make 
those directions indistinguishable from a Laue experiment. Effectively, Sr 3 Ru207 is 
orthorhombic from the crystal structure, which lacks 90° rotational symmetry, and 
hence the Fermi surface lacks that symmetry as well. Using the work of Mazin and 
Singh [37] as a guiding line, we nevertheless assumed that the a and b were almost 
equivalent, in that the Fermi surface is almost symmetrical under rotation by 90°. 
We found later that, from the work of A. Tamai et al. (described in section H~Tj) . 
this is actually the case. 

We discuss in the first section an aspect of the data independent of dHvA but of 
importance in S^R^C^, that of the phase diagram and the mapping of the metamag- 
netic transition in the field-angle plane. We will then present the angular dependence 
of the dHvA spectra, followed by a discussion of its meaning. Finally, we will describe 
beat patterns, along with an attempted simulation. 



3.3.1 Mapping of the metamagnetic transition 

The field derivative of the AC susceptibility provides information about the metamag- 
netic transition and the phase diagram of S^R^Oy. The field values at which the 
metamagnetic features appear depend strongly on field angle [38J , and it was possible 



to map out these in the present data set. Figure 13.91 shows the positions in field as 
a function of field angle for two directions of rotation, from [001] towards [110] and 
[100]. It is interesting to note that the two directions are not equivalent in terms of 
metamagnetic behaviour. At zero angle, one can see three features, the nematic phase 
lying between the top two. Increasing the angle, we observe that the top transition 
line joins the middle one at an angle of about 15°, for the [110] direction, and at 
about 36° for the [100] direction. Moreover, the middle transition moves towards 6 T 
at high angle in a different manner for both directions. Also, there is a feature arising 
at higher fields near 20° which moves fast to lower fields and splits into two, and also 
behaves differently in the two directions. 

These remarks about the direction of rotation were first observed with this data. 
It is not absolutely surprising that the metamagnetic behaviour should be different 
in these directions, as they are not equivalent! n l. In other words, the magnetism of 
Sr 3 Ru 2 07 is not isotropic in the plane. It has not been observed before, and suggests 
that a new phase diagram needs to be measured with the highest purity samples. It 



1 see section H.2.2l for details about the crystal symmetry 
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Figure 3.9: Phase diagram of Sr 3 Ru207. Mapping of the positions of the features 
in the magnetic susceptibility, taken from the second harmonic data of the rotation 
study, for both samples. Sample C698I was rotating towards [110] and sample C698A 
towards [100], where the zero of the angle is at [001]. Three features seen in the c-axis 
data are visible at zero angle, and join to a single one at higher angle. The nematic 
phase region lies between the features at 7.9 and 8.1 T near zero angle. 



was also recently found, with specific heat experiment by A. Rost o, that there are 
significant differences in the thermal properties of Sr3Ru207 with rotation direction, 
such that they are not equivalent when the field is in the plane, whether it points 
towards [100] or [110]. 



3.3.2 Angular dependence of dHvA spectra 

Figure 13.101 presents the angular dependence of the dHvA spectra for both samples 
in the low field side, one per sample, or direction of rotation. The colour in the plot 
corresponds to the square root of the amplitude of Fourier transform power spectra 
for each field angle centred at the angle value. The Fourier transforms were taken 
from 5 to 6.5 T, except at high angles, where the interval was taken from 5 up to the 
position of the metamagnetic transition^. The intensity of the data was furthermore 



12 Unpublished work which will form part of A. Rost's thesis [62] . 

13 As can be seen in figure 13.01 the metamagnetic transition enters the region from 5 to 6.5 T 
near 30°. The Fourier transform of the metamagnetic signal produces large intensities near zero 
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Figure 3.10: Angular dependence of the dHvA spectra in the low field side for sample 
C698I, top, which was rotated from [001] towards [110], and sample C698A, botttom, 
which was rotated from [001] towards [100]. 
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Figure 3.11: Angular dependence of the dHvA spectra in the high field side for each 
sample, as described in figure 13.101 
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cut off in order to provide contrast for low amplitude features. 

We observed that the frequencies possess the normal upwards curvature expected 
from two dimensional materials (described in section fl.3.4|) . Figure l3.12[ top graph, 
shows a comparison between this data and the function 



cos(0)' 

with Fo, the c-axis direction frequency values, from table 13. 1[ of 0.15, 0.43, 0.91, 1.78 
and 4.13 kT I 1 I. We conclude that Sr 3 Ru20y is quasi two-dimensional in the low field 
side, and that the slight deviations from the function F / cos(6 l ) originate from k z 
dispersion. Finally, we note that a similar result was obtained in both directions of 
rotation. 

Turning our attention to the high field side, shown in figure 13.11} we observed 
that the behaviour is completely different from the low field side, and that it does not 
simply follow the function Fq/ cos(6 i ). We mentioned in section 13. 2.21 that at c-axis, 
the peaks were multiple, a feature that we found is present at all angles. Secondly, we 
observed that the amplitude also possesses a complex angular dependence, such that 
some frequencies appear while some others disappear at specific angles. In particular, 
the 1.6 kT frequency vanishes between 5° and 10°, while a 0.9 kT peak appears near 
the same angle, and a peak near 0.5 kT appears between 15° and 20°. The large 
structure near 4 kT is also interesting. It does not have the same behaviour as that 
of the 1.6 kT peak, but still vanishes at a lower angle than other peaks and is broken 
down into more than 25 peaks. It follows an upwards curve with increasing 9, though, 
and features very complex beat patterns. Finally, the group of peaks near 7.9 kT is 
very faint and it is difficult to observe its angular dependence. Note that the spots in 
the upper right corner correspond to low frequency time-periodic noise, and contains 
no dHvA[£. 

Figure I3.12[ bottom panel, presents the high field side data overlaid with curves 
of F / cos(9), using for F the values 0.40, 0.85, 2.95, 3.7 and 4.5 kT. We observed 
that the the data did not follow this trend, but curved either significantly more or 
less than F / cos(#). These observations suggest that the system in the high field 
side has more three dimensional character than in the low field side. And finally, no 
significant differences were seen between the two directions of rotation. 

In conclusion, we stress that the electronic structure of Sr 3 Ru207 changes possibly 



frequency, which may interfere with interesting features. 
14 Note that these curves are not fits. 

15 It is due to the voltage limited mode (see section !^. 1.3p . which makes the time and 1/B frequency 
domains similar. When the field is swept exponentially as a function of time, noise at low time 
frequencies appear periodic in — log(-B/i?o), which is closer to 1/B than in a linear field sweep. 
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Figure 3.12: Angular dependence of dHvA for sample C698I (rotation towards [110]), 
in the low field side (top) and the high field side (bottom), with expected frequency 
curves overlaid on top, in white dashed lines. The dashed lines correspond to the 
function F = F / cos(0). 
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more radically than was previously thought across the metamagnetic transition. The 
correspondence between the low and high field sides is not as obvious as previously 
thought. Moreover, the phenomenon of dHvA peaks appearing and disappearing 
with angle is difficult to explain, as is the change in the angular dependence of the 
frequency values across the transition. 



3.3.3 Interference patterns in the angular dependence 

dHvA data contains more information than that which is immediately visible in its 
Fourier transform. As was described in section ll.3.4[ closely spaced frequencies nor- 
mally appear as weakly split peaks in the transform, but may have more dramatic 
consequences in the overall appearance of the raw data, as beat patterns affect their 
amplitude. Hence, in that view, it is sometimes the case that examining the amplitude 
of modulation of frequencies is more effective than analysing its Fourier transform, 
as was introduced by Bergemann and co-workers [53]. In simple 2D materials, the 
beat patterns are unique fingerprints of the specific type of corrugation of the Fermi 
cylinders forming the FS, but the interference may originate from other phenomena 
as non- linear magnetisation and spin splitting (introduced in section ll.3.5j) . 

In Sr 3 Ru207, it is well known that non-linear magnetisation occurs near 8 T, 
accompanied by changes in the dHvA frequencies in the low field side of the meta- 
magnetic transition. Consequently, we anticipated that this analysis would become 
more complex than for Sr 2 Ru04. Effectively, performing the extraction of the beat 
patterns of all the frequencies present in the low field side revealed patterns that were 
too complex for us to model except in the case of the 1.8 kT frequency. This Fermi 
cylinder has consequently become the centre of our attention in this section. 

Figure 13.1 31 shows the amplitude of the frequency at 1.8 kT in the low field side as 
a function of field and angle, for both rotation directions. This was calculated using 
the method described in section |2.2.6[ using a rounded top-hat window of width 0.3 
kT, centred at (1.8 kT]/cos#, which was chosen to include all the data even where 
the frequency changes 16 !. The regions in the top right corner, top plot, and both 
top corners, bottom plot, lie above the metamagnetic transition, and were excluded 
from the computation. These calculations revealed that the beat patterns for both 
directions of rotation are similar, and the corrugation of this part of the Fermi surface 
appears isotropic apart from the in-plane shape (of warping parameters fc^o), which is 
not detectable with dHvAl 17 !. We deduced from this fact that all warping parameters 
with fi > 0, v > are absent or very small (for example, k 2 i, /c 41 , fc 2 2> ^42, etc.), which 
leaves an isotropic warping of the type fi = 0, v > 0, most probably k i or /c 2- Our 



16 Since the frequencies evolve with magnetic field. 

17 This information was provided to us by ARPES, section H~T1 
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Figure 3.13: Amplitude of the oscillations at 1.8 kT as a function of field and angle 
in the low field side, for sample C698I, which was rotated from [001] towards [110], 
top, and C698A, which was rotated towards [100], bottom. 
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Figure 3.14: Comparison of the data and simulations of beat patterns as a function 
of field and angle, using a single parameter in the cylindrical expansion of koi = 
4.4 x 10 6 m _1 . The plots on the top presents the same data as in the previous 
figure for the rotation towards [110], left, and [100], right, overlaid with black curves 
corresponding to the nodes of the beat pattern. The plots on the bottom present the 
simulation using the single warping parameter, using the same angular span as the 
data. 



second observation was that the c-axis region appears normal, with a simple beat 
pattern that may be used to extract the value of k 0u , where effectively, we remarked 
that the beat nodes, or zeros, curve downwards with angle as expected (see, for 
example, figure fLTT] section H.3.4j) . 

However, we found that the high angle data appears anomalous, where effectively, 
we expected to find a Yamaji angle. Instead, we observed rapid beating, of which 
the period becomes shorter as field increases towards the metamagnetic transition, 
in opposition to the normal situation, where beat nodes appear at equal intervals 
of inverse field. We will show in section 14.21 that this can be partly explained by 
non-linear magnetisation and spin splitting resulting from a peak in the DOS. We 
note that in the region where the beat pattern is very rapid, the angle step of about 
1.6° is not small enough and the pattern is under-sampled, and warn the reader that 
the appearance of the data can be misleading: the beat pattern in a data set where 
the beating period in angle is comparable or smaller than the sampling rate may 
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appear in some cases as a checkerboard pattern, in others as horizontal nodes, where 
in reality the nodes may cross the vertical with a sharp angle. 

In order to find whether the warping is of type fcoi, &02, simulations were produced 
in order to compare with the data. We also used equation 11.251 for plotting the posi- 
tions of the zeros of the beat patterns. Figure 13.141 presents such an analysis, where 
data is shown on the top overlaid with beat pattern nodes in black, and simulations 
on the bottom. Although the agreement is far from perfect, the best estimate we 
can make is that k i = 4.4 x 10 6 m _1 is the dominant warping parameter, which 
reproduces the low angle part of the data. With k§ v = 4.4 x 10 6 m _1 , v > 1, we also 
reproduced properly the c-axis direction beat pattern, but the nodes fell towards zero 
field too sharply in angle compared to the data. Using a combination of more than 
one parameter of the type k Ql/ did not reproduce the data, but instead changed the 
beat pattern at c-axis. Since parameters with fi > 0, v > were already excluded, 
we conclude that no warping parameter of significant amplitude produces the addi- 
tional modulation to the pattern from simple warping ko\. We therefore associated 
the remaining modulation to non-linear magnetisation and spin splitting, as will be 
discussed in section I4T21 

The envelope extraction was also performed with the high field data, where even 
more surprising results were found, in line with the observations of figure 13.111 The 
calculation was performed using a window of width 0.4 kT, centred at (1.6 kT) / cos 9, 
over the region between 7 and 18 T, excluding the region near c-axis which lies in 
the low field side of the transition. Figure 13.151 shows the result, the shape of which 
reminds the author of a simple exclamation mark. We observed that the oscillations 
around 1.6 kT are completely suppressed from an angle of about 5° from c-axis, a res- 
ult that is very difficult to explain. We stress that only a very complex combination of 
many warping parameters could reproduce such an angular dependence of the oscil- 
lations, which is not probable. Moreover, even though non-linear magnetisation and 
spin splitting are likely to produce complex beat patterns in the high field side, they 
cannot suppress the oscillations altogether. Consequently, we found no explanation 
for this puzzling behaviour. 

We conclude this section by stating that the amplitude modulation of the 1.8 kT 
frequency does not behave like the theory predicts at all fields and angles except 
near c-axis in the low field side. In the low field side, we found a crossover from 
normal to abnormal behaviour at an angle of about 15°, above which we observed 
fast beat patterns that are not periodic in inverse field and assume that this unusual 
modulation is the result of non-linear magnetisation and spin splitting, which will 
be modelled in section I4T21 In contrast, we discovered that the high field side of the 
transition exhibits a behaviour that is not explained either by warping or non-linear 
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Figure 3.15: Amplitude of the frequency at 1.6 kT in the high field side, as a function 
of field and angle, for sample C698I, which was rotated from [001] towards [110], top, 
and C698A, which was rotated towards [100], bottom. 
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magnetisation and spin splitting and remains mysterious. 



3.4 Field dependent quasiparticle masses 

This section presents a reinversigation of the enhancement of the quasiparticle masses 
near the quantum critical end point in S^Rx^CV Since the beginning of this project, 
discrepancies have been found between the measurements of Borzi et al. [40J (see 
section ll.3.6[ figure II . 13 j) and ours. Several problems with eddy current heating and 
high noise levels prevented us producing a definitive answer to the problem until we 
performed measurements with the system in Cambridge. These prove that the conclu- 
sion in the work of Borzi that there exists an enhancement of the quasiparticle mass 
of the 1.8-1.6 kT and 4.2 - 4.0 kT peaks was erroneous. It furthermore demonstrates 
that none of the masses of the other measured peaks exhibit any increase either. 



We saw in section 13.11 that we have performed extremely careful checks on the 
thermometry during the Cambridge experiments, and showed that no eddy currents 
were present. In such a situation, we have eliminated the main source of systematic 
errors on the thermometry. Moreover, in the Cambridge data the signal to noise was 
high enough for us to be able to remove the third parameter in the LK fit *f[ which 
made our analysis much less prone to systematic errors, and thus more reliable than 
in previous work. We came to the conclusion that a combination of field varying 
low signal to noise ratio and three parameter LK fits may have led to the spurious 
continuous enhancement of the quasiparticle mass in the work of Borzi. This section 
presents the quasiparticle masses calculated with two parameter fits, for all the major 
dHvA frequencies of the low field side. We will explain in the next chapter, section 
14.4. 1 ^ how Borzi et al. probably obtained enhancements of those masses through 
systematic errors. 



3.4.1 All masses, at c-axis and 10° 

We have measured during the Cambridge series of measurements dHvA as a function 
of temperature using sample C698I with its c-axis aligned with the magnetic field, 
and simultaneously sample C698A aligned with its c-axis 10° away from the field 
towards [100], at 10 different temperatures between 90 and 500 mK. Subsequently, 
we performed a similar measurement using C698A aligned with the field and C698I 
at 10° away towards [110], at 8 different temperatures within the same boundaries. 
Measurements performed at lower temperatures were not used for this section due to 
possible thermal equilibrium issues, as explained in section in section 13.1.31 

18 see section 12.2.41 for details LK fits and the number of parameters 



Chapter^ The de Haas van Alphen experiment on Sr^R^Oj 



114 



c-axis 




|<D<D<Dl><IIO<D <D 9 ^ 
go 



□ 0.15 kT 

O 1.8- 1.6 kT 

A 4.2 - 4.0 kT 



8 10 12 

Average Field ( T ) 



14 



16 



c-axis 



<> 



D 0.43 kT 
O 0.9 kT 



4.5 



5 5.5 6 

Average Field ( T ) 



6.5 



- 1 T— 




□ 0.15 kT 
O 0.9 kT 
A 4.2 - 4.0 kT 



8 10 12 

Average Field ( T ) 



16 



7.5 



18 
16 
14 
12 

£ 10 

(A 
tfl 
10 

S 8 

□_ 

O 

6 
4 
2 




10 



□ 0.43 kT 
O 1.8 kT 



5 5.5 6 

Average Field ( T ) 



Figure 3.16: Field dependence of the quasiparticle masses in sample C698I, measured 
with the magnetic field aligned in the c-axis direction, left, and at 10° towards [110], 
right. The frequencies used were of 0.15, 0.43, 0.9, 1.8 and 4.2 kT. The inverse field 
window of 0.01 T _1 is approximately 5.5 data points wide. 



We chose to use the frequencies present in the low field side for the calculation 
of the field dependence of the quasiparticle masses, those at 0.15, 0.43, 0.9, 1.8 and 
4.2 kT (see table I3.ip . The reason for this choice lies with the simplicity of the 
low field side and the overwhelming complexity of the high field side. In the c-axis 
direction, the low field side frequencies at 0.15, 1.8 and 4.2 kT were in the high field 
side as well, at values of 0.14, 1.6 and 4.0 kT, while those at 0.43 and 0.9 were present 
only in the low field side (see figure l3T|) . In contrast, at 10°, the frequency at 1.8 kT 
was not present in the high field side, but the one at 0.9 kT was (see figures 13.111 and 
I3.10p . The ones among these that extended to the high field side in either directions 
were used over the whole field range, but the others were used only in the low field 
side. 

Figures 13. 16} for C698I, and 13. 17} for C698A present the field dependence of the 
mass of these frequencies, at c-axis in the left side graphs, and at 10° on the right. 
The calculations were performed using an inverse field window of 0.01 T" 1 , and a 
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Figure 3.17: Field dependence of the quasiparticle masses in sample C698A, measured 
with the magnetic field aligned in the c-axis direction, left, and at 10° towards [100], 
right. 



non-linear LK fit of the type (see section 12.2.4(1 



LK(X ,T) = P 1 



P 2 U.7X T 
sinh(P 2 14.7X T)' 



where P\ and P 2 are the only two fit parameters. One hundred such fits were per- 
formed for each frequency, between 18 and 4 T, and the quality of each single one was 
evaluated by eye, along with the temperature dependence of the Fourier spectrum at 
each field location, in order to verify the nature of the data that was being used, and 
all fits based on inappropriate data sets were rejected. The error bars represent the 
fit error, as defined by the least squares method. For most of the fits, the third para- 
meter, if it had been introduced, would have been small compared to the amplitude 
Pi of the LK function. Nevertheless, three parameter fits were also performed using 
only the 1.8 - 1.6 kT and 4.2 - 4.0 kT frequencies, and no difference was observed in 
the resulting field dependent masses, and are not shown here. 

We observed no significant enhancement of the quasiparticle mass in any of the 
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low field side frequencies as a function of field, either at c-axis or at 10°. The only 
variations that were seen always corresponded to regions where the signal to noise 
vanishes, for example at nodes of beat patterns^, where not only enhancements but 
also decreases in mass were observed, and all of those were of no more than 4 electron 
masses. Moreover, in contrast to the first harmonic data of the previous work by Borzi 
et al., in these data sets the signal to noise does not vanish near the metamagnetic 
transition, and the metamagnetic background only becomes significant at field values 
much closer to the transition (see [3D], figure 1, compared to !3.4|) . In that region, only 
windows covering the metamagnetic background signal were discarded^. Thus, we 
conclude that there is no enhancement of the quasiparticle mass of any of the Fermi 
surface sheets observed with dHvA near the quantum critical end point in Sr 3 Ru 2 07. 



3.5 dHvA in the nematic phase 

The discovery of a new phase of the electron liquid in SrsRu207 in 2005 [TT] has raised 
many new ideas, but also many questions about its nature [4"H 135] . The findings 
about its nematic nature in 2006 [32] attracted significant interest in the theoretical 
community [631 E3J . One of the issues that have been raised regards the behaviour of 
the FS in this nematic phase, and the most appropriate probe for gaining this type of 
information is dHvA, the only other experiment that can measure FS sheet specific 
information being ARPES, which cannot be performed in a magnetic field. This 
section presents the first successful observation of dHvA inside the nematic phase, 
which was carried out in Cambridge. 

The oscillations were first discovered in 2007 in St Andrews but exhibited a very 
faint amplitude. The experiment was repeated in Cambridge using second harmonic 
and the voltage limited mode described in section 12.1.31 and very low noise data were 
obtained in and around the nematic region. The oscillations were studied on axis as 
well as at 10° as a function of temperature. We also discovered that they possessed 
a strong angular dependence, which motivated a detailed analysis as a function of 
angle. The first part to this section presents the raw oscillations as a function of 
temperature, for both samples, and their detailed analysis. Then, the second part 
shows how they evolve with field angle. 
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Figure 3.18: Top left Raw second harmonic data for sample C698I, with the c-axis 
aligned with the magnetic field, at 15 different temperatures from 25 to 500 mK. 
Top right The background substraction scheme is shown in red on top of the raw 
data, on the left, and the signal after background substraction is presented on the 
right. Bottom left Fourier transforms of the data in the top right plot, at the various 
measured temperatures. Bottom right Square root of the integral of the peaks in 
the power spectra as a function of temperature, along with two parameter LK fits, in 
solid lines. Fit results are given in the legend. 



3.5.1 Temperature dependence of the oscillations 

Figure 13.181 presents the temperature dependence of the dHvA oscillations in the 
nematic phase for sample C698I, with its c-axis well aligned with the magnetic field. ^] 
These were obtained at 15 temperatures between 25 and 500 mK. The data are of 
high quality and the noise was hardly visible at all in comparison with the oscillations 
outside the phase. 



19 See, for instance, the 1.8 kT frequency near 5 T 

20 The Fourier transform of the metamagnetic background leads, as in figure 13.71 to a large peak 
at zero frequency, the tail of which extends to the regions of interest, and could lead to systematic 
errors. 

21 The error is of about ±1° in the [110] direction, the plane of rotation, and ±3° in the [100] 
direction. 
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Figure 3.19: Same analysis as in figure I3TT91 for the data taken with sample C698A. 



The top left plot of figure 13.181 presents the raw data for all temperatures. One 
can see the metamagnetic transition signal, first discussed in section 13.1.11 We men- 
tioned that the first feature at 7.5 T on the left is a metamagnetic crossover, and 
the second and third are the first order transitions delimiting the nematic phase. 
We observed between these a few cycles of oscillations, which repeated with identical 
phase at different temperatures, and their amplitude decreased when the temperature 
increased. The top right plot presents how the background substraction was done, 
in red, where a fifth order polynomial was adjusted to the data between appropri- 
ate field values such that it properly reproduced the magnetic background, and was 
subsequently substracted from the data. Any lower polynomial did not yield the 
appropriate shape, while higher orders reproduced partly the oscillatory component, 
which in turn resulted in removing part of the signal. The result after removing the 
background in shown on the right, using the same temperature ordering as in the 
figure on the left. 

From this data set, we carried out Fourier transforms, shown in the lower left plot. 
Two frequencies were found, at 1.0 and 2.3 kT, at all temperatures. We integrated the 
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peaks and plotted the result as a function of temperature, shown in the lower right 
graph. The distributions followed the LK function, and two parameter LK fits were 
performed, using data at temperatures above 90 mK 3 where we obtained masses 
of 6.2±0.2 and 6.8±0.1 electron masses for the 1.0 and 2.3 kT peaks respectively. 
Figure 13.191 shows the same measurements and analysis performed on sample C698A. 
Very similar results were obtained, with frequencies of 1.0 and 2.5 kT, and masses of 
5.7±0.2 and 6.8±0.6 electron masses, respectively. 

We may speculate about an origin for the measured frequencies and masses. We 



know already from the field dependence of the dHvA spectra, figure 13 .7\ section 
13. 2. 2\ that the dHvA peak at 1.8 kT increases in frequency to about 2.5 kT near 
the metamagnetic transition, but the group of peaks near 0.9 kT does not evolve 
significantly. It is therefore appropriate to associate the peak at 2.3 or 2.5 kT to 
the one at 1.8 kT in the low field side, and the one at 1.0 kT to the one at 0.9 kT. 
Although the masses do not agree extremely well, since the 0.9 kT group has a mass 
of about 8 m e and the 1.8 kT peak a mass of about 7 m e , we argue that systematic 
errors due to the difficult background substraction might have led to underestimating 
the mass within the phase. 



3.5.2 Angular dependence of the oscillations 

dHvA oscillations were measured at a base temperature of 25 mK, at 16 different 
angles between 0° and 25° for sample C698I, and between -14° and 16° for sample 
C698A, with an angle step of about 1.6°. The difference in angle spans originated from 
a misalignment of the samples with respect to the field (see section [3731 for details). 
Figure 13". 201 presents the angular dependence of dHvA obtained on both samples. The 
top left plot presents data from C698I, where the magnetic field was rotated from 
c-axis towards [110], while the top right plot relates to C698A, where the field was 
rotated towards [100]. We observed that the oscillations were only present near 0°, and 
were significantly suppressed as we rotated away the magnetic field, where effectively, 
dHvA vanished at an angle of about 5-6°. In order to emphasise this observation, 
we performed Fourier transforms using the same background substraction method 
as described previously, shown in the bottom plots. Both peaks at 1.0 and 2.3 kT 
decreased in intensity as we rotated away from c-axis. 

Although the temperature dependence and the frequency values for the dHvA 
peaks that we found in the nematic phase lead us to associate these to the low field 
side dHvA spectra of Sr 3 Ru 2 07, this strong angular dependence could be related to 
the high field dHvA, where we observed that frequencies disappear and appear as a 

22 This was done for the same considerations as in the previous section regarding low temperature 
data. 
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Figure 3.20: Top left Raw second harmonic data taken on sample C698I, as a function 
of the angle between its c-axis and the magnetic field, at a temperature of 25 mK. 
The field is rotated in the crystallographic direction of [110]. Top right Same type of 
data taken on C698A, and the rotation is performed towards [100]. Bottom Fourier 
transforms of the data in the nematic phase as a function of angle for sample C698I 
(left) and C698A (Bottom). 

function of angle. If the peak at 2.3 or 2.5 kT really corresponds to the 1.8 - 1.6 kT 
frequency, then it is consistent to find that it vanishes as a function of angle inside 
the nematic phase. It is not the same for the peak at 1.0 kT, since we know that in 
the high field side, it is not present around the c-axis, but appears at around 5° (see 
figure I3TTTI) . It is then, at this point, difficult understand what we observed inside 
the nematic phase. 



Chapter 4 



Discussion 



We present in this chapter an analysis and an interpretation of the data presented 
in the preceding chapter. We provide a complete model for the Fermi surface of 
Sr3Ru207, as well as a an interpretation of its properties near the QCEP at the meta- 
magnetic transition. With our data alone, we could not produce a unique model 
for the FS, but this has become possible by the use of a combination of our data, 
ARPES measurements that were performed by A. Tamai et al. and magnetocaloric 
oscillations experiments carried out by A. Rost (36J [62] . The model introduced by A. 
Tamai et al. suggests that metamagnetism may be produced by a peak in the DOS 
situated near a new small FS pocket which was subsequently discovered by A. Rost 
in magnetocaloric oscillations. This suggestion is consistent with our data showing 
that no quasiparticle mass enhancement is present for any of the previously known 
FS pockets. We unfortunately provide no definite proof but only partial evidence in- 
dicating that the new FS pocket is responsible for metamagnetism and the previously 
observed quantum critical properties. 

We moreover make use of various models in order to interpret some unusual prop- 
erties of our dHvA data, the anomalous beat patterns in the low field side and the 
highly split spectra in the high field side. For the first problem, we introduce a simple 
generic model of metamagnetism that features a peak in the DOS, which we use to 
simulate the field and angle dependence of the amplitude of the 1.8 kT frequency. We 
obtain only partial agreement with the data, which indicates that the physics taking 
place are more complex. For the second problem, we use the theory for magnetic 
breakdown in order to show that a large number of almost degenerate orbits can arise 
from only a few symmetrical breakdown points in the BZ. We provide there the most 
probable hot spots in the BZ where magnetic breakdown could occur, and calculate 
the size of the new expected cyclotron orbits. The two models presented are far from 
exact. They were only developed to point to probable sources for these anomalous 
properties; detailed calculations are beyond the scope of this thesis. 
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Finally, we discuss the absence of an enhancement of the quasiparticle masses 
near the QCEP and the appearance of dHvA inside the nematic phase. For the 
first of these phenomena, we provide a mechanism which was probably responsible 
for producing spurious mass enhancements through systematic errors in the non- 
linear LK fits of the previous work of Borzi et al. [40J. We give examples of an 
artificial mass enhancement, but we provide detailed evidence in appendix IF1 using 
numerical simulations. We furthermore discuss the implication of the absence of a 
mass enhancement by comparing with electronic specific heat data and measurements 
of the A coefficient of the resistivity. Regarding the nematic phase oscillations, we 
provide an interpretation for their presence, and argue that their suppression with 
angle may be related to putative nematic domains. 

4.1 Zero field Fermi surface 

In a dHvA experiment on a two-dimensional material, one cannot determine the in- 
plane shape of the FS (see section H.3.4p . In such a case, one knows the sizes of the 
various FS pockets, but not their position in k-space or their form and orientation 
(the terms with fi ^ in the expansion of eq. 11.21 j) . nor whether they consist of 
electrons or holes. In order to orient ourselves in this respect, ARPES experiments 
were highly desirable. 

Such experiments were performed by Tamai and co-workers at the Stanford Syn- 
chrotron Radiation Laboratory, using single crystals of Sr 3 Ru20y provided by the 
author, originating from the ensemble studied in section [2T41 [36J. High quality data 
were obtained, and a complete interpretation of the low field side FS was constructed, 
using both ARPES and dHvA data from this work. These results are part of this 
thesis work only inasmuch as they used low field dHvA frequencies and quasiparticle 
masses from this project, which had first been measured by Borzi et al. [10] . and 
subsequently remeasured with more precision by the author. They are presented here 
for the reason that they are essential to the interpretation of the main results of this 
project. We present in this section the photoemission data and its interpretation, 
along with a comparison to dHvA data. 

Moreover, quantum oscillations in the magnetocaloric effect were discovered in 
Sr 3 Ru207 by A. Rost|62j. featuring an additional low frequency orbit that was not 
detected by the author of this work using dHvA, due to lack of sensitivity in this 
frequency region^. This information has been important for the interpretation of 
our results and is presented as well. Therefore, we first describe the measurements 

1 Second harmonic dHvA possesses a sensitivity that vanishes at zero frequency as F 2 , while the 
magnetocaloric oscillations feature good sensitivity at low, but not at high frequencies. 
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• FFT Data 
LK: m=9.8 ±0.6 




Figure 4.1: Quantum oscillations in the magnetocaloric effect. Left Fourier transform 
of the oscillations between 4 and 7 T. Right Temperature dependence of the amplitude 
of the peak at 110 T, circles, along with a fit of dLK(T)/dT, solid line. Taken from 
A. Rost [52]. 

performed with the magnetocaloric effect, followed by those using ARPES. We then 
introduce the model constructed by A. Tamai et al. and an interpretation for the zero 
field FS of Sr 3 Ru 2 7 . 

4.1.1 Oscillations in the magnetocaloric effect 

Quantum oscillations in the magnetocaloric effect in S^R^C^ were found by A. 
Rost[62j, using samples that were grown by R. S. Perry and characterised by the 
author of this project, as described in section 12.41 Two high purity samples were 
used, labelled C698J and C698K. Oscillations were found both above and below the 
metamagnetic transition field value. In the low field side, the transition, all of the 
frequencies measured in this project were reproduced except that at 4.2 kT, and 
the extracted quasiparticle masses were in agreement as well. An additional very 
low frequency peak was discovered, which in second harmonic dHvA possessed an 
amplitude that was too low for detection, but became crucial for the interpretation 
of physics of Sr 3 Ru207. 

Figure 14.11 left panel, shows the Fourier transform of the magnetocaloric oscil- 
latory signal between 4 and 7 T, where one can see peaks at 0.11, 0.15, 0.43 and 
1.8 kT, and a faint group near 0.9 kT. Excluding the peak at 0.11 kT, the other 
ones correspond to the same FS sheets as those observed in this work (see table 13.11 
p. !9"oT) except for the peak at 4.2 kT, which was not observed in the magnetocaloric 
effect. As presented in the work of Rost [62], these oscillations are damped with in- 
creasing frequency due to phase smearing related to the heat relaxation time constant 
of the system. Figure 14.11 presents the quasiparticle mass analysis of the new peak 
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at 0.11 kT. As the magnetocaloric effect involves the temperature derivative of the 
magnetisation, the temperature dependence of the oscillations follows the temperat- 
ure derivative of the LK function rather than the LK function itself. By adjusting 
dLK(T) /dT, Rost obtained a quasiparticle mass of 9.8±0.6m e . This discovery will 
be of relevance in the next section in order to obtain a complete set of quantum 
oscillations peaks in Sr 3 Ru207. 

4.1.2 Angle Resolved Photo emission Spectroscopy 

ARPES data were obtained over several series of experiments using samples that were 
characterised by the author, using the full analysis method described in section 12.41 
The crystals used were the ones labelled C6103A, C6103ZA, C698C, C698DA and 
C698DB, and their properties can be found in tables 12.11 and 12.21 These possessed 
very low disorder and had impurity phase volume fractions of less than 1%, except 
for C698DA and C698DB, which contained 4% of Sr 2 RuO,4. For more details on the 
technical aspects of the ARPES experiments, the reader is referred to the review on 
the subject by Damascelli [6] and the article by A. Tamai et al. [36J. 

Figure I4T21 presents the measured ARPES FS. The plot labelled (a) presents the 
FS photoemission data in the undistorted zone, with the v2 x y/2 reconstructed zone 
indicated as a white diamond. From this image, FS contours have been extracted, 
shown in (b). Six FS pockets have been identified and labelled a±, «2, 7i, 72, P and 5. 
The first two were so labelled for their common origin with the a pocket of Sr 2 Ru04, 
which possesses d xz , d yz orbital character (see section [1.2.31 and figure fT3|) . The next 
two were labelled 71, 72 for their relation with the 7 pocket of Sr 2 Ru04, or the d xy 
orbital, although they are also hybridised partly with a band parent to the f3 pocket 
of Sr 2 Ru04, and it is not presently known with certainty whether 72 actually crosses 
ep or exists below the Fermi energy. The /3 pocket is, on the other hand, related to 
that with the same name in Sr 2 Ru04. Finally, the S pocket originates from a new 
band pushed down from the e g manifold above the Fermi level, of d x 2_ y 2 character, 
as inferred from LDA calculations of Mazin and Singh, [37] . The third part of figure 
14.21 (c), presents LDA calculations performed by Singh, for an octahedral rotation of 
6.8°. The blue lines represent the FS at the basal plane, k z = 0, and the black lines 
are for k z c = tt/2, showing the extent of the k z dispersion, or corrugation of the FS 
sheets. 

The bottom table in figure 14.21 shows a comparison between dHvA and ARPES 
data. The first row presents values of the dHvA frequencies F associated with the 
Fermi surface contours identified by ARPES, with the same notation as in (6), followed 
by the fraction of the BZ area each extend to, and their quasiparticle masses. The 
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Figure 4.2: (a) ARPES FS data measured by Tamai and co-workers |36j. The image 
was symmetrised from the first quadrant, (b) FS contours extracted from (a), (c) 
FS calculated by Singh, (d) Schematic structure of a single ruthenium oxide plane. 
Table dHvA FS areas and quasiparticle masses obtained by the author (see table l3~Tj) . 
compared with the same quantities obtained with ARPES. The error on the dHvA 
frequency and areas are of 10 T and 0.1%. 
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cross sectional areas in the c-axis direction from dHvA data and those from ARPES 
all agree very well, except for the 72 pocket, which possesses a small area and was not 
detected by dHvA. There is also a good agreement between the dHvA quasiparticle 
masses and those extracted from ARPES Fermi velocities. 

This data set provides us with a robust in-plane picture of the FS of Sr 3 Ru207. 
ARPES moreover tells us, from dispersion data, whether the pockets consist of elec- 
trons or hole^. It was discovered that the circular pocket near T (5) is electron-like, 
the square and cross shaped bands 02) centred on T are hole- like, the two surfaces 
at the M point (71, /3) are electron-like, and the very small pocket near the corner X 
of the zone (72) is hole-like. 

4.1.3 Luttinger sum rule and specific heat 

As was described in sections [1.3.11 (eq. II. Hj) and 11.3.21 (eq. II . 16j) . the total number 
of electrons detected with dHvA can be counted, as can be the total number of 
degrees of freedom, which may be compared with, respectively, the expected number 
of electrons per formula unit and the electronic specific heat at zero field. For this 
purpose one requires a model of the FS, with knowledge of whether the various sheets 
are composed of electrons or holes. 

One such model was proposed by A. Tamai et al., where these two sums were 
performed when using the ARPES areas and quasiparticle masses. It is the following: 
«i and «2 consist of holes and count as negative. 71 and f3 are composed of electrons, 
but from energy dispersion data, Tamai and co-workers have demonstrated that 71 is 
double, a phenomenon attributable to bilayer splitting of the d xy band, and so counts 
twice. The (3 band is not double, however, since as one can see in figure [1751 f . . the 
bilayer splitting was already taken into account and leads to another band situated 
elsewhere^. Moreover, and 71 are situated at the M point, and one is required to 
count them both twice, such that (3 is counted twice and 71 four times. Finally, the 
putative pocket 72 is very small and does not make a very large contribution to the 
total area, but possesses a mass that is not negligible, and appears four times inside 
the BZ. It is not presently known but rather expected from DFT calculations that 
this pocket should also be bilayer split, increasing its count to eight times. Finally, 
the 5 pocket was predicted to be bilayer split too, and hence will be counted twice. 

Since ARPES does not provide a definite answer as to whether 72 is part of the 
FS or not, we tried to find out through the Luttinger sum and the total specific heat. 
One should account for 16 electrons in Sr 3 Ru20y, since four electrons originate from 
four Ru atoms in the distorted lattice. We mentioned in section 11.2.31 how the FS 



2 A. Tamai et al., private communication. 

3 The same argument applies to the «i, «2 bands. 
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could be constructed from simple arguments, and presented a schematic plot of the 
result, figure lL5l In this model, the back-folding of the four outermost FS sheets seen 
in figure 11.51 (d) should give rise to eight electron filled BZ, one per spin flavour per 
sheet. With this FS model, ARPES and dHvA data, the Luttinger sum excluding 72 
is 

A tot = 2(1 - A ai + 1 - A a2 +4A 71 +2A P + 2A S ) + 8 
= 12 - .08 ± .09 electrons (ARPES), 
= 12 - .162 ± .006 electrons (dHvA), 

(4.1) 

where the factor of two refers to spin flavours and the additional eight corresponds to 
filled zones. In both results, there are four electrons missing, which could be hidden 
in additional filled Brillouin zones. In this view, the ARPES result is consistent with 
16 electrons, but the dHvA result misses almost a fifth of an electron. Including 72, 
assuming that it covers an area of 1% with an error of 100% of its value for ARPES, 
and using the frequency of 0.11 measured by A. Rost for dHvA, the sum becomes 

A tot = 2(1 - A ai + 1 - A a2 + AA 71 + 2A P + 2 - 8A J2 + 2A S ) + 8 
= 16- .24 ±.12 electrons (ARPES), 
= 16 - .291 ± .008 electrons (dHvA). 

(4.2) 

In this case, we obtain the appropriate number short of a quarter of an electron for 
both experiments. 

With the same two models, the total specific heat is counted by adding the qua- 
siparticle masses, where electrons and holes contribute positively, and one should 
reproduce the value of the zero field specific heat of 110 mJ/mol Ru K 2 measured by 
Ikeda et al. [21]. Excluding 72, we obtain 

\J m* = m ai + m a2 + 4m 71 + 2mp + 2ms 

= 91 ± 15m e 7 = 67±21mJ/molRuK 2 (ARPES), 
= 75.8±1.2m e ,^7 = 56±lmJ/molRuK 2 (dHvA), 

(4.3) 

and find that around half of the specific heat is missing. If we now include 72 into 
the sum, taking for the missing value in dHvA that measured by A. Rost with the 
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magnetocaloric effect, we obtain 




m ai + m a2 + 4m 7l + 2mp + 2ms + 8m. 



174 ± 29 m e 7 = 129 ± 21 mJ/molRuK 2 
159 ± 11 m e , 7 = 118 ± 8 mJ/molRuK 2 



(ARPES), 
(dHvA), 



(4.4) 



and recover the correct specific heat within experimental error. 

We conclude from the calculation of the specific heat that the band giving rise 
to 72 most probably crosses the Fermi level and that it should be included in the 
FS, whereas if we do not, half of the specific heat is not accounted for. However, 
the Luttinger sum does not seem to agree with this statement, though, as we do not 
obtain a integer number of electrons in either cases whether we include or not the 72 
pocket, but the difference is greater when we do include it. 

Finally, let us mention that A. Tamai et al. also discovered the existence of a 
saddle point very close to the Fermi level, in the region surrounding the 72 band. A 
dispersion of this type corresponds to a van Hove singularity (vHS), which produces 
a peak in the DOS if it is located at the Fermi level, and, if it is located within a 
few meV of the Fermi level, could be responsible for metamagnetism in S^R^C^ 
when it is reached by Zeeman splitting. Tamai and co-workers found evidence for 
such a peak in the DOS below the Fermi level by calculating a histogram of the 
number of states as a function of energy from a topographic map of the dispersion 
in k x ,k y space near the corner X of the BZ, around the 72 band. Figure 14.31 (a), 
shows the dispersion near the 72 pocket along the k y direction. One can observe 
that this band is very flat and of hole type. The (b) panel shows the extracted 
energy distribution curves, which we do not discuss here. Panel (c) presents the level 
curves in the vicinity of the corner of the BZ X, where the saddle point is located, 
between the band 72, in dark blue, and X. From these level curves was extracted a 
contribution to the DOS that features two peaks, shown in (d). Tamai and co-workers 
suggested that the energy scale of 1-3 meV makes these peaks good candidates for 
explaining metamagnetism in Sr3Ru207. This vHS would be reached with magnetic 
field as the 72 band expands towards the corner of the zone. They also suggested that 
this could explain the changes in dHvA frequencies observed by Borzi et al. at the 
metamagnetic transition, which could correspond to a Fermi surface reconstruction 
involving transfers of quasiparticles between bands. 



Chapter St Discussion 



129 




Figure 4.3: (a) Dispersion near the corner X of the BZ and the 72 pocket, in the k y (or 
k x ) direction. (6) Raw energy distribution curves extracted from (a), (c) Topographic 
map of the dispersion in the plane near X, where a saddle point exists, (d) Histogram 
of the dispersion in (c) as a function of energy. 

4.2 Low field Fermi surface 

The low field FS of SrsRi^Oy possesses a structure that is more complex than that 
predicted in a two dimensional system like Sr 2 Ru04, as we saw in section 13.31 In 
Sr 2 Ru04, the anisotropy of the crystal lattice leads to FS sheets that are connected 
in the k z direction, the shape of which can be expanded in cylindrical coordinates 
(section 11.3.41) . In the work of Bergemann [35], the beat patterns of the quantum 
oscillations in field and field angle were calculated very accurately in outstanding 
agreement with dHvA measurements. In S^R^Oy, we saw that the beat pattern of 
the 1.8 kT peak is not accounted for by any combination of terms in the expansion 
of that model. An additional contribution plays a role here, most probably related 
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to metamagnetism. 

We propose in this section a model where a peak in the total DOS, combined with 
charge transfers between bands, produces a non-linear variation of the k space area of 
the Zeeman split bands, resulting in interference in the dHvA oscillations. Moreover, 
we present evidence that the gyromagnetic factor of the quasiparticles g is spatially 
anisotropic and produces an additional angle dependence in the dHvA. We do not, 
however, reproduce exactly the complex beat pattern of the 1.8 kT peak, but this 
model explains relatively well some of the features of the data. A more elaborate and 
accurate theory proved too difficult to construct for the scope of this project. 

4.2.1 Anisotropic g factor 

In materials with strong crystalline field, the coupling in energy of the electron's 
magnetic moment in a magnetic field can be anisotropic. The energy term of Pauli 
paramagnetism is expressed with a tensorial form for the g factor, S ■ g • B, where 
g — ifeii + 39yy3 + kg^k, and the directions i,j,k are defined with respect to the 
crystal field symmetry. The g factor has three main axes with different values, and 
the shift in energy due to the Zeeman splitting changes with the direction of B. In 
spherical coordinates, it is expressed as 

g(9,4>) = sj g 2 xx sin 2 9 cos 2 <p + g 2 yy sin 2 9 sin 2 <p + g 2 zz cos 2 6, 

where 9 is the angle between B and the z direction and <j) that in the xy plane. The 
g factor respects the crystal field symmetry, and we take g zz as that in the c-axis 
direction, and, g xx ,g yy in the afe-plane. For the sake of simplicity, we assume that 
g xx and g yy are of very similar value, which we will find is approximately but not 
absolutely true in Sr 3 Ru 2 7 . With this approximation, the expression for g reduces 
to g{9) = v/^sin^ + ^Wtf. 

The consequence of the angular dependence of g is that properties driven by the 
magnetic field should exhibit an additional variation with angle. It is indeed the 
case with the metamagnetic transition, as we saw from figure 13.91 In that case, if we 
assume that the magnetisation follows an empirical function / of the type 

M{B)<xf(~v B g{9)B-e c ), 

where e c is the critical energy, then it follows that at the metamagnetic transition, 
the product g{9)B is constant and we have a critical field B* that follows 
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Figure 4.4: Out-of-plane magnetic field component as a function of the in-plane 
magnetic field component of the metamagnetic transitions and crossovers, taken from 
the positions of peaks in the second harmonic AC susceptibility. The data is well 

B 2 B 2 

modelled with elliptical curves, of the type — ^ + = 1, shown in red. The plots 
present data for rotations performed from [001] towards [110], left, and towards [100], 
right. The legends show the pairs of parameters a and b for each curve, which reveal 
a small anisotropy in the plane. 

This can be expressed with the in-plane and out of plane critical fields B c and B a b in 
the form of an ellipse, 

H'B 

From the work of Grigera0, that of Perry [27J E3 EH] and this thesis, we know that the 
critical field in the c-axis direction is between 7.9 and 8.1 T, and that in the a6-plane 
is between 5.1 and 5.3 T and that it evolves continuously between the two sets of 
values. Figure 14.41 shows the position of the c-axis component of the critical fields 
as a function of the afe-plane component, along with elliptical curves taken from this 
description. They agree remarkably well. Note that we do not presently know the 
exact value of g, except for the fact that it is smaller in the c-axis direction compared 
to the in-plane value by a factor of about 5.2/7.9 ~ 0.66, which corresponds to the 
ratio of the critical field values for these directions. 



4.2.2 Evolution of frequencies with magnetic field 

We presented in section I3T21 figure [37T| evidence for an evolution of several of the dHvA 
frequencies near the metamagnetic transition, where in particular, the 1.8-1.6 and 4.2- 

4 S. A. Grigera, private communication. Complete mapping of the metamagnetic transition in 
the field-angle plane, as seen from various physical properties of S^Ra^Ov. 

5 The in-plane anisotropy of the critical field was never investigated in detail; we know from this 
work and from that of A. Rost (in preparation) that differences exist in the a6-plane between 5.1 
and 5.4, which indicates that g xx ^ g yy . 
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Figure 4.5: Field dependence of the oci pocket. The colours correspond to the intensity 
of the Fourier transforms (see figure I3.7P and the white circles to the maximum of the 
peak near 1.8-1.6 kT for each field value, using 360 windows of width 0.005 T _1 . 
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Figure 4.6: Three scenarios that can lead to changes in measured dHvA frequencies. 
Blue denotes energies and dHvA frequencies related to spin down quasiparticles, red 
to spin up and black to zero field values. 
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4.0 kT frequencies change by 0.2 kT (shown again in figure 14.5} for the 1.8-1.6 kT 
frequency). These changes lead to anomalous interference patterns in the amplitude 
of the 1.8 kT frequency in the low field side of the transition, as we saw in section 
13.3.3} figure 13.131 Several physical phenomena may produce such changes in dHvA 
frequency, and we analyse here three of these. We conclude that, as was suggested 
by A. Tamai et al., charge transfers indeed occur in Sr3Ru207 at the metamagnetic 
transition, probably due to a peak in the DOS, which means that the Fermi levels 
changes in that region of the phase diagram. 

The first case that we study here is that of a field dependence of the g factor, 
sketched figure I4TB1 (a). In metallic systems, changes in the energy coupling between 
the paramagnetic moment of the electrons and the magnetic field may arise, which 
produces non-linear changes in the magnetisation, sketched in the top row of figure 
14.61 In such a situation, Zeeman splitting in a normal dispersion accelerates, produ- 
cing non-linear changes in the size of FS pockets, and consequently in their associated 
real dHvA frequency. However, one does not measure the real size of the FS pockets 
with dHvA but rather their back-projected frequency, as described in section 11.3.51 
Equation II .281 means that one measures the real frequency minus its linear part, which 
corresponds to the phase of the oscillations, 

dF 

Fobs = F — B——. 

dB 

One therefore observes peaks in the frequency, described in figure 11.121 In the case 
of a field dependent g factor, these peaks appear with opposite sign for the spin 
species, or in other words, one measures both a sharp peak and a sharp dip. This is 
however not what was observed in the data, as only a positive increase in frequency 
was measured. 

The second case which leads to evolutions in the frequency is that where anomalies 
are present in the dispersion, depicted in figure POU (b). In that situation, the anom- 
alies do not likely appear symmetrically above and below the Fermi level, and it is 
probable that only one spin split pocket undergoes a non-linear change. In that case 
the dHvA peak becomes split and the measured frequency features a peak, positive 
or negative, for one of the spin species, and none for the other, which remains equal 
to the zero field frequency. Moreover, the frequency featuring a peak also possesses a 
field dependent quasiparticle mass, since the gradient of the dispersion changes with 
energy value, and the phenomenon of spin dependent mass is observed. In our data, 
no evidence was found for strong spin splitting with one of the peaks featuring a field 



6 The Fermi level here means the energy at mid-point between the Fermi levels for the up and 
down spin species. 
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dependent masgB Since the change in frequency between the low and high field side 
is of about 0.2 kT, much more than the peak width, such a phenomenon would have 
been detected if it had been present. However, if anomalies are placed symmetrically 
around ep, one should observe a behaviour similar to that of the field dependent g 
factor, a possibility we have already eliminated. 

The third case, which we think is of relevance for Sr3Ru207, is where a change 
in Fermi level produces non-linearities in the frequencies, shown in figure 14.61 (c). 
According the ARPES data of A. Tamai et al. presented in the previous section, a 
vHS was found in the vicinity of the 72 pocket, which produces peak in the DOS below 
the Fermi level. If the 72 pocket suddenly increases its share of the DOS, the required 
charge will be taken from the whole system, resulting in, to first order, an increase of 
the Fermi level. When measuring the frequency of a band other than 72, one should 
observe that one of the spin split frequencies features a non-linear increase, while the 
other one exhibits a slowdown, even though the dispersion does not possess any sharp 
change in gradient. Consequently, both measured frequencies possess peaks, but not 
necessarily of the same magnitude, as this depends strongly on the band structure. 
For electron-like bands, these peaks will be negative (as in figure 14.61 (c)), while for 
hole-like bands, they will be positive, as in the data of figure 14.51 

We did not measure the field dependence of the 72 frequency, but we did for all 
the other bands. We believe that these possess normal dispersions since we measured 
that their quasiparticle mass does not change. For some of these, a\ and 0:2, non- 
linear changes in the real frequency arise for both spin species, of the same sign. This 
third model is consistent with our observations in figure 13. 1\ where only a positive 
peak in the frequency was measured for both the a\ and frequencies, at 1.8-1.6 
and 4.2-4.0 kT respectively. We conclude that electron transfers occur between these 
bands and 72. 

However, in this simple picture a change in Fermi level should affect all the bands, 
and in particular, electron bands should feature dHvA frequencies that increase across 
the metamagnetic transition. The data does not seem to suggest that, but instead 
the 71 and S pockets vanish at c-axiqf|. It is possible that we have misinterpreted high 
field side data, and that these bands simply possess very different frequencies there 
compared to the low field side. We argue that the metamagnetic transition might 
affect some bands more than others, such that a%, a<i and f3 remain essentially intact 
while the other two evolve more dramatically. 



7 No strictly horizontal lines of intensity are present in figure [3771 at 1.8 or 4.2 kT. 
8 However, 71 and S seem to appear again at higher angles, see figure 13.111 and the discussion in 
section 14.3.21 
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4.2.3 Interference patterns from non-linear magnetisation and 
spin splitting 

Anomalous beat patterns may occur in dHvA when frequencies are not constants of 
the magnetic field. In such a case, the interference occurs between spin split bands, 
and the resulting pattern is more complex than the usual "spin zero" structure (eq. 
11.271 section 11.3. 5p . In our case of study, we believe that a change in Fermi level 
leads to non-linear changes in the fc-space size of the spin split bands, and that this 
evolution produces the intricate interference that we observed in the 1.8-1.6 kT band 
(the q;i band) calculated in section I3.3.3[ shown in figure 13.131 We present in this 
section a simple model that reproduces some of the observed features, but that fails to 
reproduce some others. This indicates that the real physics is more complicated than 
that which is presented here. However, we believe that comparison of the model and 
the data, along with the analysis shown in the last section, provides partial evidence 
for theories of metamagnetism involving a peak in the DOS. 

Interference in dHvA can be written in the following manner, using subscript 
symbols of plus and minus to represent spin up and spin down: 



where Fq represents the zero field frequency, AF± the respective changes in frequency 
for the spin species, and 9 is the angle between the field and the c-axis. The sum of the 
signal for both spins results in an interference factor which involves the difference in 
frequency changes AF + — AF^ B We are interested in modelling the field dependence 
of AF±(B) which correspond to changes in electron density of the up and down spin 
sheets of specific FS pockets as the system approaches the metamagnetic transition. 
If these were simply linear functions of B as it is the case with simple spin splitting, 
one would obtain ordinary spin zeros, but we know from the data that this is not the 
case. 

In our model, the 72 band crosses a putative vHS at the metamagnetic field, 
producing a peak in the DOS, which affects all the other bands through a change 
in Fermi level in order to conserve the number of quasiparticles. Essentially, if the 
peak is located below the Fermi level (as proposed by A. Tamai et al.), all bands 
are required near the metamagnetic field to donate a certain amount of holes, raising 
the Fermi level. Consequently, hole bands like the a± and a 2 undergo a reduction 

9 Note that AF + is negative. 




(4.5) 
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of size without changes in their quasiparticle masses, since their band structure does 
not feature sharp changes in their gradient. This produces the changes in frequency 
discussed in the last section, where both spin up and down electron densities for these 
bands are affected. 

Figure W7f\ left panel, illustrates the situation presented here. In two-dimensional 
materials, a peak in the DOS due to a vHS follows a logarithmic form [66_|l°, 



D( 



M log 



.1 



e - e c 



Taking the zero at the zero field Fermi energy, in such a model, when one applies 
a magnetic field the change in the number of quasiparticles with up (AN + ) and 



down spins (AAL), shown as blue and red areas in the left graph of figure fl~Tl must 
conserve the total number of particles, AiV_ = —AN + . The Fermi level, situated 
midway between e+ and e_, cannot remain at its zero field value but must move 
upwards as the magnetic field increases, in a way which is complex to calculate. The 
number differences are therefore difficult to evaluate since the limits of integration 
are not easily found. 

However, if one does not conserve the total number of particles, one has number 
differences, denoted with a star, which are not equal in absolute value, AN^ ^ —AN^, 
but the integrals for the number differences are known, a situation shown in the right 
hand side panel of figure 14.71 One finds AN± by evaluating the following integrals of 
the DOS, 



AN! 



±y B g{e)B 



D(e)de 



^li B g(Q) log A 2 - (^(i B g{6)B ± e ( 



=F(e c log e c - e c ) 



log \-fi B g{0)B ± e c | 



(4.6) 



AN+ features a superlinear rise near .5u. B g(6)B = e c , where g{9) is the orientation 
dependent g factor discussed earlier, while ANZ does not. The magnetisation corres- 
ponds to the difference between the two, M oc AiV* — AiV* , 



M(B) = A ± J2 



-u B g(9) log A 2 



-u B g(6)B ± e ( 



log \-ix B g{9)B ± e c 



10 The form of this function is not critical, but it must meet two criteria: it diverges at e c but its 
integral does not. A function of the form Ai\e — e c \~ A2 , < A 2 < 1 could also be used, and produces 
similar results. 
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Figure 4.7: Le/t Illustration of the DOS near a van Hove singularity. Number differ- 
ences for up and down spins {AN + and AN J) calculated when conserving the total 
number of quasiparticles are obtained from integrals of the DOS, shown in blue and 
red coloured areas, and the number differences are equal, but the limits of integration 
are not easily found. The Fermi level is required by particle conservation to move 
upwards as the magnetic field increases, shown with a dashed line. Right Similar 
illustration for the case where particle conservation is not required. In such a case, 
AN+ 7^ AiV*, but the limits of integration are known, and AN± are simple to cal- 
culate. Inset Number differences calculated as functions of magnetic field from the 
integrals illustrated in the main graph. The total magnetisation is proportional to 

an* + -an*_. 



and also features a superlinear rise near e c , since only AiV*_ undergoes a sharp increase 
(see the inset of figure H~T|) . 

In order to construct a simple model, we have assumed that the 72 pocket produces 
the sharp increase in the DOS below the Fermi level, and when this is reached by 
magnetic field, the missing number of electrons are provided by all the other bands 
through a change in Fermi level. However, we have kept the total number differences 
of eq. I4.6I For these other bands (for instance, «i and c^), their frequency differences 
AF + — Ai<L are related to unknown fractions of the total number differences, AN a ±. 
We supposed that the change in Fermi level produces evolutions of these numbers 
such that 

AN a+ - AN a _ oc AiV; - AiV* . 
Such a statement means that the argument of the beat pattern factor is simply 

AF + — AF_ txM(B). 

It does not mean, however, that the share of the DOS produced by these bands 
(01,0^2) possesses a peak, as we know from field resolved mass studies that this is not 



Chapter St Discussion 



138 



the case. It is through the change in Fermi level and a share of the DOS that is not 
constant, but which evolves slowly, that these fractional number differences AN a ± 
undergo non- linear changes 

These assumptions are reasonable in regard to the data, which seems to indicate 
that the frequencies evolve in a way which is roughly consistent with this model (see 
figures 14.51 and 14.61 (c) of the last section, compared to the curves in the inset of the 
left panel of figure |4TT|) . However, this calculation does not take into account properly 
particle number conservation, and is not perfectly adequate. It also neglects particle- 
particle interactions, which could lead to a discontinuous jump across the peak in 
the DOS, corresponding to a first order jump in the magnetisation. An appropriate 
calculation is outside the scope of this thesis, and the reader is referred to the work 
of Binz and Sigrist for a model that respects particle conservation [66] . 

Figure 14.8} top graph, shows a surface of the argument 



(AF+ - AF_)/Bcos9, 



calculated with this model. The nodes (or zeros) of the interference can be calculated 
by drawing a contour plot of this expression at values of n/2 + nn, n an integer. As 
level curves do not cross but close onto themselves, the same properties are associated 
with these beat pattern nodes. Moreover, these oscillations undergo a drastic change 
at \g{6)iiBB = e c , but also, since F(B) is highly non-linear, we expect that the beat 
patterns accelerate as the system comes near the metamagnetic transition on the low 
field side, as was observed in the data. Note that the modulation produced by the 
corrugation of the Fermi cylinder is independent of that resulting from non-linear 
magnetisation since, as one can see from eq. I1.23[ it can be expressed as a separate 
factor to the oscillatory term. 

Figure 14.8} bottom left graph, shows the result of the calculation of the beat pat- 
tern due to non-linear susceptibility combined with the anisotropy of g, using values 
for A x and A 2 of 0.4 kT and 80 T, and the DOS in the form A l log \A 2 /(.5g(9)B-B c )\. 
It produces a pattern where the original spin zeros curve towards the horizontal at 
high fields, and come closer together at high angles. The region at high fields where 
the zeros change curvature corresponds to the metamagnetic transition, along the 
line .bgnsB = e c . The middle left graph shows the zeros of the calculation plotted 
on top of the data for the rotation towards [110], at 1.8kT, One can see that this 
model produces the appropriate behaviour at high angles^. However, near c-axis, 



n If the structure of these bands was perfectly parabolic, and their share of the DOS constant, 
there would be no difference between the up and down spin dHvA frequencies, and therefore no beat 
patterns, even though both would feature a non-linear change near 8 T. 

12 Unfortunately, the data is highly under-sampled. One would need measurements performed 
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Angle ( deg ) Angle ( deg ) 

Figure 4.8: Top Surface plot of (AF+— AF-)/B cos 9 as a function of field B and angle 
9. Middle left Amplitude of the 1.8 kT frequency for the rotation towards [110], in 
colour, overlaid with beat pattern nodes, solid black lines. Bottom left Simulation of 
the beat pattern produced by the model of this section. Middle right Amplitude of 
the 1.8 kT frequency for the rotation towards [110] for comparison purposes. Bottom 
right Same simulation of the beat pattern, deliberately under-sampled with similar 
angle steps as the data. 
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the agreement not good, and on axis, there is no resemblance at all between the sim- 
ulation and the data. In order to convince the reader of the qualitative agreement 
of the model and the data at high angles, the bottom right graph presents the same 
calculation as in the bottom left graph, using a sampling identical to that of the data, 
shown middle right graph. The data and simulation on the right of the red lines agree 
qualitatively, but do not agree on the left. The effect of sampling is to repeat the 
modulation periodically with angle, as in a moire pattern, which is what appears in 
the dat£0 

We know from figure l3TT4"t section l3.3.3[ that the low angle part of the data (left of 
the red lines in figure |4"T8"| top right) is well reproduced by a simple warping parameter 
k i of 4.4xl0 6 m _1 , except very near the metamagnetic transition, but not at high 
angles, where the Yamaji angle region of high intensity is not observed. Combining 
the two contributions to the beat patterns, warping and spin splitting, does not 
improve the agreement, since in that case both beat patterns are superimposed over 
the whole field and angle range. The modulation produced by simple warping, given 
by eq. 11.231 (p. [31]), originates from a factor independent of Zeeman splitting, and 
does not influence the position of the zeros produced by the effect discussed herfi 
Consequently, this simple model does not fully explain the anomalous beat pattern 
observed in the data. 

The angular dependence of the g factor in this model produces a simple displace- 
ment of the peak in the DOS as a function of field. Since the anomalous modulation 
is of less intensity near c-axis, we propose that in this scenario, the peak in the DOS 
changes its width with angle, such that the magnetisation begins its superlinear rise 
at a higher field interval from the critical field at high angles than at c-axis, where 
it becomes very smalj^J Such a model is difficult to construct, and is beyond the 
scope of this thesis. We nevertheless conclude that the combination of an anisotropic 
g factor, a peak in the DOS and non-linear spin splitting is the most probable origin 
of the anomalous beat pattern of the 1.8 kT peak in field and angle, and that the 
same can probably be said of the other frequencies in the low field side. 

with finer angular steps, which was clearly out of reach with the time that was allowed for this 
experiment. 

13 Note that the damping of the intensity in the data as a function of angle is not correctly 
reproduced in the simulations, but is related to the experimental system and is not the focus here. 

14 Even when adding Zeeman splitting into the value for fcoo used in the calculation of the simple 
warping, an unrealistically high amplitude for the splitting has to be used in order to significantly 
change the beat pattern. 

15 Less than .5 T at c-axis compared to 3 T at 45°, see figure [37131 p. 11091 
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4.3 High field Fermi surface 

The dHvA effect at fields higher than about that of the metamagnetic transition 
possesses an enhanced complexity compared to low field regions. Two phenomena 
take place: all the dHvA peaks split into a large number of frequencies closely spaced, 
but frequencies also appear and disappear as a function of field angle. We propose in 
this section a model for the explanation of the first of these two phenomena. The field 
dependence of the dHvA frequencies calculated in section I3.2.2[ figure I3.7[ suggests 
that the splitting of the peaks increase with magnetic field. This type of behaviour is 
generally indicative of the occurrence of magnetic breakdown, where electrons tunnel 
between FS sheets. However, we have not found any reasonable explanation for the 
anomalous angular dependence of the high field side frequencies, but we argue that it 
is magnetically driven structural evolutions that lie at the origin of the phenomenon. 



4.3.1 Magnetic breakdown 

Magnetic breakdown occurs in a material when the the energy spacing between 
Landau levels, Ae^, becomes comparable to the gap energy between two bands in 
a certain region of k space, e g [45J: 

4 

where ep is the Fermi energy. In such situation, electrons may jump between FS 
sheets, and the hopping probability Pjj is exponential, 

P H oc e- H */ H , 



where Hb is the magnetic breakdown field, given by 

Hs = (4.7) 
en €f 

Due to simple conservation of the number of electrons, the intensity of the signal from 
the original band decreases as magnetic breakdown increases, and the intensity from 
each orbit changes with the value of the magnetic field. As was derived initially by 
Pippard j6jj EE], if one takes p as the quantum mechanical wave function amplitude 
for an electron that hops between bands, where 

r> 2 ^a- 
r H = p oc e H i 
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Figure 4.9: Left Regions of the BZ where magnetic breakdown may be possible. 
Four different points, numbered 1-4 in the plot, are repeated by symmetry all over 
the BZ. Right Table of calculated gap energies between bands for the four magnetic 
breakdown points, and corresponding magnetic breakdown fields. The error on these 
numbers is of about ±20%. 

and q, the amplitude of an electron that does not hop, such that q 2 = 1 — p 2 , one 
can calculate the total probability of specific paths through k space by counting the 
number of inter-band hops m and non-hops n, and the amplitude reduction factor of 
dHvA amplitude for this orbit is 

R M oc d(ip) m q n , 

where d accounts for the degeneracy of the orbit, and m and n are even numbers. Since 
different paths possess different probabilities, one expects complexity to arise with 
magnetic field, until it is high enough that only the breakdown orbits are populated. 

In Sr 3 Ru207, we effectively have a situation where the complexity of the dHvA 
spectrum increases with magnetic field (see the field dependence of the spectrum, 
figure I3TT1 p. 19*81) . Figure l4~9| left panel, presents four regions in the low field side 
BZ of Sr 3 Ru207 where magnetic breakdown may potentially occur, which repeat 
all over the BZ due to symmetry. Using average Fermi velocities, calculated using 
average values of kp and m*, we evaluated the energy gaps at these points by taking 
approximate fc-space distances (Ak) and using the equation 

f ,4(V flW+ V # ,)4M(i^). (4.8) 

We converted these results into magnetic breakdown fields using eq. 14. 7\ and tabu- 
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Frequency ( kT ) 

a. 4.2, 4.26, 4.32, 4.38, 4.44, 0.064 

b. 12.1, 11.65, 11.43 

c. 1.80, 2.40, 3.00, 3.60, 0.60 

d. 4.2, 5.1, 6.0, 6.9, 7.8, 0.9 

e. 2.5 

f. 0.43, .48, .9, 1.32 




-jt ka jr -it ka jt 



Figure 4.10: Schematic possibilities of magnetic breakdown in SrsR^CV. Top 
right Most probable magnetic breakdown points, where quasiparticle could tunnel 
in between FS sheets. Top left Table of roughly calculated additional frequencies 
starting from low field side orbits (ignoring electron transfers), a. to f. Sketches of 
possible new quasiparticle orbits. 
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lated these along with the gap values in the table of figure 14.91 All the breakdown 
fields were found to be within experimental range The error value of 20% was 
obtained from the maximum and minimum possible values the Fermi velocities used 
in eq. 14.81 are expected to take. 

From the four breakdown regions of figure fl~9| a number of breakdown orbits were 
considered, shown in the main panel of figure l4TT0| using the low field side FS. Taking 
these into account, we were able to calculate approximate frequencies for the resulting 
peak splitting for various FS sheets, using geometrical considerations. Unfortunately, 
these do not take into account the electron transfers occurring at the metamagnetic 
transition, since we do not know with certainty the details of the high field side FS. 
The table in figure 14.101 shows the results. 

The first situation we considered is that of panel a., figure l4TT0| where electrons hop 
from bands a 2 to It gives rise to eight possible hopping points, and consequently 
four regions where the orbit can be expanded. If a is the original k space area and 
6 the additional area due to magnetic breakdown, we expect to find, to first order, 
that the «2 pocket gives rise to five orbits at a, a + 6, a + 26, a + 36 and a + 46, with 
amplitudes corresponding to q 8 , p 2 q 6 , p A q i ) p 6 q 2 and p 8 . However, it is possible that 
the small orbit of area 6 is observed by itself, or that it is added in degenerate fashion, 
such that areas of a + 56, for instance, are possible. The number of ways equivalent 
areas are generated can be high in some cases, for example with a + 26, where it 
can appear at probabilities of p A q 4 , p 4 q 6 and p 4 g 8 , with degeneracy of 6, 12 and 10 
respectively. This degeneracy may enhance the amplitude of some of the possibilities. 

The complete list of possibilities is of course infinite, but we were able to identify all 
orbits of area up to a + 86, and 278 different paths were counted, in 46 different values 
of probability, for 9 possible areas, excluding the orbit of area 6 and its harmonics. In 
this specific case, 6 is of size of about 0.06 kT, and frequencies between 4.2 and 4.7 kT 
are expected. If we assume that the electron transfer of the metamagnetic transition 
reduced the size of «2 to around 3.9 kT, then we expect frequencies between 3.9 and 
4.4 kT. In dHvA experiments, frequencies between 3.9 and 4.7 were observed, which 
gives credibility to this model (see figure I3T71 p. |9"5|) . 

The second situation we analysed is that depicted in panel 6. of figure I4.10[ where 
an orbit similar to the /3 pocket of Sr2Ru04 is recovered, which lies around 12 kT, 
along with small variants. These were not observed. Next, in c, we have a situation 
similar to that in a., which gives rise to multiple frequencies between 1.6 and 4.0 kT 
spaced by around 0.6 kT, which could explain the peak seen near 3.2 kT, which in 



16 The absolute values should not be considered accurate, as they were produced using average 
values for k and m*, while the absolute value of these parameters vary in /c-space. More accurate 
values could be obtained using ARPES Fermi velocities, which was not possible during this project. 
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combination with the situation in a., could furthermore explain why this peak near 
3.2 kT is multiple. We do not see, however, other peaks between 4.0 and 1.6 kT, 
except at higher angles (see figure I3.10p . Another similar situation is that in d., 
where frequencies between 4.2 and 8.9 kT were predicted, spaced by around 1 kT. 
Combined with the possibilities in a., could potentially explain the group of peaks 
near 8 kT. The situation in e. gives rise to one orbit of around 2.5 kT, which was not 
clearly observed, although it cannot be ruled out. And finally, hopping between 7x 
and f3 gives rise to frequencies between 0.43 and 1.32 kT, which are clearly observed. 

We conclude this section by stating that the situations depicted in a., a, d. and /. 
and combinations between these could potentially explain the full complexity of the 
dHvA spectrum at c-axis. In particular, the situation in a. is a very likely explanation 
to the complex splitting of the 4.0 kT pocket. Moreover, it seems quite accidental that 
magnetic breakdown appears to be only observed above the metamagnetic transition. 
We argue that it is actually not the case, and that it starts to occur in the low field 
side of the transition. One can observe in figure 13.71 that the 4.2 kT peak splits into 
several above 5 T, a phenomenon that is also visible in the spectra taken from the 
data interval between 5 to 6.5 T, shown in figure I33| top graph. 



4.3.2 Strongly angle-dependent amplitudes 

The amplitude of the dHvA signal in field and angle in the high field side possesses an 
intriguing feature that is difficult to understand, that some frequencies disappear in 
specific angular spans. In particular, we have that the 1.8 kT peak vanishes completely 
above an angle between 5 and 10°, and the 0.9 kT peak vanishes below a similar angle. 
Moreover, the 0.43 kT peak appears at an angle of around 15-20°, and other peaks 
appear briefly in the region between 2.5 and 4.0 kT at low angles. We also have that 
those for which we can observe the angular dependence do not follow the Fq/ cos 6 
form. These types of behaviour are not related to anything in the usual dHvA theory 
for 2D systems, and we do not present in this thesis any model to explain it. It 
remains to this point a rather mysterious phenomenon, but we suggest here a few 
ideas about its possible origin. 

The simplest explanation appears to be that angle changes the condition for mag- 
netic breakdown. However, its origin does not lie with a combination of Zeeman spin 
splitting and magnetic breakdown. The reason is that the effect of magnetic break- 
down is not directly related with Zeeman splitting 17 !. Magnetic breakdown stems from 
Landau diamagnetism, or the orbital motion of electrons in a magnetic field, without 
reference to their spin. In an imaginary world where charged fermions without spin 



''Many thanks to C. A. Hooley, who pointed that out. 
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existed, the dHvA effect could still be observed, whereas in a spin split band structure, 
there is no process by which the two spin species can interact except through spin 
orbit coupling. Consequently, the effect of the anisotropy of g discussed in section 14.21 
does not change with angle the condition for magnetic breakdown. 

Thus there must be another process which changes the magnetic breakdown field 
values, and we suggest two here. The first and simplest is that the gap between 
bands may possess a k z dependence, which can be smaller in certain regions than 
others, which at certain angles, would become by coincidence extremely broken down 
extremal orbits. In such a situation, one may find that new orbits appear briefly 
as a function of angle, and that original ones disappear as well, both of which have 
been observed. The second process we suggest, which appears the most probable, 
is that the strong magneto-structural coupling in Sr 3 Ru 2 07 produces small changes 
in the crystalline structure with field angle, which distorts the FS itself [621 ED]- 
Consequently, the gaps between bands may change as the FS evolves and thus change 
the condition for magnetic breakdown. 

Finally, another process has been suggested by S. R. Juliarl^l where stress induced 
changes in samples occur due to the experimental mounting procedure, an effect which 
we have not been able to rule out due to time constraints. In this scenario, the layer 
of silver epoxy used for the heat sinking contacts could induce, with magnetostriction 
effects at the metamagnetic field, inhomogeneous stress across the samples. Due 
to the strong magnetostriction properties of Sr 3 Ru 2 7 , these could deform the FS 
inhomogeneously and produce various dHvA frequencies from different parts of the 
sample. We do not believe that this is possible, due to the very thin nature of the 
epoxy layers. Moreover, as we mentioned earlier, magnetic breakdown occurs not 
only above but also below the metamagnetic transition, indicating that it is probably 
not related with the metamagnetic transition and magnetostriction occurring there. 
We could not however rule out this possibility based on the current data. 

4.4 Absence of quasiparticle mass 

enhancements as seen from dHvA 

The absence of an enhancement of the quasiparticle masses in Sr 3 Ru207 is an ex- 
tremely unexpected phenomenon. Not only since it was reported otherwise in the 
past [40], but also through the fact that other metamagnets, CeRu 2 Si 2 , CePd 2 Si 2 , 
CeRhIn 5 [56], [7TJ [72] do feature mass enhancements and Sr 3 Ru 2 7 now figures as 
an exception. Secondly, as we saw in section ll.3.2[ in 2D systems one can calculate 

18 Private communication. 
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the electronic specific heat coefficient by summing the quasiparticle masses, and this 
should be true at all fields. We described in section 11.2.41 that the electronic specific 
heat of SraRu207 is greatly enhanced near the QCEP, and so should the sum of the 
quasiparticle masses be, a fact that was confirmed again recently down to 100 mK by 
A. Rost [62J. Thirdly, in a Fermi liquid the A coefficient to the quadratic temperat- 
ure dependence of the resistivity is, as we saw in section ll.l.l[ eq. I1.2[ proportional 
to the square of the quasiparticle mass. This parameter was also reported strongly 
enhanced (see section 11.2.40 . 

We are facing an important discrepancy between measurements if we assume that 
we have determined a complete set of FS sheets and know all of the FS quasiparticle 
masses. However, we cannot make such a statement with certainty, but can only in- 
stead present very stringent constraints to the origin of the "missing" mass. Moreover, 
in order to strengthen our position, we are required to investigate how it came to be 
that the opposite conclusion to ours was inferred in the past. We thus present here 
first an explanation of how a divergence of the quasiparticle mass may be erroneously 
deduced from field dependent low signal to noise ratio and the use of non-linear LK 
fits with three parameters. We discuss here the process by which it can happen, and 
a give a detailed numerical analysis in appendix [Fl We moreover show there that our 
quasiparticle mass extraction method is able to well detect an enhancement of the 
mass as long as the signal to noise ratio is larger that 0.5, a condition which is fulfilled 
for all the bands that we analysed in section 13.41 Finally, we discuss the implication 
of our result and argue that the 72 pocket may possess the missing mass. 



4.4.1 Source of systematic errors in previous work 

Previous dHvA work on Sr 3 Ru207 by Borzi et al. featured the use of three parameter 
non-linear LK fits to the temperature dependence of the dHVA amplitude Qj. We 
suggest here how such a procedure can generate an increase of the quasiparticle mass 
when the signal to noise ratio is low and varies systematically with magnetic field. 
The signal one has in such a case is of the form of a modulated oscillating function to 
which random noise / is added. In a hypothetical case where only the LK function 
modulated the oscillations, one has that 

M(X) = LK(X, T) sin(FoX) + f(X } T). 

In the normal analysis procedure of dHvA, in order to calculate the oscillatory amp- 
litude, one performs a Fourier transform over an inverse field width AX, centred at 



19 For a definition, see eq. 12.81 page [53J This information was obtained through a private commu- 
nication with R. A. Borzi. 
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Figure 4.11: Two illustrations of a situation where a three parameter fit is used with 
noisy data and results in a systematic error for the mass. Top Realistic situation 
with only a few data points. The data without noise is shown in black squares, 
with associated LK function of mass m* = 12m e at a field of 12 T, shown with a 
solid line. In red circles appears the absolute value of noise created with a random 
number generator of gaussian distribution, and in green triangles, the absolute value 
of the sum of the LK function and the noise. In solid green appears the best three 
parameter LK fit, with a spurious enhanced mass of m* = 18.5 ± 2m e . Note that 
the same temperatures as those in the experiment of Borzi were used. Bottom left 
and right Same calculation using 500 uniformly distributed data points and the same 
noise level. On the left is shown the original LK function in black, the noise in red 
and the sum of both in blue. On the right is shown the LK function, in blue, the 
absolute value of the sum of the noise and the LK function, in black, and the three 
parameter LK fit, in red. 
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M{F) = LK(F, X , T) + f(F, X , T). 

Without noise, one integrates the power spectrum, or modulus square, of LK over F 
in the region where LK(F) takes the shape of a peak, and take a square root, leading 
to a value of LK proportional to (LK(X)), an average of LK(X) over the range AX, 
centred at X H With noise added, one has 



'J \M(F, X )\W = | (LK(X, T) + f(F, X, T))\, 

which is the average of the sum of the noise and the amplitude function, in absolute 
value. This does not correspond to the average of LK(X, T) plus the average of 
/(X, T), and even less to the sum of the absolute value of the averages: 

\{LK(X,T) + f(F,X,T))\ < \{LK(X ,T)) + (f(F ,X,T))\ 

< \(LK{X ,T))\ + \(f(F ,X,T))\. (4.9) 

We stress that the result is smaller than just the LK function added to the modulus of 
random noise. As one performs calculations using the absolute value of the average of 
the sum of the signal and the noise, one finds that at high temperatures where the LK 
function is very small, the result is positive and non-zero, of which one tries to take 
account of by adding an offset to the LK fit in the three parameter scheme. However, 
the signal of amplitude LK is not simply added to a constant, as we showed in eq. 
14. 9[ but to noise which may be positive or negative and the average and absolute 
value are taken after the sum. As noise increases in magnitude, the LK function 
becomes "buried" into noise, only the low temperature part rising above. Using a 
third parameter in the fit appears to work, but leads to sharper LK functions than 
reality, which means the quasiparticle mass, inversely proportional to the width of 
the LK function, is estimated larger. 

Moreover, (LK(X,T)) and (f(X,T)) are averages of the originals over a width 
AX (see section 12.2.21 for details), such that the noise function f(X ,T) = (f(X,T)) 
is not completely random, but has a correlation in X of length AX. It should not 
be surprising, then, to find that when the noise is large compared to the signal, the 
deviations that one measures on the mass are smooth and continuous curves as a 
function of inverse field. 

Figure 14.111 presents two numerical examples of this systematic effect, using a 
quasiparticle mass of 12 m e at a field of 12 T. The noise was created with a random 



•'Using Parseval's theorem. 
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number generator of normal distribution, and had an rms amplitude of about 20% of 
the signal at the lowest temperature. In the first example (top graph), the temper- 
ature values were the same as in the experiment by Borzi. We found that when we 
adjusted the three parameter LK function to the absolute value of the data added 
to the noise, in this specific instance we obtained a mass of 18.5 ± 2 m e . When we 
performed this calculation several times by regenerating the noise, we systematically 
obtained masses between 1 and 1.5 times larger than the original value of 12 m e . 
This is changes with differing noise amplitudes, and the systematic error on average 
becomes larger the smaller the signal to noise ratio. In the second case (bottom two 
graphs), 500 data points were used with the same scenario. An enhanced mass of 
14 m e was obtained. Hence, combined with the fact that the noise has a correlation 
length, the field dependent mass calculated with three parameter fits may feature an 
apparent continuous enhancement if the signal to noise ratio is low and has itself a 
field dependence, a condition which was fulfilled in the previous work of Borzi et al^i 
Appendix |F] presents detailed calculations that we performed in order to verify 
the self-consistency of our calculation method for the quasiparticle masses, and to 
demonstrate how the masses can feature enhancements following the argument of 
this section. The method we used was to simulate dHvA oscillations with noise 
added to them, and extracted the quasiparticle mass as a function of magnetic field 
in the same way as was done in section 13.41 We found that our calculation method is 
appropriate, and that we effectively produced a continuous mass enhancement when 
using three parameter fits applied to data where a low signal to noise ratio evolved 
with magnetic field. 



4.4.2 Interpretation of the absence of an enhancement 

We present here strong constraints to parameters of a putative band, from the view- 
point of the dHvA experiment, that possesses the missing, enhanced quasiparticle 
mass that we require in order to reconcile the field dependent specific heat and the 
sum of the dHvA quasiparticle masses. It must possess parameters that make it diffi- 
cult to measure with second harmonic dHvA, and there are four possibilities. Firstly, 
this band could possess such a large in-plane k space area that the Dingle reduction 
factor damps it enough to escape detection. Secondly, it can have a very large qua- 
siparticle mass even away from the metamagnetic transition, where it is expected to 
diverge, and in that case it is the LK reduction factor that decreases its amplitude. 
Thirdly, it can be so small in area that, due to our use of second harmonic detection, it 
is the Bessel function modulation factor that damps its low frequency amplitude. Fi- 

21 This information was obtained by an analysis of the original data of Borzi et al.. Note that we 
reproduced the reported enhancement when using three parameter fits with the original data. 
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nally, the field dependent missing mass observed in specific heat and resistivity could 
be caused by exotic non-fermionic degrees of freedom, not by the band structure, and 
so would be invisible to dHvAcf. 

In the first case, one can see in the ARPES data of figure 14.21 that it is difficult 
to place an orbit in the BZ that is much larger than 4.2 kT, which we detected with 
good signal to noise ratio. Part of this large band would need to be situated in the 
narrow region between ai and 71, which according to ARPES dispersion data along 
the TM line 



3, no evidence was found for such a band. 
In the second case, a very high, field dependent mass for a single band would be 
consistent with the scenario where 72 is not included in the FS. In that case, near zero 
field this band would be required to possess a quasiparticle mass which accounts for 
half of the total specific heat, of around 74 m e . However, such high masses have been 
observed in the past with the Cambridge system 3 and if this band was small, it 
would give rise to an extremely flat dispersion near the Fermi level, which would most 
probably have been observed by ARPES. From ARPES data, there is no evidence 
of any other band than 8 located at T within the area covered by a±. Finally, if 
such a band was to be located elsewhere than in the centre of the zone, by symmetry 
argument, it would be required to be present in in two or four equivalent locations, 
depending whether it is situated at the M point or elsewhere. If it was at the M 
point, the required mass per part of the orbit would be of 37 m e , but no evidence 
for a corresponding very flat dispersion was found there. If it was located elsewhere, 
then symmetry requires it to be present four times and the required mass per part of 
the orbit would be around 19 m e , low enough to be observed. 

Consequently, if one dismisses theories involving non-fermionic degrees of freedom, 
all evidence points by elimination towards 72 being part of the FS and possessing the 
missing mass enhancement and the mass responsible for half of the zero field electronic 
specific heat. As stated earlier, band structure calculations by Mazin and Singh [37] 
predict that this band should be bilayer split, and consequently, counted eight times. 
Moreover, the finding by A. Rost, presented in section 14.1.11 of a quantum oscillations 
peak around 110 T with a quasiparticle mass of about 9 m e accounts for nearly the 
right number of degrees of freedom, 72 m e when counted eight times. A. Tamai et 
al. also reported that a vHS is present in its vicinity, which could be reached by 
Zeeman splitting in a magnetic field, and lead to a peak in the DOS. In conclusion, 
we argue that it is the small 72 band observed in ARPES that holds the key to the 
reconciliation of our measurements to the field dependent specific heat. 



22 Such theories are beyond the scope of this thesis and are not discussed. 
23 A. Tamai et al., private communication. 

24 



See for instance the work of Taillefer et aL[75] 
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This model does not, however, explain the divergence of the A coefficient of the 
power law of the resistivity, measured by Grigera et al. [10]. Even though this para- 
meter is proportional to the square of the quasiparticle mass, in a multi-band mater- 
ial, one obtains the resistivity by adding the conductivity of the various bands and 
inverting the result. In other words, the total resistivity is added as in a parallel 
circuit, 



where k refers to specific bands. A divergence of A for one band does not show as a 
divergence of A in the total resistivity. It is therefore required that the A coefficients 
for more than one band diverge simultaneously. We do not hold any reasonable 
solution to this conflict between measurements. 

4.5 dHvA in the nematic phase 
4.5.1 Interpretation 

The observation of quantum oscillations inside the nematic phase was an extremely 
interesting discovery which provided us with important new information. The amount 
of established facts about this phase of the electronic system in S^R^O^ is still small 
at this point, and the observation of dHvA may help moving towards an explanation 
of the phenomenon. 

The most important fact that this measurement established is that the phase is 
most probably composed of Landau quasiparticles. Essentially, quantum oscillations 
with an LK temperature dependence stem from basic Fermi liquid theory, where 
fermionic charged particles follow the the Fermi-Dirac distribution function, and the 
LK function originates from the Fourier transform of the derivative of this distribution 
(see section H.3.2[ where the calculation is given). 

These oscillations also told us more about the nature of the electronic structure 
of this phase. We observed two quasiparticle orbits, at 1.0 and 2.5 kT, with masses 
of around 6-7 m e . We have argued in section [3751 that these frequencies, due to back- 
projection issues, corresponded to 0.9 and 1.8 kT, two of the low field side orbits, of 
which the masses were not entirely inconsistent. This suggested that although we did 
not observe any other frequencies in the nematic phase, its Fermi surface appeared 
to possess some similarity with that of the low field side. We therefore argue that 
the nematic behaviour in S^R^O^ is produced by a specific subset of the low field 
side FS. Such a statement is in accordance with the suggestion by A. Tamai et al. 
from ARPES (section 14. 1 .3j) measurements that the metamagnetic behaviour could 
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be produced by the 72 pocket or the saddle point in the dispersion situated nearby 
[36]. 

Moreover, one can extract a length scale from these frequencies, since the area 
within a k space orbit is proportional to that of the real space orbit through the 
factor (eH/h) 2 . From eq. I1.9[ the frequency of 1.8 kT corresponds to a fc space 
orbit of average radius of 3.57xl0 9 m -1 , which through eq. I1.8[ in a field of 8 T, 
equals to an average real space diameter of approximately 580 nm. This value may 
constitute an upper boundary for the size of any real space orbit, and a lower boundary 
for the distance between nematic phase domain walls. Moreover, the interaction of 
the domain walls with the cyclotron orbits may explain the extreme reduction in 
amplitude of the oscillations between the low field side and the nematic phase. 

Finally, the strong angular dependence of the amplitude of the quantum oscilla- 
tions in the nematic phase is an aspect that we do not understand at this point. The 
similarity with the behaviour at higher fields than the metamagnetic region reminds 
us that we face the same problem in these two situations: we do not know of a process 
that damps the amplitude as a function of angle as dramatically as what is observed. 
But we may still suggest a process that could have this result, related to the aniso- 
tropic transport properties that are triggered in Sr3Ru207 by the application of a 
small in-plane magnetic field. If, as was suggested by Grigera et al. the nematic 
phase domain walls are aligned in some way by the in-plane magnetic field, resulting 
in anisotropic transport, the distance between walls could also be affected. If that 
distance is reduced, it may be more difficult for the system to fit quasiparticle orbits 
anywhere in real space, and lead to suppression of dHvA. 



Conclusion 



We review in this conclusion the facts that have been established in this thesis regard- 
ing the FS near the QCEP in S^R^Oy. We had set out to study three aspects of 
this problem: the low and high field FS, that in the nematic phase and the evolution 
of the quasiparticle properties in the region surrounding the QCEP. We have, in col- 
laboration with ARPES experiments by A. Tamai et al. [36J, determined with some 
degree of confidence the topology of the FS at low fields. We discovered that dHvA 
oscillations exist inside the nematic phase, and finally, we established that none of the 
orbits for which we were able to obtain high resolution dHvA signals undergo quasi- 
particle mass enhancements near the QCEP. We present here a brief unified picture 
of the physics of the FS of Sr 3 Ru20y based on what is known at the present time. 

The dHvA experiments at low fields revealed the proof that the FS of SrsRu207 
is quasi-two-dimensional. We showed that all the FS pockets except one exhibited 
a behaviour consistent with a cylindrical topology connected to the top of the BZ. 
Following this, A. Tamai et al. suggested a model of the in-plane FS that is consistent 
both with ARPES and our dHvA data, respecting measured cross-sectional areas and 
quasiparticle masses. One pocket was missing in our dHvA data, but was found by 
A. Rost in magneto-caloric quantum oscillations [62]. This model was extremely 
successful in reproducing the zero field electronic specific heat measured by Ikeda et 
al. [2T], but the Luttinger counting revealed a quarter of an electron missing over the 
16 required by the crystal structure. 

We moreover determined that a non-linear change in Fermi level occurs at the 
metamagnetic transition, such that some FS pockets change size abruptly, as inferred 
from dHvA frequencies. Through this frequency evolution, we found that the meta- 
magnetic transition influenced the dHvA oscillations at low fields, starting from 3 T 
away from the critical region at high angles and from 0.5 T at c-axis. The gradual 
difference arising in the number of electrons between the spin species of the ot\ pocket 
produced an additional interference pattern to that inherent to its corrugation along 
the k z direction. We were not able to reproduce fully this pattern, but we performed 
a simulation that agreed reasonably well with high angle data, suggesting that the 
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effect gradually changed magnitude with angle, reflecting the fact that the new pat- 
tern modulated the data from a higher field interval down to lower fields at high 
angles compared to the c-axis direction. We interpreted that as a change of width 
of the peak in the density of states responsible for metamagnetism, which could be 
due to small evolutions of the crystal structure with field angle. This is effectively 
expected from the known magneto-structural coupling. We assume that the other 
frequencies present in the low field side possess the same property as the ol\ pocket, 
but their field and angular amplitude modulation was too complex to analyse due to 
the close proximity of different peaks, related to bilayer splitting or possible magnetic 
breakdown even at such low fields. 

We studied the FS at high fields as well, above the metamagnetic transition, and 
found surprising complexity in the spectra, at all angles. The frequencies were all split 
into a high number, much more than what is predicted by spin splitting due to super- 
linear changes in the magnetisation. As this complexity increased with magnetic 
field, we associated it with magnetic breakdown, which starts at fields as low as 5 
T in one case. We discovered another unexpected phenomenon, which is not fully 
explained yet, that the frequencies possessed highly unusual angular dependence of 
their amplitude, where some disappeared and others appeared at specific angles. We 
suggested that it could be related with small evolutions in the crystal structure with 
the value of the in-plane field, resulting in changes in the gap value between bands, 
modifying the condition for magnetic breakdown. In such a model, new orbits appear 
while others lose intensity, controlled by the different hopping probabilities included 
in each /c-space path. The frequency values also increased with angle, but were not 
completely consistent with the normal quasi-two-dimensional angular dependence. 

Turning our attention to the nematic phase, we made the interesting discovery 
that dHvA orbits exist within the first order phase transition boundaries. We found 
a few cycles of two frequencies which were consistent with the a± and 71 pockets. 
They possessed normal LK temperature dependence, and their quasiparticle masses 
did not agree perfectly with those of a 1: 71 but possessed nearly the right value. We 
moreover discovered that these oscillations exist only within 5° of the c-axis, and are 
rapidly suppressed at higher angles. 

The existence of quantum oscillations inside the nematic phase is a surprising 
result since the resistivity is highly enhanced and one expects high scattering of charge 
carriers. However, if the scattering is produced by domain walls, it is possible that 
only orbits enclosing real space areas of average diameter larger than the average 
distance between walls are suppressed by that process. We calculated that such a 
distance is of the order of 600 nm, comparable to the value of the mean free path of 
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300 nm. Moreover, we reflected on the fact that with an in-plane field, as the domains 
become aligned, the average distance between domain walls may decrease significantly 
and suppress all cyclotron orbits altogether, producing the observed strong angular 
dependence of their amplitude. 

The most important result of this thesis is related to quantum criticality and 
metamagnetism. We calculated the field dependence of the effective quasiparticle 
mass of five of the six orbits observed by ARPES, and discovered no enhancement 
near the QCEP, in complete contradiction with previous dHvA experiments [40J. We 
infer that quantum critical properties are not present in any of the bands observed 
in dHvA, putting strong constraints on many theoretical models. Moreover, since 
in two-dimensional materials, the sum of the quasiparticle masses must correspond 
to the electronic specific heat, this conclusion is in complete contradiction with past 
measurements [27J El] except if the enhancement is produced by the missing pocket, 
labelled 72 by A. Tamai et al.. It is also in contradiction with existing transport 
measurements for the same reason, since the A coefficient of the quadratic temper- 
ature dependence of the resistivity is proportional to the square of the quasiparticle 
mass, and features a strong enhancement [10J. 

We proposed that the specific heat enhancement is either produced by the 72 
pocket or by non-fermionic degrees of freedom. In the first case, we suggested that 
metamagnetism is produced by the van Hove singularity identified by A. Tamai et al. 
near the corner of the orthorhombic BZ, reached with magnetic field by the 72 pocket, 
which already possesses half the total electronic specific heat. It would also explain 
the electron transfers that we detected in the a.\ and pockets, as moving through 
a large peak in the density of states with field should force a change in Fermi level, 
as well as a non-linear increase in magnetisation. We hold, however, only indirect 
proofs of this inference, and cannot completely rule out other explanations related to 
non-fermionic modes. Direct quasiparticle mass measurements of the 72 pocket are 
required in order to decide which explanation is the right one, and they have not been 
possible during this project, the signal being too small to detect with dHvA. 

We end this conclusion by mentioning two types of experiments that could be 
done in the future to complete this work. The first class corresponds to those which 
are required to close all the uncertainties in the work to date, and the second to new 
experiments which could help gather more information about the 72 pocket. 

Three measurements should be done using the current dHvA system in order to 
clear a number of issues in the interpretation of the data that was obtained in this 
project. The first regards questions about the sample mounting procedure, and to 
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whether, as was suggested to us by S. R. Julian, it produces inhomogeneous strain in 
the samples, leading to the high peak splitting of the oscillations in the high field side. 
The second is to improve the knowledge of the thermometry or thermal equilibrium 
of our mounting procedure used in the Cambridge system, for instance by using 
another known Fermi liquid system using the same procedure. This should enable 
us to understand better the low temperature deviations of our data compared to the 
LK function. The third is to verify into even more detail whether a modulation field 
higher than that which was used by Borzi et al. did not suppress an enhancement 
of the quasiparticle mass, and would involve measuring the first harmonic dHvA at 
very low modulation fields. 

A number of new experiments could be done in order to determine with certainty 
whether 72 is responsible for the specific heat enhancement and metamagnetism in 
Sr 3 Ru207. Unfortunately, we have already used most of the methods that provide 
information related to specific branches of the FS, as only ARPES and quantum 
oscillations can resolve individual FS pockets. There exists still one type of quantum 
oscillations experiment that has not been performed yet, using torque magnetometry. 
However, the signal from such experiments possesses an amplitude that is proportional 
to the derivative of the dHvA frequency F as a function of angle 6, which equals zero 
at c-axis, the direction of interest. Consequently, one cannot study the 72 pocket at 
the QCEP, but valuable information could still be obtained near the metamagnetic 
transition, since it exists at all angles. 

However, more experiments could be done using dHvA. All our measurements 
were performed with second harmonic, the amplitude of which is proportional to the 
square of the dHvA frequency. The first harmonic could be used, since is it linear in 
dHvA frequency and better signal could be obtained for F < 1 kT. Time constraints 
have not allowed us to perform these experiments in this project. One could also 
perform dHvA using samples under pressure, which could tune the metamagnetic 
transition to a different field value. However, in both of these cases, the hope to 
resolve the field dependence of the quasiparticle mass is slim, as good signal to noise 
will be difficult to achieve. In that respect, experiments using the magnetocaloric 
effect could produce better signal than dHvA for small FS pockets if performed at 
lower temperatures than the recent measurements by A. Rost [62], using for instance 
the Cambridge dilution refrigerator. 

Other types of experiments could also provide more information about 72. In 
particular, ARPES experiments performed with doped samples of Sr 3 Ru20y may 
allow the observation of the behaviour of 72 as one crosses the peak in the density 
of states. But according to our model, the band responsible for metamagnetism 
is of hole type, which means that hole doping is necessary. This was confirmed 
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by J. Farrell, who doped electrons in S^Ro^C^ and observed an increase in the 
^magnetic fle.d vaiufl Hole doping has not been performed in Sr 3 Ru207 to 
our knowledge, and is probably difficult to achieve 26 !, but would provide an excellent 
tool to study metamagnetism and the QCEP, as it would bring the metamagnetic 
transition towards zero field. 



Private communication, to be published. 

It involves replacing the Sr atoms by one of the alkaline elements, which are extremely reactive. 
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Fourier transform of the derivative 
of the Fermi function 



This Fourier transform is more easily done backwards, by applying the back Fourier 
transform on the LK function, as a function of A, with reciprocal variable ip — <f>/\. 
The LK function is 



sinh(7rA) 
The back Fourier transform is 

This integral converges, since sinh(A) diverges faster than A for A — > ±00. One applies 
Cauchy's residue theorem to solve it, by setting A = Re 10 and performing an integral 
over a semicircle contour in the complex plane, the linear part on the real axis. The 
arc is chosen to lie either on the positive or negative half complex plane, depending 
on the sign of tp. The integral involves two terms, one for the section on the real axis, 
one for the arc in the complex plane. By taking the limit of an infinite radius R, the 
sum of these integrals becomes equal to 2ni the sum of the residues at the poles, 

roc nXe iX^ f 2 " nRe*?**"* n .v-/ ■ , \ 

/ — — — —-d\+ hm / — — — ^ ... Rat) — 2m > (residues), 
sinh(TrA) R-+ooJ smh{nRe ld ) ^ K h 

where the second term is zero everywhere. 

One is required to locate the poles and determine their order. The integrand 
diverges on the imaginary axis, where the smh(iirR) function becomes — sin(7ri?), 
except at R — 0, where the integrand is zero. The poles are evenly spaced, appearing 
at 

R = n, 
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with n an integer. Their order corresponds to that of the first term in a Laurent 
series, 

1 1 

+ ..., 



sin(7nR) nR 

and is equal to one. The residues are calculated in the following manner 

— iirRe~ R ^ 

lim (inR - irm) - — = -{-lYin-KeT^ 

R->n sin(7ri?) 

As indicated earlier, one chooses the side to take the integral along the arc that 
defines the sign of n, with that of ip, such that the argument of the exponential in 
the numerator of the integrand is always negative when on the imaginary axis. The 
result of the integral is 



-dX = -27rz^(-l) n zn7re-^ n| = 2vr^(-l) n r 



sinh(7rA) 

v ' n=l n=l 

Since this is an analytic function of ip except at zero, one uses its derivative in order 
to have a more tractable form: 

, oo 

/ = 27r^-y(-l)"e-^ n| . 
dib ^ v 

r n=l 

This is a geometric series, where setting 



results in 



n=0 

The integral becomes 



r = — e 



oo 

Y- = —. 

^^^^^ iyTl ly 



d 1 -7re~M 1 

/ = 2tt — rn = 2- rnTTT = 2vr- 



d^l + e-l^l (1 + e-l^l) 2 1 + cosh|^|' 
and the result of the Fourier transform is 

1 1 



9W 



V27T 1 + cosh 
the derivative of the Fermi function. 
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Arbitrary cross-sectional area of a 
cylindrical expansion 



We carry out in this appendix the calculation from the point it was left off in section 
ll.3.4[ where we required to calculate the cross-sectional area A(k z , 9, 0o) in order to 
perform the integral 



with a cylindrical expansion of A with parameters /i and v. The in-plane angle 0o 
was added in order to orient the in-plane shape with the direction of rotation of the 
magnetic field. The expansion of the Fermi wave vector is 



where k z corresponds to k z x 2ir/hf )Z , and ht> z is the height of the Brillouin zone 47r/c, 
with c the c-axis lattice parameter. 

We define two cylindrical coordinate systems, one with k z direction parallel to the 
c-axis, the other rotated by 9 in the direction of <po, as shown in figure IB.lt using 
prime notations for the tilted system. Using a cartesian system with the direction 
k x in the axis of rotation of the tilted system, we have a set of equations for the 
transformation between the two coordinate systems: 




(B.l) 




(B.2) 



/i=i,f=i 



k y cos 9 + k z sin#, 
—k y sin 9 + k z cos 9, 
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V 



Figure B.l: Schematic Fermi surface, with a plane intersecting at an arbitrary angle 
9 and height in k z . Two systems of cylindrical coordinates are shown, the normal 
one with k z , p and as variables, and the tilted one with primed variables k' z , p' and 
(f)'. Two infinitesimal areas are drawn, one in the normal cylindrical coordinates, the 
other in the tilted system. The infinitesimal units of area are pdpdcj) and p'dp'dip' . In 
order to calculate the area of the tilted ellipse using p' and 0, one has to rescale the 
area unit by cos 9. 

In relation to the cylindrical coordinates 

p 2 = k 2 x + ky, k x = p cos 0, k y = p sin 0, 
= *£ + *?, K = ^ = p'sin0', 

one can put these equations together to get a new set for k z and k' z : 

k z = p' sin 0' sin 6* + k' z cos9, 
k' z = — psin0sin# + k z cos 9. 

As shown in figure IB.1[ we are required to integrate the ellipse that corresponds to 
the intersection of the plane at constant k' z and the arbitrary shape of kp{k z ). We 
then evaluate the infinitesimal unit of area on the tilted ellipse that lies at constant 
k' z , pdpd<ft, using variables that are easier to work with. For that, we differentiate the 
last set of equations on this plane, where dk' z = 0: 

dk z = dp' sin 0' sin 9 + p' cos 0' sin 9d<p' 
= — <ip sin sin # + dk z cos 9 — p cos sin #<i0, 
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which, by eliminating dk z leads to 

dp sin tan # + d<pp cos tan = dp' sin 0' sin 9 + p' cos 0' sin 0d0'. 

The radius p' only varies when p varies as well, and when only the angles change, 
using the fact that k' x = k x = pcos0 = p'cos0', we obtain 

cos J ' 

which results in an area unit of 



p'dp'd<j>' - P ' df> 



cos 6* 

and enables us to perform the calculation using the variable p' . We now calculate the 
area of the ellipse: 

2tt fk F (k' z ) 

p'dp' 



A = f da' = f p'dp'dcf)' = — !— f d<p f 
J J cos 9 J J 

d(j)k F {k' z ) 2 



2cos# 

We replace k' z for its value in terms of variables of the normal system, substituting p 
by its average value k o, and multiplying by 27i/h BZ , 

k' z = k z cos 9 — Koo sin sin 9, 

and normalising k' z to a value such that k z varies between — 7r to ir, as in equation 
EH i.e. by cos 9% 

1 f 2n 

A = / d<pk 2 F (n z — psin0tan#) 

2 cos 9 Jq 

From here, we may check our calculation with a circular cylindrical Fermi surface 
kp = k Q, which appropriately yields A = nk 2 Q / cos 9. 

We carry on with the expansion of the square of kp to first order (dropping terms 
in k^ u k aT , and the orientation angle 0o for now): 

k 2 F (k z ,(j)) = (k 00 + ... + A;^cosz//t z cos/i0) 2 

= klo + ■■■ + 2&ooA;^i/ cos vk z cos p0 + ••• 
~ k^ + 2^oo k^ u cos i>k z cos p0. 



1 The calculation is at the same point as in Bergemann, eq 9, p. 655, 
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We substitute this in the area integral: 

A=- J d<p ( kl + 2kook tlu cos/i0 cos(vk z — vkqq sin 0tan#) J . 

The cosine of a sine in an integral leads to Bessel functions, where 

i r 

JuM — — / cos(fix — \ sm(x))dx. 

We transform our integral in such a form. Expanding and regrouping, using trigono- 
metric identities, we have 

cos cos(y/^ — Koo sin tan 9) 
= -( cos(ii0 + vk z — ukqq tan 9 sin 0) + cos(/i0 — z//t 2 + ^k 00 tan 9 sin <; 

- ( cos(ii0 — z//t 00 tan 9 sin 0) cos z/k z — sin(/i0 — z//t o tan 6 1 sin 0) sin z//t 2 



(B-3) 



2 

+ cos(/i0 + z/kqo tan 6 1 sin 0) cos vk z + sin(/i0 + vkqq tan # sin 0) sin uk 



The integral of the second and last terms will be zero, and we are left with two 
Bessel functions: 

dtp cos(/i0 =F z/kqo tan 9 sin 0) cos z//t 2 = /^(ii/^oo tan 9) cos z//« 2 

The Bessel functions are either symmetric or antisymmetric, depending whether /x is 
respectively even or odd. We rewrite the sum of the integrals 

^t(^Kootan^) + J At (— z//too tan #) = 2 J M (z/ k o tan 0), /j even, 

= /i odd, 

Finally, the area is 

A = — ^ + 27r 00 ^ J„(z/ /t 00 tan 0) cos z//t z , u et>en 
cosb 1 ' COSb 1 

This is the result when imposing an orientation angle 0o of zero. When using 0o ^ 0, 
the calculation is more complicated but essentially the same and yields the result of 
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Bergemann et al., eq. 13, p. 657 

— ^ + 27r 00 ^ Jn(vK 00 tanfl) cos/i0o(— lW 2 cos i/n 2 , even u, 
cos# cosfe 1 

lMi + 27 ry ^^J M (^ otan^)cos/i</»o(-l) { ^ 1)/2 sinz/^, odd fi, 
cos 9 / — ' cos 9 

(B.4) 

With this equation, adding all the different A^ u to the area will give the total area 
of the intersection of a plane tilted by 9 and the Fermi surface with a shape defined 
by the k^ u . We cannot go any further unless we choose one k^ v and set all the others 
to zero, except /coo- 

We are finally required to calculate dHvA oscillations from a single warping para- 



meter. We replace the last equation with 0o — into the equation IB.lt an d obtain 

M oc f dk z sin ( — ^- + 2n'S~ y 00 ^ J»(z//t 00 tanfl) cosvk z ] 

(i/c^ sin (a + b^ u cos uk z ) 

dk z sin a cos(6^ cos i/ac^) + cos a sin(6^ cos uk z ), 



where the second term is zero by symmetry. We have once more an integral of the 
form of a Bessel function, of order zero. Putting back the symbols into a and b, we 
obtain (with \i even): 



hX 7T^ \ 1 fluhX kookf, 
e cos 9 J v °\ e cos 9 
2tiF X\ 1 f2nAF^ u X 



M oc sin ( — ^ ]— J ( — ^-J^(^K 00 tani 



= sin — -J — J^vkqq tan 6) , 

\ cos# / z/ \ cost' / 

where F is the average dHvA frequency, and AF^ is the difference which produces 
a beat pattern. 
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Fourier series of the functions 

sin (AsinfW)) and cos (Asin(cjt)) 

We begin with the function cos (Asin(wt)). The Fourier series requires to calculate 

oo 

cos (A sin(u;t)) = ao + a k cos(Axjt) + b k sin(Axo>£), 



with 



1 r 

a k = — cos (A sin(a;t)) cos(kut)d(tut), 

^ J -IT 

i r 

bk = — / cos (Asin(u;i)) sm(kut)d(ut), 

i r 

ao = — J cos (\sm(u>i))d(u>t). 

Using trigonometric properties, the first two equations can be rewritten 

a*; = — J ^ cos (kut + A sin(wt)) + cos (hut — A sin(u;i))^ d(uit), 
b k = 0, 

the second being zero by symmetry. The equation for a k corresponds to Bessel's first 
integral 

i r 

<^fc(A) = — / cos(kx — Xsm(x))dx, 

271_ J-7T 

and thus 

Ofc = Jfe(A) + Jfc(-A). 
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The constant offset ao also possesses the form of Bessel's first integral, with k = 0: 

ao = ^o(A). 

A similar calculation can be done with the function sin (Asin(u;£)). In this case, a k 
and ao equal zero, by symmetry, and the remaining term is 

b k = — J ^ cos {kut + A sin(W)) — cos (kut — A sin(o;t))J d(ut), 
= J k (X)-J k (X). 

The Bessel functions are either symmetric or antisymmetric, depending whether k is 
respectively even or odd. One can rewrite these results as 

Jk{X) + Jfc(-A) = 2J k (X), k even, 
= k odd, 



and 



J k (X)-J k (-X) = 2J k (X), k odd, 
= k even, 

The Fourier series are then: 

oo 

cos (A sin(wt)) = Jo(A) + 2 J 2k cos(2kut) 

k =i 

oo 

sin (A sin(wt)) = 2 J-zk-x cos ((2fc — \)uj£) 

k =i 

Finally, one may approximate J k (X) for A < 1. The Bessel function is defined as 

7(X) - V ^ ( X \ 
Jk ^> ~ 2^ n u n + k + 1 y\ 2 ) 

n=0 



2n+ k 



2/ ^ n\(n + k)\\2 

7 n=0 v 



When A is 1 or smaller, the terms of the sum become very rapidly small with increasing 
n, and keeping only the first one, n = 0, which is equal to l/k\, one obtains 

\ k 

J k (X)^ 1 



2 k k\' 
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Cooling power of the adiabatic 
demagnetisation refrigerator 



The principle of adiabatic demagnetisation refrigeration uses the entropy of an isol- 
ated paramagnetic material under a magnetic field. If the spins have a total angular 
momentum of J = L + S, their interaction with the field is e = —g(J, L, S)fi B 3 ■ H. 
From this, one can calculate the free energy, and derive the magnetisation and the 
entropy 0: 



M 
S 



OF 

— ) = Ngfi B JBj(gJx) 



x 



k B T 



dF_\ 
df) 



- I — - I = Xk 



B 



sinh g(J + | 



x 



1U ^nhf *— ~ 9J*Bj(gJx) 



k B In TV!, 



(D.l) 



with Bj the Brillouin function, 



Bj(y) = (l + jj)coth(l + ±j)y 



1 ^ y 

2J COth 2J- 



Both the entropy and the magnetisation are functions of H and T through x = 
only. Consequently, for an adiabatic process where H varies, T evolves such that 
x remains constant. In such a scenario, reducing H to zero should produce zero 
temperature. It is not the case, since below a certain field, H , the intrinsic internal 
field of the paramagnetic material prevents reaching zero field, from very small but 
non- vanishing spin interactions. When the system reaches the minimum field of Hq, 
the cooling power vanishes and one obtains a base temperature T B t- 

This analysis works when only the spin system is held in adiabatic conditions 
in the cryostat and no significant load of heat capacity is thermally connected. In 



L Ashcroft and Mermin [12 , chapter 3 



168 



Appendix^ 



169 



the ADR system, the heat capacity of the LTS is not negligible. As temperature 
is reduced, heat is transferred from the LTS to the paramagnetic material, and the 
entropy of the spin system S is not constant. As one reduces the magnetic field by 
AH, the temperature reduces by AT, but their ratio also changes by Ax, as does the 
entropy of the spin system by AS, and 



where C^f is the heat capacity load of the LTS. The LTS is made of non- magnetic 
metal: 



where N is a number of moles. One takes a path during demagnetisation where the 
rate of temperature decrease is constant. Integrating eq. ID. 21 over this path: 



Figure ID. II shows curves of T versus H when appropriate constants are used, cal- 
culated by equating the entropy values. It consists of a contour plot of the total 
conserved entropy of the paramagnetic system and LTS. One observes that the rate 
of decrease of the LTS temperature accelerates during cooldown with a constant rate 
of demagnetisation, which corresponds to the behaviour of the ADR, and is due to 
the decrease in heat capacity of the LTS with T. 



- TAS = C 



(T)AT, 



(D.2) 



C L J S {T) = N{ 1 T + I3T^), 
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Magnetic field ( T ) 

Figure D.l: Contour lines of the total entropy of the LTS of the ADR while in the 
adiabatic state. The red lines correspond to high entropy, which decreases towards 
the blue curves. The dashed line represents the intrinsic ordering magnetic field of 
the paramagnetic material, at which the cooling process stops. 



Appendix E 

Raw data from the rotation study 



This appendix presents raw dHvA data for the reader to appreciate and get a feeling 
for. The complexity of the data can be seen just by looking at the patterns formed 
by the raw data, before even performing any Fourier transforms. The data which is 
presented here is only the rotation study of dHvA measured in Cambridge, which was 
performed at a base temperature 30 ± 5 mK, during around 20 days. 

The angle was measured using a rotating knob potentiometer, which was turned 
with the mechanism by a computer controlled electrical motor. The resistance is 
then simply proportional to the angle, and was measured in kll The conversion from 
resistance to angle values was later performed by simply taking photos of the coil 
system at various potentiometer values, method that turned out to be the simplest 
We obtained conversion values of 8.23° /kfl x R — 8.23° for coil A, which held sample 
C698I, and 8.45°//cf2 x R— 17.8° for the other one, determined with a linear regression 
of angles measured on photos clS cL function of _d resistances! The rotation 
was performed in steps of 0.2 kf2, which corresponds to approximately 1.6°. Moreover, 
it was found during the experiments that the first few data sets had rotation steps of 
less than that value, probably due to torsion building within the rotation mechanism 
before the rotation actually started at the sample stage. 

The graphs presented in this section show data at resistance values of 0.6 to 7.6kfi, 
although the measurements started at zero resistance, and this is due to the fact that 
the first four were almost identical. The angles extend between -3.3° to 54.3° for 
C698I, and from -12.7° to 46.4° for sample C698A. The difference in those angle 
intervals, for the same position in rotator, is due to the an error in the installation of 
the coil onto the bobbin system (see figure l2~2| section 12. 1 .3j) . and the error in cutting 
the sample with respect to their c-axis. 

: I thank Michael Sutherland and Swee K. Goh for this careful measurement, which was done in 
my absence. 

2 The conversion numbers come from a series of photos, but the zeros are not from the linear 
regression, but rather from the symmetry of the data, i.e. the F = F / cos 9 law. 
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Figure E.l: Raw data from the rotation study of sample C698I. This first half of 
the measurements were performed at potentiometer values between 0.6 and 4 kQ, 
corresponding to angles of -3.3° to 24.7°. 




Figure E.2: Raw data from the rotation study of sample C698I. The second half of 
the measurements were performed at potentiometer values between 4.2 and 7.6 kQ, 
corresponding to angles of 26.3° to 54.3°. 




Figure E.3: Raw data from the rotation study of sample C698A. This first half of 
the measurements were performed at potentiometer values between 0.6 and 4 kQ, 
corresponding to angles of -12.7° to 16.0°. 




Figure E.4: Raw data from the rotation study of sample C698A. The second half of 
the measurements were performed at potentiometer values between 4.2 and 7.6 kQ, 
corresponding to angles of 17.7° to 46.4°. 
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Numerical self-consistency 
verifications for calculations of the 
quasiparticle mass 

We present in this appendix simulations of dHvA data sets along with their analysis 
for the quasiparticle mass. We have performed two type of calculations, presented in 
two separate sections. The first type was done with a simulated mass enhancement, 
using the same calculation procedure as was used in section 13.41 with two parameter 
fits of the LK relation, in order to verify the validity of the method. The second 
was done with a constant mass using a three parameter fit, in order to demonstrate 
how systematic errors can arise when the signal to noise ratio is poor and a free 
offset is used in the LK fit. Gaussian noise of various amplitudes was deliberately 
added to the signal, such that the simulations imitate real life experiments. We 
observed that the method we used to analyse the dHvA data from the Cambridge 
experiments produced the appropriate results, and that spurious mass enhancements 
can be produced when using the same method as that which was used in the work of 
Borzi et al. [40J, confirming our discussion of section 14.4.11 



F.0.2 Calculation of quasiparticle mass enhancements with 
two parameters fits 



The top part of figure IF. II presents simulated dHvA oscillations, produced by evalu- 
ating the following function, 



M(X) oc LK(X,T,m{X))D(£,F,X) 



cos(2ttFiX) + cos(27rF 2 X) 
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Magnetic Field ( T ) 

Figure F.l: Simulations of dHvA data incorporating a strong mass enhancement 
around 7.9 T, and quasiparticle mass values extracted from the simulations using the 
normal method. Gaussian white noise was added to the second and third simulations, 
of standard deviation 10 and 20 mVrms. The data shown were calculated with a 
temperature of 90 mK. 
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where LK(X) is the LK function, D(X) the Dingle factor (see section [2.2. ip and f{X) 
represents white noise with a Gaussian distribution. The data was calculated for a 
field range of 5 to 15 T, using two slightly different frequencies, Fx = 4.20/cT and 
F 2 = A.25kT, in order to produce ordinary interference patterns, a mean free path of 
I = I/mi and a field dependent quasiparticle mass m* = m*(X). The function was 
evaluated at ten temperatures between 90 and 500 mK Q. The mass was defined as 
constant and equal to 12 m e away from 7.9 T and increases up to 60 m e following 
a Lorentzian function, centred at 7.9 T, of half-width at half-height of 0.3 T. No 
particular reason lies behind this choice of this function, the only requirement being 
its reproduction when calculating the mass using the simulations. One may observe 
that at 7.9 T, the oscillations amplitude is very close to zero. 

The second and third data sets presented in figure IF. II feature noise, which was 
obtained using a random number generator, and subsequently weighed such that its 
probability distribution was Gaussian. The noise amplitudes quoted in the graphs, 10 
and 20 mVrms, refer to the standard deviation a, even though the visible amplitude 
appears larger. The noise in real experiments may not be necessarily of that type. 
However, from lack of better knowledge over its distribution and its correlation and 
due to time constraints, we used the simplest realistic type of noise. 

The method used to analyse the simulated dHvA data was exactly the same as 
that employed for the extraction of the quasiparticle masses presented in section 
13.41 figures 13.161 and 13.171 which we describe again here. The data was processed 
using 100 sections of length of 0.01 T _1 , which were allowed to overlap, their centres 
evenly spaced in inverse field. Fourier transforms were performed on these for all 
temperatures, and the peak in the power spectrum corresponding to the oscillations 
was integrated. The result was fitted using the normal two parameter LK function, 
eq. 12.71 (see section [2.2. 41) . Each fit was inspected by eye, and the results for which 
the data distribution was too noisy to obtain a reasonable fit were discarded. 

The lower part of figure IF. II shows the results of the mass analysis for the three 
situations described, with two different levels of noise and without noise at all, and 
the error bars correspond to the fit error given by the least squares method. In red 
is plotted the field dependent mass m*(B) which was used in the creation of the 
oscillations. Without noise, the calculation procedure reproduces the original mass 
properly except where its value is very high (note that error bars are much smaller 
than the size of the symbols). This is partly due to the width of the data sets used 
in the Fourier transform, which leads to a convolution of the mass function with the 
Fourier transform window of width of about half the window size, 0.01 T _1 . But 

lr These were the same as those in the real experiments carried out in Cambridge: T = 90, 120, 
150, 180, 220, 280, 330, 380, 430 and 500 mK. 
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also, it is due to systematic errors related to the lowest temperature point, which was 
deliberately chosen high, 90 mK, in order to reflect the real data. Very high masses 
lead to temperature distributions in which only the lowest data point possesses a non- 
zero amplitude, leading to errors much larger than those given by the fitting algorithm. 
Nevertheless, these results indicate that without noise, any mass enhancement should 
be properly detected by our data analysis procedure. 

With 30 mVrms of noise added to the oscillations, we performed the same analysis 
and obtained a similar mass enhancement, but cut off near 26 electron masses. The 
gap in the data reflects the relative size of the oscillatory amplitude and the noise in 
the region where the mass is very high, and corresponds to inverse field windows where 
the LK fit was not possible at all. In the low field side of the mass enhancement, it 
is recovered, but the results are noisier due to the exponential decrease of the Dingle 
factor. Effectively, large underestimates arise at nodes of the beat pattern, where the 
signal to noise vanishes. We conclude however that with such a signal to noise ratio, 
a mass enhancement can still be detected. 

Finally, using a noise amplitude of 60 mVrms, we performed the same analysis 
again, and different results were obtained. In the enhanced mass region, LK fits were 
much poorer and systematic errors were observed, especially in the low field side. The 
mass was mostly underestimated at every node of the beat pattern, but was accurate 
otherwise. However, the mass enhancement was cut off at around 18 electron masses, 
reflecting how the noise amplitude may hinder its detection. 

Quantitatively, we found that it is when the signal to noise ratio of the 90 mK 
temperature point decreases below 0.5 that the mass analysis becomes difficult. We 
then compared this with the signal to noise ratios for all the frequencies of which 
masses were extracted and presented in figures I3.16l and r3.17l for the c-axis direction. 
For the low field side frequencies at 4.2, 1.8, 0.9, 0.43 and 0.15 kT, the lowest signal 
to noise ratios recorded very near the transition (away from any beat pattern nodes) 
were respectively of 8, 20, 1.8 1.4 and 0.8. In the high field side, we presented 
data only for the 4.0 and 1.6 kT peaks, and they possessed signal to noise ratios of 
respectively 4 and 12. Consequently, none of the frequencies analysed were in the 
situation depicted here, where the noise hinders the mass analysis. We conclude that 
if a mass enhancement had been present for any of the pockets we analysed, we should 
have detected it. 

F.0.3 Systematic errors with three parameters fits 

We discussed in section 14.4.11 how the use of three parameter fits in order to accom- 
modate data of poor signal to noise ratio could lead to spurious continuous mass 
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12.5 - No noise 3 parameter fits 




60 ■ ° Noise 150 mVrms 3 parameter fits 

o 




Q| 1 1 1 1 1— 

4 6 8 10 12 14 

Average Field ( T ) 



Figure F.2: Simulations of dHvA data using a constant quasiparticle mass, without 
noise in the first case, and with different noise levels for the other two. Quasiparticle 
mass values extracted from the simulations using three parameter fits produce strong 
spurious mass enhancements where the signal to noise ratio is low. The data shown 
were calculated with a temperature of 90 mK. 
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enhancements as the signal to noise ratio evolve with magnetic field. Such a phe- 
nomenon is relatively easy to simulate, and we present such a calculation in this 
section, giving further support to our suggestion that it is the phenomenon which led 
to the mass enhancement measured by Borzi et al. |4Uj . 

Simulations similar to those of the previous section were performed, without an 
enhancement of the quasiparticle mass, which was defined as constant and equal to 
12 m e . All other parameters were the same, except for higher noise levels, which were 
of 60 and 150 mVrms. The quasiparticle mass analysis was performed in exactly the 
same way as previously, with the exception that the non-linear LK fits were performed 
using three parameters, as defined in section 12.81 and discussed in section 14.4.11 

The top three graphs in figure IF.2I present the simulated data, the first featuring 
data without noise, while the other two possess noise levels of 60 and 150 mVrms. 
The bottom three graphs present the mass extracted with the three parameter LK 
fits. In the first case, without noise, we observed that the quasiparticle mass oscillated 
around the original value of 12, deviating at each beat pattern node. As noise was 
added to the oscillations, we found that these deviations evolved into strong mass 
enhancements as the signal amplitude became comparable with the noise level. With 
a noise level of 150 mVrms, stronger mass enhancements were calculated, of up to 60 
at the lowest signal to noise ratio. Note that again here, all the error bars that we 
present correspond to the fit errors given by the algorithm. These seem inconsistent 
in places, being small where strong enhancements were observed. We argue that this 
points strongly to how misleading three parameter fits can be: very good fits can be 
obtained where noise is the dominant component of the data. 

The mass enhancements appear at beat pattern nodes, where the signal to noise 
ratio vanishes, as we discussed in section 14.4.11 However, here, one can observe that 
the increase is more and more pronounced the lower the signal to noise ratio is away 
from beat pattern nodes, and that the system becomes unstable to small changes in 
signal amplitude. From the mass extracted from data with 60 mVrms of noise, we 
determined that it is when the signal to noise ratio away from beat pattern nodes 
is of around 3 or less that the system becomes susceptible to large mass changes. 
This essentially means that for a system where the signal to noise ratio is not good 
where the signal amplitude is at its maximum, strong mass enhancements can be 
expected if one uses three parameter fits in regions where beat pattern nodes arise. 
It is possible that in such a case, two parameter fits are not possible either, due to 
the bad quality of the data. Such data should then simply not be used for a field 
resolved mass analysis. 

We believe that it is the assignation of an additive constant to a random value 
added to the LK function that lies behind this phenomenon. Although in many 
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fits, this seems very appropriate, and small error bars are obtained, the information 
extracted can be very inaccurate. We conclude by suggesting that only LK fits with 
two parameters should be used, and only with data that possesses a signal to noise 
ratio higher than 0.5. 
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